2017-04-11 86 views
0

我構建使用piecewiseSEM包R(Lefcheck - https://cran.r-project.org/web/packages/piecewiseSEM/vignettes/piecewiseSEM.html)分段結構方程模型分段提取SEM(結構方程模型)的路徑係數

我已經創建的模型設定,我可以評估模型適合,所以模型本身的作品。此外,數據符合模型(p = 0.528)。

但我沒有成功提取路徑係數。 這是錯誤我得到:Error in cbind(Xlarge, Xsmall) : number of rows of matrices must match (see arg 2)

我已經嘗試過(但是這沒有工作):

  • 規範,因爲警告我的數據:Some predictor variables are on very different scales: consider rescaling

  • 適應我的數據(扔一些NA值)

這是我的造型師:

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) 
) 

與 「predatie」 是一個二元變量(是或否)和所有其餘的連續變量(gapfraction,plantgrootte,olsen_P & piek1)

提前感謝!

+0

使用'sem.coefs'來提取路徑係數! – jslefche

回答

0

嘗試安裝開發版:

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

更換listpsem和運行coefssummary功能。它可能會擺脫你的錯誤。如果沒有,請在Github上打開一個錯誤!

警告:這會從CRAN覆蓋您當前的版本。您需要從CRAN重新安裝才能獲得版本1.4。

+0

NA值出現問題。刪除它們之後,我成功地提取了路徑係數並使用「舊」包製作了圖形。 – geum

0

嘗試使用lme(出nlme庫)ilstead glmer。據我所知,lmer不提供p值(而lme確實)的事實似乎是這裏的問題。 希望這個作品。