Question

我使用r-package nlme的函数lme估计了两个截距的混合多级模型。

然后，我使用plot（ACF）函数通过肉眼检查来检查自相关。

我第一次指定maxlag = 16进行绘图。

现在我有两个问题：首先，maxlag参数似乎被卡住了，即，即使将maxlag设置为其他值，也都以maxlag = 16绘制了更多图。 2.即使滞后0的值明显为1，也会在y = 0.8处裁剪图。

下面，我分享各自的观点，以期希望得到有关如何解决这两个问题的答案或投入。

Link to the dataset and if prefered to copy-paste to the following code-script as well:

#read.dataset:

datafclr <-read.csv("datafclr.csv", header = TRUE, sep = ",", dec = ".", fill = TRUE)

#required packages:

library("Matrix")
library("nlme")

#model-estimation:

tim2 <- lme(fixed=EERTmn ~ male + female + 
              (male:time7c)          + (female:time7c)          +                           
              (male:IERT_Cp)         + (female:IERT_Cp)         + 
              (male:IERT_Cp_Partner) + (female:IERT_Cp_Partner)-1,

            control=list(maxIter=100000), data=datafclr,                                             

            random=~male + female -1|dyade/female, correlation=corAR1(), na.action=na.omit)

summary(tim2)

#checking for autocorrelation:

plot(ACF(tim2, maxlag = 16), alpha = 0.01)

结果如下图所示：

This results in thin plot

当我更改最大时差：

plot(ACF(tim2, maxlag = 10), alpha = 0.01)

It results in the same plot

非常感谢！

最好，帕特里克

Answer 1

Joes Schwartz帮助我解决了R-Studio社区中的这些问题。对于这种情况，有人会遇到与我相同的困难，我在这里分享他的答案：

第一个问题：需要将maxlag键入maxLag，并且该函数可以正常运行。

第二个问题：以下链接下的详细帮助：

https://community.rstudio.com/t/resetting-plotting-settings-plot-acf-data/19441

如何在R-studio的plot（ACF）设备中更改或重置参数

1 个答案: