R: Error in eval(predvars, Daten, env) : numeric 'envir' argument nicht von der Länge eines inneren.pars() in multgee Paket

Habe ich einen Daten-frame dfmit den folgenden Daten:

family_id_cass   twin_order  ls75cat
  10026          1           2
  10026          2           2
  10071          1           3
  10071          2           3
  10132          1           1
  10132          2           2

Ich würde gerne intrinsic.pars(), was ich so mache, ist die folgende:

intrinsic.pars(df$ls75cat, df$family_id_cass, df$twin_order, rscale="ordinal") aber ich bekomme die folgende Fehlermeldung:

Error in eval(predvars, data, env) :
numeric 'envir' arg not of length one

Was mache ich falsch?

  • Paket multgee
InformationsquelleAutor roland | 2015-02-21
  1. 1

    intrinsic.pars() zu sein scheint, geschrieben mit einem sehr temperamentvollen Schnittstelle, die effektiv ignoriert die meisten Parameter. (Ich vermute das gemeint war, zu einer internen Funktion.) In jedem Fall müssen Sie, um eine Daten.Rahmen, um das argument 'Daten', die die Variablen enthält y,id, und time:

    df <- read.csv(text=
"family_id_cass,twin_order,ls75cat
10026,1,2
10026,2,2
10071,1,3
10071,2,3
10132,1,1
10132,2,2")

newData<-data.frame(y=df$ls75cat,
                    id=df$twin_order,
                    time=df$family_id_cass)

# either of these calls will do:
intrinsic.pars(y,newData,id,time,rscale="ordinal")
intrinsic.pars(data=newData,rscale="ordinal")

#> [1] -3.498327e-07 -1.143481e-10 -5.162827e+00
#> Warning message:
#> glm.fit: fitted rates numerically 0 occurred
    • Der code läuft immer noch so weit, obwohl ich denke, es ist eine lange Zeit. Meine Stichprobe ist nur 5,306 Beobachtungen.
    • Ich habe auch eine Frage bezüglich derid und time Variablen. Die Dokumentation auf multgee Staaten, dass id ist der Vektor, identifiziert die Cluster. In meinem Fall jedes Zwillingspaar ist ein cluster, so sollte nicht id = df$family_id_cass wo die Zwillinge haben die gleichen family_id_cass? Dies würde bedeuten, dass time = twin_order.
    • Allerdings, wenn ich intrinsic.pars() diese Weise bekomme ich eine Fehlermeldung: Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : contrasts can be applied only to factors with 2 or more levels. Meinst du meine interpretation der Parameter ist falsch?
    InformationsquelleAutor Jthorpe
Schreibe einen Kommentar