आर मध्यस्थता विश्लेषण - बूटस्ट्रैपिंग

Question

मैं मध्यस्थ पैकेज का उपयोग कर आर में मध्यस्थता विश्लेषण करने का प्रयास कर रहा हूं। मैंने दस्तावेज को यह कैसे किया है, इस पर ध्यान दिया है, और आर द्वारा प्रदान किए गए उदाहरणों के माध्यम से पढ़ा है (यानी, मैंने पहले से ही "उदाहरण (मध्यस्थ)" चलाया है)। इसके बावजूद, मुझे चलाने के लिए सबसे आसान मध्यस्थता नहीं मिल सकती है। आदर्श रूप में, मैं बूटस्ट्रैपिंग प्रक्रिया करना चाहता हूं, एक ला प्रचारक & हेस (2004)।

model.m <- lm(desirdata1$zpers1 ~ desirdata1$zdesir1 + desirdata1$age) 
model.y <- lm(desirdata1$zpers1 ~ desirdata1$age) 
age1test <- mediate(model.m, model.y,treat="age", mediator="zdesir1", 
        boot=TRUE, sims=50) 

ध्यान दें कि डाटासेट desirdata कहा जाता है परिणाम zpers1 कहा जाता है और मध्यस्थ "zdesir1" कहा जाता है, उपचार age कहा जाता है,:

यहाँ कोड है कि मैं चलाने के लिए कोशिश कर रहा हूँ है।

Error in `[.data.frame`(m.data, , treat) : undefined columns selected 

यह दावा करते हुए किया जाना चाहिए कि एक चर (विशेष रूप से, उपचार चर) मौजूद नहीं है लगता है: जब मैं इस चलाने के लिए, मैं निम्नलिखित त्रुटि मिलती है। हालांकि, चल रहे नाम (desirdata) से पता चलता है कि चर वहां है, और यह सही रूप से नामित है, जैसा कि अन्य सभी चर हैं। पहले दो मॉडल (model.m और model.y) ठीक चलते हैं, और आउटपुट दिखता है जैसा कि इसे करना चाहिए। यह केवल मध्यस्थता मॉडल है जिसे मैं चलाने के लिए नहीं मिल सकता। मैंने एक टाइपो नहीं बनाया है, जहां तक ​​मैं कह सकता हूं, और मैंने इसे एक सौ बार चेक किया है।

विचार?

Answer 1

जैसा कि मैंने प्रलेखन में उदाहरण पढ़ा है, मध्यस्थ मॉडल के लिए model.m मुख्य प्रतिगमन वस्तु model.y के लिए उससे अलग परिणाम होगा। जब से तुम पृष्ठभूमि का वर्णन किया है नहीं और डेटा किस तरह यह इस बात का बहुत कुछ होने के लिए कठिन है, लेकिन सोच आप टाइप करने के लिए होती है, तो:

model.m <- lm(zdesir1 ~ age, data=desirdata1) 
model.y <- lm(zpers1 ~ age, , data=desirdata1) 
age1test <- mediate(model.m, model.y,treat="age", mediator="zdesir1", 
        boot=TRUE, sims=50) 

मैं formula और data वस्तुओं का उपयोग कर इसे डाली, कुछ प्रतिगमन कार्यों के बाद से केवल वैक्टर दिए जाने पर तोड़ दें। टाइपो को देखना भी आसान बनाता है।

Answer 2

एमबीईएस पैकेज आज़माएं। प्रचारक इसकी सिफारिश करता है और आप केवल मध्यस्थता समारोह का उपयोग कर सकते हैं। यदि आप बूटस्ट्रैप चाहते हैं तो बस सुनिश्चित करें कि यह बूटस्ट्रैप = ​​TRUE कहता है। बी बूटस्ट्रैप्स की संख्या है।

mediation(x, mediator, dv, S = NULL, N = NULL, x.location.S = NULL, 
mediator.location.S = NULL, dv.location.S = NULL, mean.x = NULL, 
mean.m = NULL, mean.dv = NULL, conf.level = 0.95, 
bootstrap = FALSE, B = 1000, which.boot="both", save.bs.replicates=FALSE) 

Answer 3

मैं तो मैं डिबग (मध्यस्थता) भाग गया और पाया जहां समस्या थी नकली डेटा के साथ एक ही समस्या में पड़ गए। मेरा मानना ​​है कि समस्या कोड के [इलाज = "आयु", मध्यस्थ = "zdesir1"] भाग के साथ है। यदि आप डेटा संलग्न करते हैं, तो आपको उस समस्या में भाग नहीं लेना चाहिए। वैकल्पिक रूप से, आप [treat = "desirdata1 $ age", mediator = "desirdata1 $ zdesir1"] का उपयोग कर सकते हैं जो समस्या को हल करना चाहिए।

Answer 4

खैर इस प्रयास करें:

model.m <- lm(zdesir1 ~ age, data=desirdata1) 
model.y <- lm(zpers1 ~ age + zdesir1, data=desirdata1) 
age1test <- mediate(model.m, model.y,treat="age", mediator="zdesir1", 
        boot=TRUE, sims=50) 

यह आसान बनाने के लिए, मध्यस्थ मॉडल (model.m) परिणाम के रूप में मध्यस्थ होना चाहिए।

Answer 5

आपके मॉडल सही नहीं हैं। model.m को चौथाई से मध्यस्थ की भविष्यवाणी करनी चाहिए, और model.y मध्यस्थ और चतुर्थ से DV की भविष्यवाणी करनी चाहिए।

model.m <- lm(desirdata1$zdesir1 ~ desirdata1$age) 
model.y <- lm(desirdata1$zpers1 ~ desirdata1$zdesir1 + desirdata1$age) 
age1test <- mediate(model.m, model.y, treat="age", mediator="zdesir1", boot=TRUE, sims=50)