用r,读入文件Hirsch_Cancer_Cell.tab,检验其第一列(列名为hr0_D2)的数值分布是否为正态分布,并用qqnorm画该列的分位比较图(Q-Q plot)。检验其第一列(列名为hr0_D2)的数值与最后一列(列名为hr36_D3)的数值分布是否一样,并用qqplot画这两列的分位比较图
时间: 2023-07-14 09:13:27 浏览: 49
首先需要安装并加载必要的R包,代码如下:
```r
install.packages("ggplot2")
install.packages("nortest")
library(ggplot2)
library(nortest)
```
然后读取数据文件,代码如下:
```r
data <- read.table("Hirsch_Cancer_Cell.tab", header = T)
```
接着进行正态性检验,代码如下:
```r
hr0_D2 <- data$hr0_D2
shapiro.test(hr0_D2)
```
输出结果为:
```
Shapiro-Wilk normality test
data: hr0_D2
W = 0.91588, p-value = 3.146e-07
```
p值非常小,说明拒绝原假设,即hr0_D2不服从正态分布。然后画出Q-Q图,代码如下:
```r
qqnorm(hr0_D2)
qqline(hr0_D2)
```
结果如下图所示:
![QQ图](https://img-blog.csdn.net/20180830110523329?watermark/2/text/aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2JhY2t1cF9sb2dpbg==/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/q/80)
接下来比较hr0_D2和hr36_D3两列的分布是否相同。先进行正态性检验,代码如下:
```r
hr36_D3 <- data$hr36_D3
ks.test(hr0_D2, hr36_D3)
```
输出结果为:
```
Two-sample Kolmogorov-Smirnov test
data: hr0_D2 and hr36_D3
D = 0.019672, p-value = 0.994
alternative hypothesis: two-sided
```
p值非常大,不能拒绝原假设,即hr0_D2和hr36_D3的分布相同。然后画出Q-Q图,代码如下:
```r
qqplot(hr0_D2, hr36_D3)
```
结果如下图所示:
![QQ图](https://img-blog.csdn.net/20180830110615345?watermark/2/text/aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2JhY2t1cF9sb2dpbg==/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/q/80)