R-आउटपुट के साथ काम करना (2)

वैज्ञानिक शोध प्रश्न पर सोचते हुए, यदि IQ केवल जेनेटिक्स से प्रभावित होता है, तो हम उम्मीद करेंगे कि दोनों जुड़वाँ समूहों के बीच रिग्रेशन रेखा का स्लोप 1 हो. परिकल्पित स्लोप मान 1 का परीक्षण इस तरह किया जा सकता है कि एक नया टेस्ट स्टैटिस्टिक बनाया जाए, जो मापे कि प्रेक्षित स्लोप परिकल्पित मान 1 से कितना दूर है.

$$new_t = \frac{slope - 1}{SE}$$

यदि यह परिकल्पना सही है कि स्लोप 1 के बराबर है, तो नया टेस्ट स्टैटिस्टिक t-वितरण का अनुसरण करेगा, जिसका उपयोग हम p-value निकालने के लिए कर सकते हैं.

मॉडल से जैविक टर्म biological_term के रूप में उपलब्ध है.

यह अभ्यास पाठ्यक्रम का हिस्सा है

R में Linear Regression के लिए Inference

पाठ्यक्रम देखें

अभ्यास निर्देश

twins डेटासेट के degrees of freedom की गणना करें.
वैकल्पिक परिकल्पना के लिए दो-तरफ़ा p-value निकालें कि वास्तविक स्लोप 1 से भिन्न है. गणना को चरणों में बनाइए.
- टेस्ट स्टैटिस्टिक की गणना इस प्रकार करें: स्लोप estimate में से 1 घटाएँ, और पूरे को standard error से भाग दें.
- t-वितरण के cumulative distribution function pt() का उपयोग करते हुए, अभी निकाले गए degrees of freedom के साथ, test_statistic पर टेस्ट स्टैटिस्टिक की एक-तरफ़ा p-value निकालें.
- दो-तरफ़ा p-value को एक-तरफ़ा p-value का दोगुना लें.

इंटरैक्टिव व्यावहारिक अभ्यास

इस अभ्यास को इस नमूना कोड को पूरा करके आज़माएँ।

# Calculate the degrees of freedom of twins
degrees_of_freedom <- nrow(twins) - 2

biological_term %>%
  mutate(
    # Calculate the test statistic
    test_statistic = ___,
    # Calculate its one-sided p-value
    one_sided_p_value_of_test_statistic = ___,
    # ... and its two-sided p-value
    two_sided_p_value_of_test_statistic = ___
  )

कोड संपादित करें और चलाएँ