similar to: Stats package

Hi further to my upgrade question thanks for pointing me in the right direction any ideas how I get round the following problem: I'm trying to install on a Tru64 alpha ecs1h[jps]69: make `Makedeps' is up to date. `libappl.a' is up to date. `Makedeps' is up to date. `libnmath.a' is up to date. `Makedeps' is up to date. `libunix.a' is up to date. `Makedeps' is up

Hi All Not sure if this a bioconductor question or general R mailing list so apologies if this has gone to the wrong one................. When plotting dendrograms created by hclust you can "identify" clusters by clicking on the graphics and returning a list of what is contained in each cluster. However I'd like to be able to "zoom in" on specific clusters and plot

HI All sorry if this question has already been asked but I couldn't find anything that answered my question I have 24 columns of data that I'm trying to plot in heatmap.2 (gregmisc) and I'm having difficulty ordering them except in numerical sequence: I have transposed my matrix so it will appear with the dendrogram I want appearing at the top of the heatmap If I use either of

Hi I know this subject has been mentioned before but from the mail archives I'm under the impression that this is not possible ? I'm trying to carryout the equivalent of ls -l in R to give some date/time label to each of my objects If the case is that there is no equivalent, is it possible to list all objects in an environment that share a common component ? So that the common

Hi, I would like some extra information on the 'cgh' package in R. I noticed that there isn't much activity regarding this package on the R and BioC mailing list (I googled it). I started using this package and I have few questions: 1/ As I have a custom tiling like array @8um features resolution (affy), I have a lot of probes to work with. I'm assuming it is correct to

Hi R friends, I've been attempting to create plots with multiple alpha values using Cairo to save them on a windows (32b XP) platform as it doesn't support more than 3 alpha values. This worked well until I wanted a postscript file (unsupported) and as a attempted work around I installed RGtk2. So far so good, however now when I try to use a >CairoPDF("alpha.pdf", 6, 6,

First, apologies for no example data but I don't think it's needed in this case, Q: Can (and if so how ) the line end style be changed for 'cloud' plots? I've tried par(lend=2), trellis.par.set(add.line = list(lend=2)) and much googling but to no avail Thanks in advance Dan P.S. the reason for this is that the round end looks bad at lwd=3 or more Daniel Alcock Malaria

Hi - I am interested in solving variance components for the data below with respect to the response variable, Expression within R. However, the covariates aren't independent and they also have a class (of which the total variance explained by covariates in that class I am most interested in). Very naively, I have tried to look at each individual covariates variance like this >

Hi I hope someone can shed some light on this: For some reason when I read.table("bfx.txt") R decides to only give back the first character from each column in each row as one single column. Like this: V1 1 ÿþr 2 \n 3 r 4 1 5 0 6 A 7 G 8 \n 9 r 10 1 11 0 12 T 13 C 14 \n The data should be:

Hello Please find bug report attached. Thanks Matthew -- Matthew Gillman Team 105: Variation Informatics Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK Tel: 01223 834244 Ext. 4922 -- The Wellcome Trust Sanger Institute is operated by Genome Research Limited, a charity registered in England with number 1021457 and a company registered

Dear All, Apologies for sending this email to both list, but at this point I'm not sure which one could help me the most. I have 4 sets of data, 1 test and 3 different sets of controls. The measurements are binary, with a matrix of 0 and 1 I'm measuring across time (rows, ~815) the behaviour of organelles in the cell by microscopy in response to different stimuli (several measurements

Hi, First my apologies for a non-working piece of code in a previous submission, I have corrected this error. I'm doing is individual based modelling of a pathogen and it's host. The way I've thought of doing this is with two dataframes, one of the pathogen and it's genes and effector genes, and one of the host and it's resistance genes. During the simulation, these things

Dear R users & Experts, This is just a curiousity, I was wondering why the dominant eigenvetor and eigenvalue of the following matrix is given as the third. I guess this could complicate automatic selection procedures. 0 0 0 0 0 5 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 Please

Dear all, I would like to know if it is possible to fit in R a Cox ph model with time-dependent covariates and to account for hierarchical effects at the same time. Additionally, I'd like also to know if it would be possible to perform any feature selection on this model fit. I have a data set that is composed by multiple marker measurements (and hundreds of covariates) at different time

Hi, I'm a sysadmin who's been tasked with installing R on our 1000-node compute cluster. I have licences for the Intel C and FORTRAN compilers, so I'm using the following to compile: CFLAGS="-O2 -axWK" FFLAGS=$CFLAGS CXXFLAGS=$CFLAGS CC=icc F77=ifort CXX=icc FPICFLAGS=-fpic ./configure --without-x --without-tcltk The compilation seems to go OK, with a few warnings.

Hi All, I'm trying to run a simulation of host-pathogen evolution based around individuals. What I need to have is a dataframe or table of some description - describing all the individuals of a pathogen population (so far I've implemented this as a matrix): ID No_of_Effectors Effectors (Sequences) [1,] 0001 3 ## 3

Hello, I would like to know if there is a function in R that will test for normality and handle censored data sets. Currently, I evaluate each censored data set by the extent to which a normal scores plot approximate a straight line. For complete data sets I use shapiro.test(). Below is an example of a censored data set. data1<-c(0.00, 0.00, 0.00, 5.86, 5.17, 8.17, 5.12, 4.92, 7.08,

Thorsten is right. There is a direct formula for computing the Moore-Penrose inverse using the singular value composition of a matrix. This is incorporated in the following: mpinv <- function(A, eps = 1e-13) { s <- svd(A) e <- s$d e[e > eps] <- 1/e[e > eps] return(s$v %*% diag(e) %*% t(s$u)) } Hope it helps. Dietrich

I am relatively new to R and am looking for advice, ideas or both... I have a data set that consists of pathogen population sizes on individual plant units in an experimental field plot. However, in order to estimate the pathogen population sizes I had to destroy the plant unit and could not determine if that plant unit became diseased or to what extent it would have become diseased. I

Hi Sanjiv, 2009/6/15 Sanjiv Gupta <sanjiv.gupta at microchip.com>: > Hi Mikhail, > How do you build mcc16 executable? This should work: $ cd $LLVM_DIR/tools/llvmc/examples/mcc16 $ make If you're building from some other dir, you'll need to update mcc16/Makefile, so it knows where Makefile.common is located. > There are so many confusing things there: driver, plugins,