Displaying 2 results from an estimated 2 matches for "phenofile".
2012 Aug 24
0
A question about GRAMMAR calculations in the FAM_MDR algorithm
...;simulation.raw")
file.remove("permuted.txt")
file.remove("permchisq.txt")
file.remove("rawoutput.txt")
# loading data and bringing in GenABEL format
rawfile="simulation.raw"
convert.snp.ped(pedfile, mapfile, rawfile)
simulation.GenABEL = load.gwaa.data(phenofile = phenofil, genofile =
rawfile, force=F,makemap=F,sort=F)
pedigree=read.table(pedfile)
pedsize=nrow(pedigree)
nsnps=(ncol(pedigree)-6)/2
# minor allele count and handling missing genotype data
allelic = function(k){
geno=pedigree[,(5+2*k):(6+2*k)]
allelic=rowSums(geno==2)-(geno[,1]==0 & geno[,...
2010 Feb 24
6
R error- "more columns than column names"
Hi all! I am desperately trying to figure out the solution to this error, but
nothing as of yet is working.
As noted in an earlier post I am using GenABEL. In an attempt to read in
the phenotype file, in the format .dat, R keeps giving me the error "more
columns than column names"
I have tried to read in the data without the headers; I have also tried to
trim the data to remove any