Revision c52a091b2759b18afd9f7383ec23eee1d6aa9b84 authored by Hervé Seitz on 18 June 2020, 16:35:53 UTC, committed by GitHub on 18 June 2020, 16:35:53 UTC
1 parent 83736c9
R_measures_editing
args <- commandArgs(TRUE)
CUTOFF=20
data=read.table(paste('Quality_in_',args[1],'.dat',sep=''),header=T)
genomic_As=read.table('Adenines_in_genomic_sequence.dat')
memo=c()
for (x in genomic_As$V1)
{
nt=data[,1+x]
q=data[,151+x]
memo=rbind(memo,c(x,length(nt[q>=CUTOFF]),length(nt[q>=CUTOFF & nt=='A']),length(nt[q>=CUTOFF & nt=='C']),length(nt[q>=CUTOFF & nt=='G']),length(nt[q>=CUTOFF & nt=='T'])))
}
write.csv(memo,paste('Editing_data_',args[1],'.csv',sep=''))
Computing file changes ...
