Фрейм данных, состоящий из символов, преобразуется в факторы после команды pdata.frame в R. Как заставить работать OLS в пуле?

Question

Я использую набор данных панели:

• y - моя независимая переменная, равная 0 или 1 -> цифра
• x1 мои индивидуумы -> цифра
• x2 - мои индикаторы времени -> числовые
• x3, x4, ..., x65 - мои независимые переменные -> символ

В приведенном ниже коде Iпреобразовать все переменные в символы и сообщить R, что я использую данные панели с помощью команды pdata.frame в последней строке.Однако теперь проблема заключается в том, что команда pdata.frame преобразует переменные x2 и x3 (отдельные лица и индикатор времени) в факторы, даже когда stringsAsFactors = FALSE.

#Regressions
df=read_excel("C:/Users/Luuk/Desktop/Master Thesis EME/Data/indep_dep_indlevel.xlsx")
df_dep=data.frame(df[,79])
count=as.data.frame(rep(1:3669, times=1, each=3))
df=cbind(count,df[,3:79])
df_indep=data.frame(df[,c(1:5,8,10:15,17:25,27:44,45,53:77)])
dflm=cbind(df_dep,df_indep)
dflm1 <- data.frame(lapply(dflm, as.character), stringsAsFactors=FALSE)

names(dflm1)[c(2:66)] <- c(paste("x", 1:65, sep=""))
names(dflm1)[1] <- "y"
dflm2=pdata.frame(dflm1,index=c("x1","x2"),stringsAsFactors=FALSE)

Следовательно, следующая оценка модели объединенного OLS дает ошибку:

Ошибка в классе (x) <- setdiff (class (x), "pseries"): <br>добавление класса "factor" к недопустимому объекту. Дополнительно: Предупреждение: в model.response (mf, "numeric"): использование type = "numeric" с ответом фактора будет проигнорировано

xnam <- paste("x", 3:65, sep="")
Formula <- formula(paste("y ~ ", paste(xnam, collapse=" + ")))
fit=plm(Formula, data=dflm2,model="pooling")

Howя могу заставить мою процедуру оценки OLS в пуле работать?

Answer 1

Запуск модели панели с возможностью объединения в пул эквивалентен простому запуску простого OLS lm (y ~ x). Я не понимаю, зачем вам делать больше, если вы пытаетесь получить объединенные оценки. Я не могу воспроизвести вашу ошибку, используя шаги, аналогичные вашим (см. Мой код), используя файл xlsx, который я создал для этой цели. Пожалуйста, включите минимальный рабочий пример, который воспроизводит вашу ошибку.

df=read_excel("~/Downloads/strtest.xlsx")
df_dep=data.frame(df[,1])
df2=data.frame(lapply(df, as.character), stringsAsFactors=FALSE)
names(df2)[c(1:3)] <- c(paste("x", 1:3, sep=""))
df2=cbind.data.frame(c(1,2,3,80),df2)
names(df2)[1]='y'
df3=pdata.frame(df2,index=c("x1","x2"),stringsAsFactors=FALSE)
plm(y~x1+x2+x3, data=df3,model="pooling")