我用的是lavvan package。網上搜索之後發現lavvan沒有自帶關於潛變量的交互作用測量。一開始查到了nlsem包的lav2nlsem命令,sample code比較複雜,是兩個潛變量之間的交互作用,我試圖改寫了下
MOD<-'
C=~ X1 + X2
Y~ C + X3 + C:X3
'
model <- lav2nlsem(MOD)
但在運行的時候報錯
Error in dimnames(x) <- dn :
length of 'dimnames' [2] not equal to array extent
但我確定我的dataset沒有錯,因為可以在lavvan下順利跑其他的regression model.
我的問題是:1)lav2nlsem是否能實現我的需求?如果是,code應該怎麼寫?
2)如果lav2nlsem不能實現我的需求,還有其他什麼package或者function可以實現我的需求?假設我的變量有:因變量Y,自變量X1,X2,X3, 然後X1,X2屬於潛變量C,交互作用是潛變量X和自變量X3之間,那麼code應該怎樣寫?
非常感謝!