使用从循环(R)中选择的变量的Ancova

时间:2014-06-14 19:50:21

标签: r loops static-analysis mixed-models

我在Ancova中使用了1,000个因变量(y)。所有模型中的独立变量(x)都是相同的。

Yvar1, Yvar2, …, Yvar1000, x1, x2
1       2          5       11  16
2       3          6       18  23

我需要一个R代​​码来替换每个Yvar的example$Y.to.use值。

我试过了:

for (i in 1:1000){ 
example$Y.to.use<- example$paste("Yvar",i, sep="")

# Error: attempt to apply non-function

paste("fitted.model",i)<-lme(log(Y.to.use+1)  ~ x1*x2, data=example, random=~1| Exp/Person)
} 

# Error in paste("fitted.model", i, sep = "") <- lme(log(Y.to.use+1)  ~ x1*x2  : target of assignment expands to non-language object

我将创建一个表格,其中包含每个Yvar的拟合模型的系数。

0 个答案:

没有答案