大家好,又见面了,我是你们的朋友全栈君。
题目:一名研究者用光子吸收法测量了妇女骨骼中无机物含量,对三根骨头主侧和非主侧记录了测量值,数据框“T1bones.txt”中的第2至第7列记录了相应数据。对各个变量做qq图,在同一个界面画出所有的qq图,不同的qq图用不同颜色表示。
文件:T1bones.txt
受试者编号 主侧桡骨 桡骨 主侧肱骨 肱骨 主侧尺骨 尺骨 1 1.103 1.052 2.139 2.238 0.873 0.872 2 0.842 0.859 1.873 1.741 0.590 0.744 3 0.925 0.873 1.887 1.809 0.767 0.713 4 0.857 0.744 1.739 1.547 0.706 0.674 5 0.795 0.809 1.734 1.715 0.549 0.654 6 0.787 0.779 1.509 1.474 0.782 0.571 7 0.933 0.880 1.695 1.656 0.737 0.803 8 0.799 0.851 1.740 1.777 0.618 0.682 9 0.945 0.876 1.811 1.759 0.853 0.777 10 0.921 0.906 1.954 2.009 0.823 0.765 11 0.792 0.825 1.624 1.657 0.686 0.668 12 0.815 0.751 2.204 1.846 0.678 0.546 13 0.755 0.724 1.508 1.458 0.662 0.595 14 0.880 0.866 1.786 1.811 0.810 0.819 15 0.900 0.838 1.902 1.606 0.723 0.677 16 0.764 0.757 1.743 1.794 0.586 0.541 17 0.733 0.748 1.863 1.869 0.672 0.752 18 0.932 0.898 2.028 2.032 0.836 0.805 19 0.856 0.786 1.390 1.324 0.578 0.610 20 0.890 0.950 2.187 2.087 0.758 0.718 21 0.688 0.532 1.650 1.378 0.533 0.482 22 0.940 0.850 2.334 2.225 0.757 0.731 23 0.493 0.616 1.037 1.268 0.546 0.615 24 0.835 0.752 1.509 1.422 0.618 0.664 25 0.915 0.936 1.971 1.869 0.869 0.868
R语言代码:
data <-read.table("T1bones.txt",header=TRUE);data
d1 <- data[1:25,2];d1
d2 <- data[1:25,3];d2
d3 <- data[1:25,4];d3
d4 <- data[1:25,5];d4
d5 <- data[1:25,6];d5
d6 <- data[1:25,7];d6
qqnorm(d1)
qqline(d1)
par(new=T)
qqnorm(d2,col="green",xaxt="n",yaxt="n")
qqline(d2,col="green")
par(new=T)
qqnorm(d3,col="blue",xaxt="n",yaxt="n")
qqline(d3,col="blue")
par(new=T)
qqnorm(d4,col="yellow",xaxt="n",yaxt="n")
qqline(d4,col="yellow")
par(new=T)
qqnorm(d5,col="red",xaxt="n",yaxt="n")
qqline(d5,col="red")
par(new=T)
qqnorm(d6,col="orange",xaxt="n",yaxt="n")
qqline(d6,col="orange")
运行结果:
R语言代码:
par(mfcol=c(2,3))
qqnorm(d1,main="d1")
qqline(d1)
qqnorm(d2,col="green",main="d2")
qqline(d2,col="green")
qqnorm(d3,col="blue",main="d3")
qqline(d3,col="blue")
qqnorm(d4,col="yellow",main="d4")
qqline(d4,col="yellow")
qqnorm(d5,col="red",main="d5")
qqline(d5,col="red")
qqnorm(d6,col="orange",main="d6")
qqline(d6,col="orange")
运行结果:
相关资料:数据分析与R软件第二版(李素兰著)每章例题代码和数据
发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/169348.html原文链接:https://javaforall.cn