search for: genabel

Dear all, since yesterday, I have been experiencing problems with the package GenABEL. When I try to load the package (library(GenABEL)) I get the following error message: Loading required package: MASS Error : .onLoad failed in loadNamespace() for 'GenABEL', details: call: stringSplit[[1]] error: subscript out of bounds Error: package/namespace load failed for ?GenABEL...

Hi all! I am trying to install package - GenABEL in a Linux (RedHat) environment using R 2.6.2. The CRAN Mirror starts up, I select the site, but I end up with the following error: ****************************************************************** Warning: unable to access index for repository http://cran.stat.ucla.edu/src/contrib Warning messa...

Does anyone use the R library GenABEL? I am using it to calculate SNP interactions. I have a list of 100 SNPs, I need to look at the interaction between each of two SNPs among the list. my question is how to perform this in GenABEL. I want to use the "lm" function, but don't know how to use the SNP information. for exampl...

Hi all! I am using GenABEL on R for GWAS analysis. I am having a couple of issues: First, I am having a problem reading files (.map, & .ped, size 900Mb, using windows 32-bit) onto R in the "convert.snp.ped" statement. I am thinking this problem is likely due to the large size of the files & my version of...

Hi, I'm using kinship2 to calculate heritabilty, but I would like calculate in GenABEL too. I trying the code: > require(kinship2) > require(GenABEL) > pedig = with(Dados, pedigree(id=IID, dadid=PAT, momid=MAT, sex=SEX, famid=FID, missid=0)) > kmat = kinship(pedig) > (mod1 = polygenic(altura ~ SEX + idade, data=Dados, kin=kmat)) Erro em intI(i,...

Hi, I had originally posted regarding an error when trying to install package - GenABEL - it has now become clear that R is not able to connect to the net. Below are a couple of things I've tried with the resulting errors I am running R 2.6.2 on Linux e15 X86 How can I fix this? Thank you! 1) download.file("http://www.nytimes.com", destfile<-tempfile()) Error:...

Hi, I've been trying to use the impute2mach function in the GenABEL library. I unfortunately cannot make a working example as when I scale down the files I no longer get an error message. The error message that I get on the full file is: Error in checkSlotAssignment(object, name, value) : assignment of an object of class "NULL" is not valid for slo...

...# fixed number of permutations, 1000 is recommended alpha = 0.05 # type-I error Pval = 0.1 # significance level at first step of mbmdr set.seed(1982) # for reproducibility # loading libraries and code library(kinship2) library(GenABEL) source("mbmdr_gaussian.r") # removing all possible existing files first as a safety measure file.remove("permoutput.txt") file.remove("simulation.raw") file.remove("permuted.txt") file.remove("permchisq.txt") file.remove("rawoutput.txt")...

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 trailing tabs or spaces but it doesn't solve the prob...

Hi, I am using a function 'polygenic' in the package called 'GenABEL'. -------------------------------------------------------------------------------- > h2.gs <- polygenic(b.dbp~age+age2+age3+sex+b.bmi, kin=gs.gkin, data=gs) ------------------------------------------------------------------------------- ## and I got the following error message: Error in s...

I have a few questions about GenABEL, gwaa data. 1) is there a universal way that most GenABEL people use to add more individuals into a 'gwaa' data? For example, I have a 'gwaa' data, but I need to add some dummy parents, for 'gwaa at phdata', it's easy to add these rows, but for 'gwaa at gtdata',...

...plenty!) but I am still unable to solve the error. ERROR: cannot allocate vector of size 455 Mb I am using R 2.6.2 - x86_64 on a Linux x86_64 Redhat cluster system. When I log in, based on the specs I provide [qsub -I -X -l arch=x86_64] I am randomly assigned to a x86_64 node. I am using package GenABEL. My data (~ 650,000 SNPs, 3,000 people) loads in okay and I am able to look at the data using basic commands [nids, nsnps, names(phdata)] The problem occurs when I try to run the extended analysis: xs <- mlreg(GASurv(age,dm2)~sex,dta) ****************** 1) I have looked through the memory li...

I a using plink on a large SNP dataset with a .map and .ped file. I want to get some sort of file say a list of all the SNPs that plink is saying that I have. ANyideas on how to do this? -- Thanks, Jim. [[alternative HTML version deleted]]

Dear Sir/Madam, Since a few day now I try to use the command "polygenic" from the GenAbel package. However, I keep bumping up against an error message: "Error in polygenic(Testo, kin = kinship, data = data1) : dimension of outcome and kinship.matrix do not match". My data exists of 1240 individuals with 74 markers. It mainly consists of small families (2 or more brothers,...

Hi, Following my earlier posts about having problems performing a PCA, I have worked out what the problem is. The problem lies within the PLINK to gds conversion. It seems as though the SNPs are imported as "samples" and in turn, the samples are recognised as SNPs: >snpsgdsSummary("chr2L") Some values of snp.position are invalid (should be > 0)! Some values of

...kages reside on just one of the many Web sites devoted to R, and the number of packages has grown exponentially. One package, called BiodiversityR, offers a graphical interface aimed at making calculations of environmental trends easier. Another package, called Emu, analyzes speech patterns, while GenABEL is used to study the human genome. The financial services community has demonstrated a particular affinity for R; dozens of packages exist for derivatives analysis alone. "The great beauty of R is that you can modify it to do all sorts of things," said Hal Varian, chief economist at Goo...

...kages reside on just one of the many Web sites devoted to R, and the number of packages has grown exponentially. One package, called BiodiversityR, offers a graphical interface aimed at making calculations of environmental trends easier. Another package, called Emu, analyzes speech patterns, while GenABEL is used to study the human genome. The financial services community has demonstrated a particular affinity for R; dozens of packages exist for derivatives analysis alone. ?The great beauty of R is that you can modify it to do all sorts of things,? said Hal Varian, chief economist at Google. ?And...

I'm in the process of converting some S3 methods to S4 methods. I have this function : setGeneric("enrichmentCalc", function(rs, organism, seqLen, ...){standardGeneric("enrichmentCalc")}) setMethod("enrichmentCalc", c("GenomeDataList", "BSgenome"), function(rs, organism, seqLen, ...) { ... ... ... })

Hello, I want to know if there are some functions or packages to solve differential and integral equation using R. Thanks. Shao chunxuan. [[alternative HTML version deleted]]

...s/DatABEL a package providing interface to C++ FILEVECTOR library facilitating analysis of large (giga- to tera-bytes) matrices; matrix storage is organized in a way that either columns or rows are quickly accessible; primarily aimed to support genome-wide association analyzes e.g. using GenABEL and ProbABEL * DeducerPlugInExample (0.1) Ian Fellows http://crantastic.org/packages/DeducerPlugInExample A example GUI plug-in package to serve as a template. * DiceDesign (1.0) D. Dupuy http://crantastic.org/packages/DiceDesign Space-Filling Designs and Uniformity Criteria. * Dic...