Differences

This shows you the differences between two versions of the page.

--- econometricsr [2025/06/02 13:07] – [Regressions] glebelg
+++ econometricsr [2025/06/04 11:44] (current) – glebelg
@@ Line 1: / Line 1: @@
-====Summary Stats====
+=====Start a script =====
   library(stargazer)
+  rm(list = ls())
+set working directory
+  setwd("/path/to/my/directory")
+data from excel
+  dataset <- readXL("data.xlsx", rownames=FALSE, header=TRUE, na="", sheet="Sheet1", stringsAsFactors=TRUE)
+data from Rdata
+=====Summary Stats=====
   stargazer(dataset, type = "text", title= "Summary Statistics", digits = 2, out="output.csv")
@@ Line 9: / Line 22: @@
-====Regressions====
+=====Regressions=====
+linear
   reg1 <- lm(depvar ~ var1 + var2 + var3, data=dataset)
+use rlm for robust regressions
+probit
   reg2 <- glm(dummyvar ~ var1 + var2 + var3, family=binomial(link=probit), data = dataset)
+table
   stargazer(reg1, reg2, header=FALSE, title="Regression Table", type='html', digits=2 , out="output.html")
+=====Exemples=====
+  install.packages("MASS")
+  install.packages("readODS")
+  install.packages("tidyverse")
+  install.packages("devtools")
+  devtools::install_github("gvelasq/tidytab")
+#load libraries
+  #import ods
+  library(readODS)
+  #import xls/xlsx
+  library(readxl)
+  #robust lm
+  library(MASS)
+  #for tables
+  library(stargazer)
+  library(tidytab)
+  # for data manipulation
+  library(tidyverse)
+# set work directory
+  setwd("/home/glebelg/nextcloud/Work/time to code/")
+# import ods
+  db1<-read_ods(path = "dbtest.ods", sheet = 1)
+  db2<-read_excel("dbtest.xlsx")
+#generate variables in dataset
+  db1$niveauscore <- db1$niveau * db1$score
+#option2 - si beaucoup d'opérations!
+    attach(db2)
+    db2$niveauscore <- niveau * score
+    detach(db2)
+    attach(db1)
+    db1$niveauminscore <- niveau - score
+    detach(db1)
+# OLS regressions
+    model1 <- rlm(score ~ age, data=db1)
+    model2 <- rlm(score ~ niveau, data=db1)
+    model3 <- rlm(score ~ age + niveau, data=db1)
+    stargazer(model1, model2, model3, type="text",
+              title="Regression Results",p.auto=TRUE, single.row=FALSE,
+              ci=FALSE, ci.level=0.9, omit.stat=c("ser"))
+#summary statistics table
+  stargazer(as.data.frame(db1), type="text", summary=TRUE)
+# tabulate
+  table ("Age"=db1$age,"Sexe"= db1$sexe)
+  db1 %>% tab(age, sexe)
+  db1 %>% tab1(age, sexe)
+  db1 %>% ftab( sexe, age)
+=====Import/Export=====
+import excel
+  library(readxl)
+  db <- read_excel("db.xlsx")
+  View(db)
+import stata
+  library(haven)
+  db <- read_dta("db.dta")
+  View(db)
+export excel
+  install.packages("xlsx")
+  library(xlsx)
+  write.xlsx(db, file = "db.xlsx")
+=====Data Manipulation=====
+====drop publicates====
+  newdata <- distinct(dataset,varname, .keep_all = TRUE)