R: Wie zu extrahieren r Quadrat aus lmekin in Coxme-Paket

Question

Ich versuche herauszufinden, wie gut mein gemischtes Modell mit Familie Effekt passt die Daten. Ist es möglich, r-Quadrat-Werte aus lmekin-Funktionen zu extrahieren? Und wenn ja, ist es möglich, partielle r-Quadrat-Werte für jede der Kovariablen zu extrahieren? Beispiel:

model= lmekin(formula = height ~ score + sex + age + (1 | IID), data = phenotype_df, varlist = kinship_matrix) 

Ich habe das Mumin-Paket versucht, aber es scheint nicht mit lmekin Modellen zu arbeiten. Vielen Dank.

Answer 1

ich in der Lage bin, die r.squaredLR() Funktion zu nutzen,

library(coxme) 
library(MuMIn) 
data(ergoStool, package="nlme") # use a data set from nlme 
fit1 <- lmekin(effort ~ Type + (1|Subject), data=ergoStool) 
r.squaredLR(fit1) 

(Ich bin ziemlich sicher, dass das funktioniert, aber eine Sache, die groß ist, ist zu tun, ein reproduzierbares Beispiel zu erstellen, so kann ich Ihren Code ausführen zu überprüfen Sie zum Beispiel, ich bin nicht genau sicher, was phenotype_df aussieht, und ich bin nicht in der Lage, Ihren Code zu laufen, wie es ist, eine große Ressource dafür ist das reprex Paket).