كيفية رسم توزيع t في r

لرسم دالة الكثافة الاحتمالية للتوزيع في R، يمكننا استخدام الوظائف التالية:

• dt(x, df) لإنشاء دالة كثافة الاحتمال
• منحنى (وظيفة، من = NULL، إلى = NULL) لرسم دالة كثافة الاحتمال

لرسم دالة الكثافة الاحتمالية، نحتاج إلى تحديد df (درجات الحرية) في الدالة dt() بالإضافة إلى القيم من وإلى في الدالة المنحنى() .

على سبيل المثال، يوضح الكود التالي كيفية رسم دالة كثافة الاحتمال لتوزيع بـ 10 درجات حرية حيث يكون المحور السيني للرسم بين -4 و4:

 curve(dt(x, df=10), from=-4, to=4)

على غرار التوزيع الطبيعي ، يكون توزيع t متماثلًا حول متوسط قدره 0.

يمكننا إضافة عنوان، وتغيير تسمية المحور Y، وزيادة عرض الصف، وحتى تغيير لون الصف لجعل الحبكة أكثر جمالية:

 curve(dt(x, df=10), from=-4, to=4, 
    main = 't Distribution (df = 10)', #add title
    ylab = 'Density', #change y-axis label
    lwd = 2, #increase line width to 2
    col = 'steelblue') #change line color to steelblue 

يمكننا أيضًا إضافة منحنيات متعددة إلى الرسم البياني لمقارنة توزيعات t بدرجات مختلفة من الحرية. على سبيل المثال، يقوم التعليمة البرمجية التالية بإنشاء مخططات توزيع t مع df = 6، وdf = 10، وdf = 30:

 curve(dt(x, df=6), from=-4, to=4, col='blue') 
curve(dt(x, df=10), from=-4, to=4, col='red', add=TRUE)
curve(dt(x, df=30), from=-4, to=4, col='green', add=TRUE)

يمكننا إضافة وسيلة إيضاح إلى المخطط باستخدام وظيفة legend() ، والتي تأخذ الصيغة التالية:

وسيلة الإيضاح (x، y=NULL، legend، fill، col، bg، lty، cex)

ذهب:

• x وy: إحداثيات x وy المستخدمة لتحديد موضع وسيلة الإيضاح
• وسيلة الإيضاح: النص الذي سيتم وضعه في وسيلة الإيضاح
• ملء: ملء اللون داخل وسيلة الإيضاح
• col: قائمة الألوان التي سيتم استخدامها للخطوط الموجودة داخل وسيلة الإيضاح
• bg: لون خلفية وسيلة الإيضاح
• lty: نمط الخط
• cex: حجم النص في وسيلة الإيضاح

في مثالنا سوف نستخدم بناء الجملة التالي لإنشاء وسيلة إيضاح:

 #create density plots
curve(dt(x, df=6), from=-4, to=4, col='blue') 
curve(dt(x, df=10), from=-4, to=4, col='red', add=TRUE)
curve(dt(x, df=30), from=-4, to=4, col='green', add=TRUE)

#add legend
legend(-4, .3, legend=c("df=6", "df=10", "df=30"),
       col=c("blue", "red", "green"), lty=1, cex=1.2)