Я довольно новичок в R, и я пытаюсь запустить пакет twinlm с mets. Перед попыткой сбора данных я выделил красный https://rdrr.io/rforge/mets/src/inst/misc/twinbmi.R Вот код моих данных: lt <- twinlm (Temp1 ~ Возраст + пол, id = "Family_No", zyg = "zyg", DZ = "DZ", data = temp, type = "sat") </p>
но я получаю эту ошибку Ошибка: необходимо использовать вектор в [, не являющийся объектом матрицы классов. Запустите rlang::last_error(), чтобы увидеть, где произошла ошибка.
, если I Backtrace: x 1. -mets :: twinlm (...) 2. -base :: as.factor (zygstat) 3. -base :: factor (x) 4. -base :: order (y) 5. -base :: lapply (z, function (x), если (is.object (x)) as.vector (xtfrm (x)), иначе x ) 6. -base ::: FUN (X [[i]], ...) 7. + -base :: as.vector (xtfrm (x)) 8. + -base :: xtfrm (x) 9. -base :: xtfrm.default (x) 10. + -base :: as.vector (rank (x, ties.method = "min", na.last = "keep")) 11. -base :: rank ( x, ties.method = "min", na.last = "keep") 12. + -x [! nas] 13. -tibble ::: [.tbl_df (x,! nas) 14. -tibble ::: check_names_df (i, x)
Я читал, что некоторые inputa требуют data.frame и не работают с tibble, и это может произойти, когда есть подмножество dataframe.
Infact вот мой заголовок данных (временная шкала)
Тиббл: 6 x 12
УчастникID Temp1 Temp2 Temp3 Family_No Family_No.1 YEAR_BIRTH SEX zygosity
1 21 36,9 36,7 36,7 2 2 1940 F DZ
2 22 36,8 36,6 36,9 2 2 1940 F DZ
3 131 36,1 36,2 36,3 13 13 1952 F MZ
4 132 36,2 36,3 36,4 13 13 1952 F MZ
5 241 37,1 37,1 37 24 24 1925 F DZ
6 242 36,7 36,9 37 24 24 1925 F DZ
... с 3 другие переменные: пол, zyg, возраст
Поэтому я попытался преобразовать его в data.frame
temp <- as.data.frame (temp) head (temp) ParticipantID Temp1 Temp2 Temp3 Family_No Family_No.1 YEAR_BIRTH SEX zygosity 1 21 36,9 36,7 36,7 2 2 1940 F DZ 2 22 36,8 36,6 36,9 2 2 1940 F DZ 3 131 36,1 36,2 36,3 13 13 1952 F MZ 4 132 36,2 36,3 36,4 13 13 1952 F MZ 5 241 37,1 37,1 37,0 24 24 1925 F DZ 6 242 36,7 36,9 37,0 24 24 1925 F DZ пол zyg Возраст 1 0 1 80 2 0 1 80 3 0 2 68 4 0 2 68 5 0 1 95 6 0 1 95 </p>
Но теперь я получаю эту ошибку: Ошибка в procdata.lvm (x, data = data, na.method = na.method): Нет наблюдений
если бы кто-нибудь мог мне помочь, я бы действительно ценить это!