R: Ошибка в eval(predvars, data, env): числовой аргумент 'envir' не длиннее единицы в intrinsic.pars() в пакете multgee

Question

У меня есть фрейм данных dfсо следующими данными:

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

Я хотел бы бежать intrinsic.pars()так что я делаю следующее:

intrinsic.pars(df$ls75cat, df$family_id_cass, df$twin_order, rscale="ordinal") но я получаю следующую ошибку:

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

Что я делаю неправильно?

Answer 1

intrinsic.pars() кажется, написано с очень темпераментным интерфейсом, который эффективно игнорирует большинство параметров. (Я подозреваю, что это должна была быть внутренняя функция.) В любом случае вы должны предоставить data.frame для аргумента 'data', который включает переменные y,id, а также 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