我有一個 for 回圈，我只想改變一個引數p，同時保持另一個引數q不變。p并且q是 . 中的“自動回歸”和“移動平均”訂單gls()。如何做到這一點？這是我的代碼（不作業）：
cor.results <- NULL
for(p in 0:3) {
cor.temp <- gls(di.posaf~di.index factor(subject):day factor(subject), data=PMBCall2, method='ML', correlation=corARMA(p=p, q=0, form = ~ day | subject)) #I want to vary p only and keep q constant
cor.results <- rbind(cor.results, c(p, q, logLik(cor.temp), AIC(cor.temp)))
}
uj5u.com熱心網友回復：
事實證明，corARMA()不接受p=0 & q=0 同時，所以我將起始值p增加到 1
cor.results <- NULL
for(p in 1:3) {
cor.temp <- gls(di.posaf~di.index factor(subject):day factor(subject), data=PMBCall2, method='ML', correlation=corARMA(p=p, q=0, form = ~ day | subject)) #p should start from 1
cor.results <- rbind(cor.results, c(p, q, logLik(cor.temp), AIC(cor.temp)))
}
轉載請註明出處，本文鏈接：https://www.uj5u.com/qukuanlian/411464.html
標籤：
上一篇：如何遍歷二進制數字的嵌套串列并將其十進制數字附加到串列中？例如：[[0,0,1,0,1],[1,0,0,0,1]]==[5,17]