Question

我是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”必须为数字

有人可以帮我解决这个问题吗？

Answer 1

我不确定，但是也许应该用data $ x2替换x2，然后将b转换为这样的数字值：

    data=request.data
    j_data = json.loads(data)

    factory = {
        'id': factories[-1]['id'] + 1,
        'fac_name': j_data['fac_name'],
        'lat': j_data['lat'],
        'lng': j_data['lng']
    }

如何在进行置换测试时解决错误

1 个答案: