ইন Rআমি ব্যবহার heteroscedasticity জন্য Breusch-পৌত্তলিক পরীক্ষা সম্পাদন করতে পারবেন ncvTestফাংশন carপ্যাকেজ। একটি ব্রুশ – পৌত্তলিক পরীক্ষা হ'ল এক ধরণের চি স্কোয়ার পরীক্ষা।
আমি কীভাবে এই ফলাফলগুলি ব্যাখ্যা করব:
> require(car)
> set.seed(100)
> x1 = runif(100, -1, 1)
> x2 = runif(100, -1, 1)
> ncvTest(lm(x1 ~ x2))
Non-constant Variance Score Test
Variance formula: ~ fitted.values
Chisquare = 0.2343406 Df = 1 p = 0.6283239
> y1 = cumsum(runif(100, -1, 1))
> y2 = runif(100, -1, 1)
> ncvTest(lm(y1 ~ y2))
Non-constant Variance Score Test
Variance formula: ~ fitted.values
Chisquare = 1.191635 Df = 1 p = 0.2750001