Я новичок в R. Я хотел бы использовать тест перестановки, чтобы увидеть связь между y и x1.В этом случае я хочу, чтобы y был фиксированным, тогда как x1 переставлен, а другой регрессор, x2, является постоянным.Поскольку у меня есть продольный набор данных, а моя зависимая переменная y - это данные подсчета, я использую фиксированный эффект
. Я запускаю следующий код
setwd("D:/trial")
getwd()
trial <- read.csv("trial.csv")
b <- c()
r <- c()
data = trial[,c(-1,-6)]
for(i in 1:500){
data$x1 <- sample(data$x1)
placebo <- glm(y ~ data$x1 + x2+ as.factor(code) + as.factor(year),
data=data, family=poisson)
b <- append(b, placebo$coef[4])
r <- append(r,cor(data$x1, data$y))
coef(summary(placebo))
}
hist(b)
и получаю сообщение об ошибке:
Ошибка в eval (predvars, data, env): объект 'x2' не найден hist (b) Ошибка в hist.default (b): 'x' должно быть числовым
Может кто-нибудь помочь мне решить эту проблему?