मिश्रित मॉडल की भविष्यवाणियों को एक पंक्ति के रूप में आर में प्लॉट करना

Nov 24 2020

मैं predict()अपने नीचे के मिश्रित मॉडल की भविष्यवाणियों ( ) की साजिश करने की कोशिश कर रहा हूं, ताकि मैं अपनी अवधारणा को नीचे की रेखा के रूप में वांछित साजिश प्राप्त कर सकूं।

मैंने अपने मॉडल की भविष्यवाणियों को प्लॉट करने की कोशिश की है, लेकिन मैं अपना वांछित प्लॉट हासिल नहीं कर पाया। क्या परिभाषित करने का एक बेहतर तरीका है predict()ताकि मैं अपना वांछित भूखंड प्राप्त कर सकूं?

library(lme4)

dat3 <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/dat3.csv')

m4 <- lmer(math~pc1+pc2+discon+(pc1+pc2+discon|id), data=dat3)

newdata <- with(dat3, expand.grid(pc1=unique(pc1), pc2=unique(pc2), discon=unique(discon)))

y <- predict(m4, newdata=newdata, re.form=NA)

plot(newdata$pc1+newdata$pc2, y)

जवाब

2 MagnusNordmo Nov 24 2020 at 20:21

अधिक sjPlotgridएक प्लॉट में कई भविष्यवाणियों को लपेटने के लिए पैरामीटर देखें ।

library(lme4)
library(sjPlot)
library(patchwork)

dat3 <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/dat3.csv')

m4 <- lmer(math~pc1+pc2+discon+(pc1+pc2+discon|id), data=dat3) # Does not converge

m4 <- lmer(math~pc1+pc2+discon+(1|id), data=dat3) # Converges

# To remove discon
a <- plot_model(m4,type = 'pred')[[1]]
b <- plot_model(m4,type = 'pred',title = '')[[2]]

a + b 

संपादित करें 1: मुझे ढांचे के dropconभीतर शब्द को हटाने में कुछ परेशानी थी sjPlot। मैं हार मान गया और वापस लौट आया patchwork। मुझे यकीन है कि डैनियल सही तरीका जान सकता था।

Oliver Nov 24 2020 at 20:14

जैसा कि मैग्नस नॉर्डमो सुझाव देते हैं, यह बहुत सरल है sjPlotजिसके साथ इस प्रकार के कथानक के लिए कुछ पूर्वनिर्धारित कार्य हैं।

library(lme4)

dat3 <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/dat3.csv')

m4 <- lmer(math~pc1+pc2+discon+(pc1+pc2+discon|id), data=dat3)

plot_model(m4, type = 'pred', terms = c('pc1', 'pc2'),
           ci.lvl = 0)

जो निम्नलिखित परिणाम देता है।

इस साजिश में दूसरे कार्यकाल के विभिन्न quantiles शामिल करने के लिए डिज़ाइन किया गया है termsकी कुल्हाड़ियों से अधिक pc1और pred। आप इन भूखंडों को विभाजित कर सकते हैं और उन्हें जोड़ patchworkसकते हैं और अंतराल को वर्ग कोष्ठक का उपयोग करके बदला जा सकता है terms(जैसे pc1 [-10:1]-10 और 1 के बीच अंतराल के लिए)।