用r,读入文件Hirsch_Cancer_Cell.tab,检验其第一列(列名为hr0_D2)的数值分布是否为正态分布,并用qqnorm画该列的分位比较图(Q-Q plot)。检验其第一列(列名为hr0_D2)的数值与最后一列(列名为hr36_D3)的数值分布是否一样,并用qqplot画这两列的分位比较图
时间: 2023-07-14 16:13:20 浏览: 92
MK_mannkendall_manner-kendall(m-k)_K._m-k检验法_Mk非参数检验_
以下是用R语言进行所述操作的代码:
首先,读入文件Hirsch_Cancer_Cell.tab:
```r
data <- read.delim("Hirsch_Cancer_Cell.tab")
```
然后,检验第一列的数值分布是否为正态分布:
```r
library(nortest)
x <- data$hr0_D2
ad.test(x) # 进行 Anderson-Darling 正态性检验
```
如果p值小于0.05,则可以拒绝原假设,即该分布不是正态分布。
接下来,用qqnorm画该列的分位比较图:
```r
qqnorm(x)
qqline(x)
```
如果该列的分布接近正态分布,则点应该大致分布在一条直线上。
最后,检验第一列和最后一列的数值分布是否一样:
```r
y <- data$hr36_D3
ks.test(x, y) # 进行 Kolmogorov-Smirnov 检验
```
如果p值小于0.05,则可以拒绝原假设,即两个分布不一样。
然后,用qqplot画这两列的分位比较图:
```r
qqplot(x, y)
abline(0, 1, col = "red")
```
如果这两列的分布相似,则点应该大致分布在一条直线上。
阅读全文