data.table に次のデータがあります。
h x1 y1 swNx11
1: 1 39.075565717 0 1.03317231703408
2: 1 40.445951251 0 7.14418755725832
3: 1 37.800722944 0 0.435946586361557
4: 1 41.085221504 0 0.381347141150498
5: 1 36.318077491 0 0.497077163135359
---
24996: 25 39.110138193 0 0.942922612158002
24997: 25 39.331940413 0 1.42227399208458
24998: 25 37.479473784 0 0.390657876415799
24999: 25 35.892044242 0 0.599937357458247
25000: 25 40.699588303 0 0.486486760245521
svyglm でそれらを分析する関数を作成しました。
msmMC <- function(y, x, sw, name){
msm <- svyglm(y ~ x,family=quasibinomial(link="logit"),design = svydesign(~ 1, weights = ~ sw))
out <- cbind("name",coef(summary(msm))[2,1],coef(summary(msm))[2,2])
return(out)
}
msmswNx1<-dt2[,list(dtmsm=list(msmMC(y1, x1, swNx1, Nx1))),by="h"]
outNx1 <- unlist(dt.lm[,msmswNx1])
この関数を実行すると、次のエラーが発生します。
[.data.table(dt2, , list(dtmsm = list(msmMC(y1, x1, swNx1, : 'by' または 'keyby' の列または式 1 は型リストです。列名を引用しないでください。使用法: DT[,sum(colC),by=list(colA,month(colB))]
ただし、glm や pollr などの別のモデルでも問題なく動作します。それで、ここで何が起こっているのですか?svyglm が data.table を使用したグループごとの処理にうるさいのはなぜですか?