similar to: Using PCA to correct p-values from snpMatrix

The windows replacement is R-library-grDevices-libs-winCairo_20150122.zip in http://sourceforge.net/projects/outmodedbonsai/files/R/ The cairo replacements are just two bundles, windows or OS X - containing replacements for R 2.15.3, 3.0.3 and 3.1.2. You just extract the relevant one for your R version. The windows builds of R packages snpStats x.x.x.8 and snpMatrix 1.19.0.20 were done with

R.framework-Versions-Resources-library-grDevices-libs-cairo_20150120.tgz in http://sourceforge.net/projects/outmodedbonsai/files/R/ are dropped in replacement to the cairo.so's in the official R binaries (2.15.3, 3.0.3, 3.1.2). updated to cairo-1.12.18 and freetype-2.5.4. The official R binaries' were built with early freetype 2.4.x and cairo 1.11(?) and had a number of issues with some

Dear R helpers I have a table and i need to make new table table1: sire snp1 snp2 snp3 snp4 snp5 snp6 snp7 snp8 snp9 snp10 snp11 snp12 snp13 snp14 snp15 8877 -1 -1 -1 -1 0 0 -1 -1 -1 0 1 1 1 -1 -1 7765 1 1 1 0 0 0 -1 1 1 1 0 0 0 1 0 8766 1 1 -1 0 -1 -1 0 -1 0 -1 -1 -1 0 1 0 6756 0 1 0 -1 1 -1 -1 0 0 0 0 -1 0 1 1 5644 -1 0 1 -1 0 0 0 0 -1 -1 0 0 0 0 1 I have table2 sire

Hi, I'm looking for an efficient code that will enable me to reshuffle data (phenotype) for certain number of individuals and creating a loop that will randomly simulate it for 10000 times *(permutation)*. I also need to find how I keep the information (p value for each SNP) gathered for all the 10000 iterations. My data set looks like this (n=500): Individual # Phenotype SNP1 SNP2

Hi Mr Tierney, I have noticed an error message from R 1.12.x's CMD check for a while (apparently prof Ripley completely rewrote CMD check in R 1.12+) e.g.: http://bioconductor.org/checkResults/2.7/bioc-LATEST/snpMatrix/lamb2-checksrc.html ---------------- * checking R code for possible problems ... NOTE Warning: non-unique value when setting 'row.names': ?new? Error in

FYI. There is a Thai vignette - and it went a lot further doing some Thai text processing in R, than the earlier Chinese/Tibetan/LiangshanYi/Arabic vignette, which was in reality just Chinese + Cairo graphics. Managed to cross-compile an R package for Mac OS X from Linux; and it seems to be working. See screenshots below. I'd be interested to know if there are less obvious bugs; however

This is somewhat a summary/continuation of an R bug report: (https://bugs.r-project.org/bugzilla3/show_bug.cgi?id=14645) Illumina's cluster definition files (*.egt) are one of the proprietary and undocumented file formats used by their GenomeStudio line of products for genomic studies. snpMatrix 1.17.0.7 onwards (http://sourceforge.net/projects/outmodedbonsai/files/snpMatrix%20next/)

Current core/Recommended Matrix package (0.999375-48) has been segfaulting against R 2.13-alpha/2.14-trunk for the last week or so (since R-2.13 was branched, when I started trying) when "run with R CMD check --use-gct": -------------- > pkgname <- "Matrix" > source(file.path(R.home("share"), "R", "examples-header.R")) >

I finally got round to look at a little problem I have - most of the time when I use R I would be using snpMatrix (now part of bioconductor, I wrote a substantial part of), so I had $HOME/.Rprofile to save some typing. Upgrading, switching versions used to work fine with R 2.8 then it broke with 2.9.. Apparently it is because most of SHLIB, INSTALL, etc used to be fairly extensive shell scripts

(e-mailing to R-bugs is intentional - the web itnerface seems to be down) > a<- cbind(c(1,2), c(3,4)) > a [,1] [,2] [1,] 1 3 [2,] 2 4 > a[cbind(c(2,2), c(2,1))] [1] 4 2 > a[cbind(c(2,3), c(2,1))] Error: subscript out of bounds > a[cbind(c(2,-1), c(2,1))] Error: negative values are not allowed in a matrix subscript > a[cbind(c(2,0), c(2,1))] [1] 4 Am

I have had some fun in the last few days trying to put together an annotated map of China with R and some public GIS data: http://sourceforge.net/projects/outmodedbonsai/files/snpMatrix%20next/1.17.7.11/China_Choropleth_Maps.pdf/download It is done, and rather nice... there are a few issues: - the default pdf() device cannot do CJK with embedded fonts - and cairo_pdf() is not hooked up to

Hi, I have got a curious issue with an R package which uses zlib, against the official binary here: http://cran.r-project.org/bin/linux/redhat/el4/i386/R-2.4.1-1.rh4AS.i386.rpm on a Redhat EL4 i686 system. The problem is that at the end of reading a gzip'ed file within my C code, gzgetc() returns -1 (no more to read or error) but gzeof() doesn't return true. Now one can probably jump to

snpMatrix package is quite nice (read.plink())

It was as easy as I thought it was half a day ago - here is a patch against R trunk to add cairo support to the Sweave driver, an example Sweave input, and the resulting output. A few more notes: - obviously the documentation needs to be updated... a bit more work to do. - some check to make sure "cairo" and "pdf" are not both set would be nice, as well as checking

Hi, (somebody would probably yell at me for not checking 2.6.0rc, for which I can only apologize...) Our R package (snpMatrix in http://www-gene.cimr.cam.ac.uk/clayton/software/) is broken rather badly in 2.6.0 ; I have fixed most of it now so a new release is imminent; but I'd like to mention a few things, mostly to summarize my experience and hopefully the 'writing R extensions'

--- On Thu, 3/9/09, Vinh Nguyen <vinhdizzo at gmail.com> wrote: > hi hin-tak, > > i'm trying to build r packages for windows on a > mac/linux. i guess > this used to possible and supported, but is no longer > supported. i > ran into this post of yours, > https://stat.ethz.ch/pipermail/r-devel/2009-July/053971.html, > and hope > u don't mind me

Dear All I apologise for such a simplistic question. However I am having immense trouble getting my S-Plus data set into R. I have attempted data.dump/data.restore though I am unable to get it to work for me. The error I am getting in R is : Error in seek(s) : no applicable method for "seek" Is this the best way to import S-Plus datasets or is there another suggestion that I could

Hi, I am trying to run PCA on a set of data with dimension 115*300,000. The columns represnt the snps and the row represent the individuals. so this is what i did. #load the data code<-read.table("code.txt", sep='\t', header=F, nrows=300000) # do PCA # pr<-prcomp(code, retx=T, center=T) I am getting the following error message "Error: cannot allocate vector of

Hello, I did a pca with over 200000 snps for 340 observations (ids). If I plot the eigenvectors (called rotation in prcomp) 2,3 and 4 (e.g. plot (rotation[,2]) I see a strange "column" in my data (see attachment). I suggest it is an artefact (but of what?). Suggestion: I used prcomp this way: prcomp (mat), where mat is a matrix with the column means already substracted followed by a

Dear, I would like to deal with microarray data, it can run when i deal with little data. However, the amount number of SNP data are 45181, amount numbers of animal are 3081,it can not be allocated 1000Mb memory when i importing them to R Procedure sentence show: m<-matrix(scan("D:/SNPdata.txt"),ncol=nmarkers,byrow=TRUE) Error show: Error: cannot allocate vector of size 1000.0 Mb