0

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 を使用したグループごとの処理にうるさいのはなぜですか?

4

1 に答える 1

0

エラーが一致する引数であるため、それが機能したlm glmかどうかは疑問です。polr

全体をラップする必要がありますlist

dt2[,list(dtmsm=list(msmMC(y1, x1, swNx1, Nx1))),by="h"]

または、data.frame、list、または data.table の可能性があるオブジェクトを返すように見えるlist呼び出しを誤って配置しただけかもしれません。msmMC

dt2[,list(dtmsm=msmMC(y1, x1, swNx1, Nx1)),by="h"]
于 2013-03-12T00:50:45.207 に答える