cd /seppdata/sepp/linkage/release/ /system/apps/biosoft/R-2.14.0/bin/R library(fabia) source("funct.R") indit <- read.table("phase1_integrated_calls.20101123.ALL.panel", header = FALSE, sep = "\t", quote = "",as.is=TRUE) indi1 <- as.vector(unlist(rbind(indit[,1],indit[,1]))) indi2 <- as.vector(unlist(rbind(indit[,2],indit[,2]))) indi3 <- as.vector(unlist(rbind(indit[,3],indit[,3]))) indi4 <- as.vector(unlist(rbind(indit[,4],indit[,4]))) indi4 <- gsub(",",";",indi4) indi <- cbind(indi1,indi2,indi3,indi4) save(indi,file="individuals.Rda") # EUeuropeans <- c(which(indi[,2]=="CEU"),which(indi[,2]=="FIN"),which(indi[,2]=="GBR"),which(indi[,2]=="IBS"),which(indi[,2]=="TSI")) EUeuropeans <- c(which(indi[,3]=="EUR")) EUeurop <- as.integer(sort.int(as.integer(unique(EUeuropeans)))) indi_EUeurop <- indi[EUeurop,] save(indi_EUeurop,file="individualsEU.Rda") annot <- read.table("ALL.chr1.merged_beagle_mach.20101123.snps_indels_svs.genotypes_annot.txt",header = FALSE, sep = " ", quote = "",as.is = TRUE,skip=2) for (i in 3:9) { annot[[i]] <- gsub(",",";",annot[[i]]) } save(annot,file="annotation_chr1.Rda") sPF <- samplesPerFeature(X="ALL.chr1.merged_beagle_mach.20101123.snps_indels_svs.genotypes_mat",samples=EUeurop,lowerB=1.5,upperB=0.05*length(EUeurop)) save(sPF,file="sPF_15_005_chr1.Rda")