如何從sempls {semPLS}（或從densityplot {semPLS}）中提取殘差？

Question

我通過semPLS包中的函數sempls運行偏最小二乘路徑模型（混合了形成性和反射性潛變量）。我想檢查模型殘差（由潛在變量之間的普通最小二乘迴歸產生）中的空間自相關，但必須首先弄清楚如何提取它們。 link to data file

land2011g <- read.table("percover2011.csv", header=T, sep=',', row.names=1) 
land2011<-land2011g[c(1:2,4:16),] 
topaphids<-land2011[c('A.cracc','R.padi','A.nerii','T.trifolii','R.maidis', 
'Tetr.sp.','Pemphigus.sp.')] 
r1 <- land2011[c(1:3,128:131)] 
r1.s <- as.data.frame(scale(r1, center=T, scale=T)) 
r2 <- land2011[c(11:13,132:135)] 
r2.s <- as.data.frame(scale(r2, center=T, scale=T)) 
r3 <- land2011[c(21:23,136:139)] 
r3.s <- as.data.frame(scale(r3, center=T, scale=T)) 
r4 <- land2011[c(31:33,140:143)] 
r4.s <- as.data.frame(scale(r4, center=T, scale=T)) 
r5 <- land2011[c(41:43,144:147)] 
r5.s <- as.data.frame(scale(r5, center=T, scale=T)) 
weediness<-land2011[c('largest.total.cover.m2')] 
noncols<-land2011[c('total.noncols')] 

ifrom <- c("r1","r2","r3","r4","r5","weed.cover","aphids","aphids") 
ito<-c("aphids","aphids","aphids","aphids","aphids",'aphids',"prsv","wmv") 
inner.mod<-data.frame(ifrom,ito) 
inner.mat<-as.matrix(inner.mod) 

ofrom<-c('r1.corn.wheat','r1.soy','r1a.pasture.hay.grass', 
'r1a.urban','r1a.forest.shrub','r2.corn.wheat','r2.soy', 
'r2a.pasture.hay.grass','r2a.urban','r2a.forest.shrub',  
'r3.corn.wheat','r3.soy','r3a.pasture.hay.grass','r3a.urban', 
'r3a.forest.shrub','r4.corn.wheat','r4.soy','r4a.pasture.hay.grass', 
'r4a.urban','r4a.forest.shrub', 'r5.corn.wheat','r5.soy', 
'r5a.pasture.hay.grass','r5a.urban','r5a.forest.shrub', 
'weed.cover','aphids','prsv','wmv') 
oto<-c('r1','r1','r1','r1','r1','r2','r2','r2','r2','r2','r3','r3','r3', 
'r3','r3','r4','r4','r4','r4','r4','r5','r5','r5','r5','r5', 
'largest.total.cover.m2','total.noncols','arcsin.sqrt.PRSV', 
'arcsin.sqrt.WMV') 
outer.mod<-data.frame(ofrom,oto) 
outer.mat<-as.matrix(outer.mod) 

prsv<-land2011[c("arcsin.sqrt.PRSV")] 
wmv<-land2011[c("arcsin.sqrt.WMV")] 
subsets <- cbind(r1.s,r2.s,r3.s,r4.s,r5.s,weediness,topaphids, 
noncols,prsv,wmv) 

plspm11<-plsm(data=subsets, strucmod=inner.mat, measuremod=outer.mat) 

pmodp11<-sempls(plspm11, data=subsets, method="pearson", scaled=TRUE, 
maxit=100, convCrit="square") 

殘差（pmodp11）函數返回 「NULL」 值。

semPLS包中的densityplot函數將爲我繪製殘差，因此我假設它們必須存儲在sempls對象中，但我不確定如何訪問它們。

densityplot(pmodp11, subsets, use=c("residuals")) 

density plot of sempls object pmodp11

我還發現，在semPLS相關殘差預測下面的代碼，但嘗試使用它給了我這個錯誤：

Error in UseMethod("predict") : no applicable method for 'predict' applied to an object of class "sempls" .

residuals.sempls <- function(object, what=c("LVs", "MVs"), scale=c("scaled", "original"), total=FALSE){ 
    what <- match.arg(what) 
    scale <- match.arg(scale) 
    model <- object$model 
    data <- object$data 
    # LVs 
    if(what=="LVs"){ 
     res <- object$factor_scores - predict(object, what=what, scale=scale, total=total) 
    } 
    # MVs 
    else{ 
     if(scale=="scaled"){ 
     pdata <- predict(object, what=what, scale=scale, total=total) 
     res <- data - pdata 
     } 
     else{ 
     data <- rescale(data) 
     pdata <- predict(object, what=what, scale=scale, total=total) 
     res <- data - pdata 
     } 
    } 
    return(res) 
} 

Answer 1

功能後預測。 sempls（在R-Forge_predict.sempls上定義）插入到residuals.sempls函數中，residuals.sempls現在似乎可以與sempls對象一起使用以返回m奧德爾殘差。

residuals.sempls <- function(object, what=c("LVs", "MVs"), scale=c("scaled", "original"), total=FALSE){ 
    what <- match.arg(what) 
    scale <- match.arg(scale) 
    model <- object$model 
    data <- object$data 
    # LVs 
    if(what=="LVs"){ 
     res <- object$factor_scores - predict.sempls(object, what=what, scale=scale, total=total) 
    } 
    # MVs 
    else{ 
     if(scale=="scaled"){ 
     pdata <- predict.sempls(object, what=what, scale=scale, total=total) 
     res <- data - pdata 
     } 
     else{ 
     data <- rescale(data) 
     pdata <- predict.sempls(object, what=what, scale=scale, total=total) 
     res <- data - pdata 
     } 
    } 
    return(res) 
}