User Tools

Site Tools


21.05.08

e-mail: 21.05.08

Morge Heidi,


Hier noch ein R-Script mit einem Converter von AFLPDat zu dem Selection Programm von Matthieu.


Gruss
Tinu


R-Script:

convert_aflp<-function(input,output)
{
data=read.table(input,header=T)

outfile=output

nb_loci=ncol(data)-2

pops_names=levels(data[,2])
nb_pops=length(pops_names)

cat("[loci]=",nb_loci,"\n\n",file=outfile)

cat("[populations]=",nb_pops,"\n\n",file=outfile,append=T)

for (pop in 1:nb_pops)
{
  cur_pop=data[data[,2]==pops_names[pop],]
  nb_bands=rep(0,nb_loci)
  nb_individuals=rep(0,nb_loci)
  cat("[pop]=",pop,"\n",file=outfile,append=T)
  for (i in 1:nb_loci)
  {
    nb_bands[i]=sum(cur_pop[,i+2])
    nb_individuals[i]=sum(cur_pop[,i+2]==0)+sum(cur_pop[,i+2]==1)
    cat(i,nb_individuals[i],nb_bands[i],"\n",file=outfile,append=T)
  }
	cat("\n",file=outfile,append=T)
}
}
21.05.08.txt · Last modified: 2008/07/22 13:31 by 127.0.0.1