Я запускаю пятнадцать ppml
регрессий, используя for loop
(код ниже): (d
- список моей зависимой переменной.) ppml
из пакета gravity
.for (я в 1:15) {
fit <- ppml(
dependent_variable = d[i],
distance = "y",
additional_regressors = c("x1", "x2" , "x3"),
robust = FALSE,
data = mydata
)
print(d[i])
print(fit$coefficients[1:7])
#print(summary(fit))
}
Вывод ppml странный, поскольку это не просто один класс:
> class(fit)
[1] "glm" "lm" "ppml"
Я хочу создать 15 фреймов данных исохранить каждый из результатов в одном кадре данных.
Есть предложения?