Question

我有两个问题？

 > data<-read.table("UC.txt",header=TRUE, sep="\t")
 > data$tot<-data$P1+data$P2+data$P3+data$P4
 > head(data, 5)
      geno    P1    P2    P3    P4   tot
    1   G1 0.015 0.007 0.026 0.951 0.999
    2   G2 0.008 0.006 0.015 0.970 0.999
    3   G3 0.009 0.006 0.017 0.968 1.000
    4   G4 0.011 0.007 0.017 0.965 1.000
    5   G5 0.013 0.005 0.021 0.961 1.000

问题＃1：有时候，列数会有所不同，所以，如何将column2加到最后一列。类似数据[2]：data [n]

library("plyr")
> VD<-function(P4, tot){
   if(tot > 1) {return(P4-0.01)}
   if(tot < 1) {return(P4+0.01)}
   if(tot == 1) {return(P4)}
 }
> minu<-ddply(data, 'geno', summarize, Result=VD(P4, tot))
> v <- data$geno==minu$geno
> data[v, "P4"] <- minu[v, "Result"]
> data <- subset(data, select = -tot)
> data$tot<-data$P1+data$P2+data$P3+data$P4
> head(data, 5)
  geno   P1   P2   P3   P4 tot
1   G1 0.02 0.01 0.03 0.94   1
2   G2 0.01 0.01 0.02 0.96   1
3   G3 0.01 0.01 0.02 0.96   1
4   G4 0.01 0.01 0.02 0.96   1
5   G5 0.01 0.01 0.02 0.96   1

问题2：在这里，我需要通过将P1调整为P4来将'tot'舍入为1 条件：
1）我应该在P1到P4之间调整最大值 2）调整值可能不同，如0.01,0.001,0.0001。（它基于1-tot）
怎么做？
提前致谢

Answer 1

对于question1，要对除第一列以外的所有列求和：

dat$tot <- rowSums(dat[,-1])

R中的sum，roundoff和replace值

1 个答案: