大家好，又见面了，我是你们的朋友全栈君。
题目：一名研究者用光子吸收法测量了妇女骨骼中无机物含量，对三根骨头主侧和非主侧记录了测量值，数据框“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
方法一：在一个坐标系里画出多个QQ图
R语言代码：
代码语言：javascript复制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")运行结果：
方法二：使用多个坐标系，把每个QQ图分布在一个界面中
R语言代码：
代码语言：javascript复制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
