2017-04-11 1 views
0

Je construction d'un modèle d'équation structurelle par morceaux en utilisant le paquet piecewiseSEM dans R (Lefcheck - https://cran.r-project.org/web/packages/piecewiseSEM/vignettes/piecewiseSEM.html)extraire des coefficients de chemin de morceaux SEM (modèle d'équation structurelle)

I déjà créé le modèle défini et je puisse évaluer le modèle ajustement, de sorte que le modèle lui-même fonctionne. De plus, les données correspondent au modèle (p = 0,528).

Mais je n'arrive pas à extraire les coefficients de chemin . Ceci est l'erreur que je reçois: Error in cbind(Xlarge, Xsmall) : number of rows of matrices must match (see arg 2)

Je l'ai déjà essayé (mais cela ne fonctionne pas):

  • normaliser mes données en raison de l'avertissement: Some predictor variables are on very different scales: consider rescaling

  • adapté mes données (jeté certaines valeurs NA) loin

Ceci est mon modellist:

predatielijst = list(

    lmer(plantgrootte ~ gapfraction + olsen_P + (1|plot_ID), data = d), 

    glmer(piek1 ~ gapfraction + olsen_P + plantgrootte + (1|plot_ID), 
     family = poisson, data = d), 

     glmer(predatie ~ piek1 + (1|plot_ID), family = binomial, data = d) 
) 

avec "predatie" étant une variable binaire (oui ou non) et toutes les autres variables continues (gapfraction, plantgrootte, olsen_P & piek1)

Merci à l'avance!

+0

utilisez 'sem.coefs' pour extraire les coefficients de chemin! – jslefche

Répondre

0

Essayez d'installer la version de développement:

library(devtools) install_github("jslefche/[email protected]")

Remplacer list avec psem et exécuter la fonction coefs ou summary. Il va probablement se débarrasser de votre erreur. Sinon, ouvrez un bug sur Github! AVERTISSEMENT: cette opération écrase votre version actuelle de CRAN. Vous devrez réinstaller à partir de CRAN pour récupérer la version 1.4.

+0

Un problème est survenu avec les valeurs NA. Après les avoir supprimés, j'ai réussi à extraire les coefficients de chemin et à faire la figure avec le paquet 'old'. – geum

0

essayez d'utiliser lme (hors de la bibliothèque nlme) ilstead de glmer. Pour autant que je comprenne, le fait que lmer ne fournisse pas de valeurs p (alors que c'est le cas) semble être le problème ici. Espérons que cela fonctionne.