similar to: New package: noia

Dear all, I am running qtl mapping. I have 75 RI lines with some residual heterogeneous loci. The loci are code A, B or H(heterogeneous). Questions: 1) R/qtl determine the data is F2 intercross. 2) Warning message about strange genotype pattern > library(qtl) > dat=read.cross("csv", file="rqtl_trt.csv") --Read the following data: 75 individuals

Hello R Forum users, I was hoping someone could help me with the following problem. Consider the following "toy" dataset: Accession SNP_CRY2 SNP_FLC Phenotype 1 NA A 0.783143079 2 BQ A 0.881714811 3 BQ A 0.886619488 4 AQ B 0.416893034 5 AQ B 0.621392903 6 AS B 0.031719125 7 AS NA 0.652375037 "Accession"

Dear experts, I am trying to perform an association using snpStats. I have a snp matrix called 'plink' which contains my genotype data (as a list of $genotypes, $map, $fam), and a phenotype data frame which contains the outcomes (outcome1, outcome2,...) I would like to associate with the genotype. My question is, how do I loop through the outcomes? This type of data seems different from

Dear all,   I have data from the following experimental design and trying to fit a mixed model with lme function according to following steps but struggling. Any help is deeply appreciated.   1) Experimental design: I have 40 plants each of which has 4 clones. Each clone planted to one of 4 blocks. Phenotypes were collected from each clone for 3 consecutive years. I have genotypes of plants. I

Dear UseR, I get an error when I run "R CMD check" on my .tar.gz file package, and I don't understand why since I don't obtain any error with "R CMD check" on the package directory. Do you have any idea ? $ sudo ./R-2.11.0/bin/R CMD check eqtl_1.1.tar.gz and $ sudo ./R-2.11.0/bin/R CMD --check-subdirs=no eqtl_1.1.tar.gz return an Error * checking for working

Dear R-experts, My problem is how to handle a 10GB data file containing genotype data. The file is in a particular format (Illumina final report) and needs to be altered and merged with phenotype data for further analysis. PERL seems to be an frequently used solution for this type of work, however I am inclined to think it should be doable with R. How do I open a text-file, line by line,

Hi, I'm using R for a QTL analysis of SNP data. I was wondering if anyone had any advice on fitting a dominance effect into the following function; > myfun4 function (x) { x <- scan(con, nmax=169) y <- unique(x[which(!is.na(x))]) if(length(y)>1) { summary(lme(Ad ~ x, random= ~1|sire, na.action="na.omit")) } else {print("no.infomation")} } Con is the

I want to specify a two-factor model in lme, which should be easy? Here's what I have: factor 1 - treatment FIXED (two levels) factor 2 - genotype RANDOM (160 genotypes in total) I need a model that tells me whether the treatment, genotype and interaction terms are significant. I have been reading 'Mixed effects models in S' but in all examples the random factor is not in the main

Hi, I am conducting an association analysis of genotype and a phenotype such as cholesterol level as an outcome and the genotype as a regressor using multiple linear regression. There are 3 possibilities for the genotype AA, AG, GG. There are 5 people with the AA genotype, 100 with the AG genotype and 900 with the GG genotype. I coded GG genotype as 1, AG as 2 and AA as 3 and the p-value for the

Hi, My query is regarding permutation test and reshuffling of genotype/phenotype data I have been using the haplo.stats package of R. for haplotype analysis and I would like to perform an analysis which I'm requesting your advice. I have a data set of individuals genotyped for 12 SNP and a dichotomous phenotype. At first, I have tested each of those SNP independently in order to bypass

Hello, I have data in a dataframe with 139104 rows which is multiple of 96x1449. i have a phenotype file which contains the phenotype information for the 96 samples. the snp name is repeated 1449X96 samples. I haveto merge the two dataframes based on sid and sen. this is how my two dataframes look like dat<-data.frame(snpname=rep(letters[1:12],12),sid=rep(1:12,each=12),

Dear Mme/Mr. Hope you are doing well. I am doing some genetic analysis using The R software and I have difficulties to find how I can perform an Interaction/epistasis analysis using 3 or more SNPs (=markers) ? (In the instructive manual, there is only an interaction/epistasis analysis with 2 markers). In addition can you please inform me how I can perform Haplotype analysis and if there is an

Could you help me in getting a package to compile from source in Windows? I'm running Windows NT 4.0 on a pentium 450 with 128 Mb of Ram. I've got R 1.0.0. I am developing a package for R which includes a bit of C code. I have things working in Unix just fine, and now I want to port it to Windows NT. I've been following the instructions in the R complements to Venables and

Dear all, we seek a talented and motivated post-doc to develop computational methods for inferring the molecular basis of genetic diseases by integration of personal omics data. Research topics include: identifying causal mutations of rare disease patients by meta-analysis; inferring disease-causing molecular pathways from genotype, phenotypes, and omics profile of patient-derived cell lines

Dear R-Helpers, I am using qtl package to analyze qtl data from QTL cartographer. I have the map file and cro file from QTL cartographer. I was trying to import these two files in R using qtl package. data=read.cross("qtlcart", ".", "crofile.txt", "mapfile.txt") ### I have matched the file structure with the one on the website of qtl package - It matches

On Mon, 8 Nov 2004 ssim at lic.co.nz wrote: > > Is there an add-on package in R for QTL interval mapping for outbred > population, eg. Haley-Knott regression method ? > > Stella I believe Karl Broman's R/QTL (http://biosun01.biostat.jhsph.edu/~kbroman/software/) and Brian Yandell's bim (http://www.stat.wisc.edu/~yandell/qtl/software/), Richard Mott's happyR

Hi R, I tried to install one local source package---"qtl" into R, but failed with the following message. I am using Macbook OS 10.8.3 with Xcode 4.6.1. ld: warning: directory not found for option '-L/usr/local/lib' ld: library not found for -lgfortran collect2: ld returned 1 exit status make: *** [qtl.so] Error 1 ERROR: compilation failed for package ¡®qtl¡¯ * removing

Hello R-users! I have a style question. I know that for loops are somewhat frowned upon in R, and I was trying to figure out a nice way to do something without using loops, but figured that i could get it done quickly using them. I am now looking to see what kind of tricks I can use to make this code a bit more aesthetically appealing to other R users (and learn something about R along the

Hello... I tried to compile an R addon package, R/qtl on a Solaris machine and I get this error. I know that R/qtl complies correctly on other platforms (Windows, Mac, Linuxes, etc), so I am wondering what the problem might be. Any ideas? Thanks, Saunak ----------------------------------- /export/home/sen/tmp <howrahbridge> <pts/5> R INSTALL

Hi, If you have multiple players installed on your PC/Mac .m3u will always open the last media player who are the default in charge of the extension and mime m3u. On your web site you want maybe to force a link to open real media player or quicktime/itune. You need to create a .pls to force winamp loading the streaming because windows media player won't open .pls etc. If you add a .pls