我比较新，awk所以我有一个关于做除法并将结果打印在新列中的简单问题。例如：
head data
1 13273 . G C 563 5 . 25 128
1 202259 . G T 675 8 . 12 130
1 598934 . C C 756 9 . 17 231
1 634112 . T C 125 1 . 32 89
1 779762 . G A 675 5 . 28 187
我想将第 9 列除以第 10 列并将结果打印在新的第 11 列中，最好将新结果从高到低排序。例如：
1 634112 . T C 125 1 . 32 89 0.360
1 13273 . G C 563 5 …我有一个简单的问题.如何将数据帧转换为Fisher精确测试的列联表？
我有data大约19000行:
head(data)
R_T1 R_T2 NR_T1 NR_T2
GMNN 14 60 70 157
GORASP2 7 67 39 188
TTC34 5 69 41 186
ZXDC 8 66 37 190
ASAH2 9 65 46 181
我想将每一行转换为列联表来执行Fisher精确检验.例如,对于GMNN:
R NR
T1 14 70
T2 60 157
fisher.test(GMNN, alternative="two.sided")
Fisher's Exact Test for Count Data
data: GMNN
p-value = 0.05273
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
0.2531445 1.0280271
sample estimates:
odds ratio
0.5243787
由于我有19000行数据,我宁愿输出为: …
我有一个my.list看起来像这样的列表
$S1
A B C D
1 101027 NA 0.48 NA
2 101031 1.50 1.30 0.8666667
3 101032 1.40 0.78 0.5571429
4 101127 NA NA NA
5 101220 9.30 7.30 0.7849462
$S2
A B C D
1 102142 NA 0.45 NA
2 102143 0.70 1.20 1.7142857
3 102144 NA 0.44 NA
4 102148 0.45 NA NA
5 102151 0.91 0.64 0.7032967
6 102152 0.78 NA NA
我想删除NA列表数据框中的所有行，使其看起来像
$S1
A B C D
2 101031 1.50 1.30 …