我有一个让我发疯的问题。我写了一段代码,它可以在我的计算机上运行(以及其他多个),但是当我的主管在他的计算机上运行它时,除了drop1()命令外,一切都正常。
我的数据和代码如下:
Data frame structure (condensed)
> str(Df)
'data.frame': 96 obs. of 13 variables:
$ BWA : Factor w/ 3 levels "L","M","H": 3 3 3 3 1 1 1 1 2 2 ...
$ SBW : int 25 50 75 100 25 50 75 100 25 50 ...
$ LS : num -0.16 0.137 0.028 0.044 -0.26 ...
$ Def : num 0.1316 0.4304 0.0338 0.1668 -0.1398 ...
$ fSite : Factor w/ 8 levels "1","2","3","4",..: 1 1 1 1 1 1 1 1 1 1 ...
M1<-gls(LS ~ 1+SBW+BWA,method='ML')
M2<-lme(LS~ 1 + SBW+BWA, random= ~1|fSite,method='ML')
drop1(M1,test='Chisq')
drop1(M2,test='Chisq')
我得到了
drop1(M1,test ='Chisq')
>Single term deletions
>Model:
> LS ~ 1 + SBW + BWA
> Df AIC LRT Pr(>Chi)
> <none> -156.00
> SBW 1 -151.48 6.5166 0.010687 *
> BWA 2 -148.43 11.5711 0.003072 **
---
> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(drop1(M2)的类似输出)
然而,我的主管得到了:drop1(M1,test ='Chisq')
> Error in terms.default(object) : no terms component nor attribute
要解决这个问题,我已经尝试使用各种版本的R的多台计算机(Mac和PC)上的代码文件。每次drop1()命令工作正常。所以我知道这是一台特定的计算机的问题,而不是我的代码的问题。
在麻烦的计算机上,我厌倦了解决问题:
无论我做什么,我都无法让drop1()工作。
有没有人知道为什么像drop1()这样的基本命令在这台计算机上无效。如果有人能在这里看到我遗漏的内容或提出另一种方法来解决这个问题,请告诉我。