similar to: Problem indexing into array

In my package, I create a new method for plot with the following signature: setMethod("plot", signature(x="marrayNorm", y="formula"), plot.ma) where marrayNorm is a class defined in the marray package. After building and installing my package, I get the following warnings when I load my package (with library(maVis)): Warning messages: 1: In the method signature

Hello, For a simple problem with 1 predictor (x) and 2 classes (0 and 1), the linear discriminant function should be something like 2(mu_0 - mu_1)/var x + x-independent-terms where var is the common variance.

Hello, Is there a built-in test in R for hypothesis testing with samples of known variance? For example, I've got a set of data, a mean to compare against, and a known variance, and I want to determine the p-value for which I can reject the null hypothesis (mu_1 = mu_0) and accept the alternative (mu_1 > mu_0). I've found that JMP and Minitab can both do this (in JMP, it's a

Dear R users, I'am using marray and Limma packages to analyze genepix output. 1) how can I filter bad spots from my data (data contains 3 types of bad spots). my experiment contains 12 samples and the bad spot are not associated to the same probes 2) how can I remove control probes from my data ? I'm sorry, i'm new with R and I can't find answer in packages doc. best regards,

Dear all, Trying to produce 4 maImage plots (marray package) on the same device (2 on the top and 2 on the bottom) with the layout() function or the split.screen() function, we are facing the following problem: it seems that maImage() does nt care about any of these 2 functions, and plots only one image at a time. Maybe this is inherent to this maImage() function, but we did not find anything

To whom it may concern: I have a question about how to appropriately conduct an lmer analysis for negative binomially distributed data. I am using R 2.2.1 on a windows machine. I am trying to conduct an analysis using lmer (for non-normally distributed data and both random and fixed effects) for negative binomially distributed data. To do this, I have been using maximum likelihood,

Neither biocLite nor the GUI menus can install packages on my system. Here is relevant output: > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 11.1 year 2010 month 05 day 31 svn rev 52157 language R version.string R version 2.11.1 (2010-05-31) > source("http://bioconductor.org/biocLite.R") BioC_mirror =

I'm working on some heatmaps, and the person I'm working with would prefer a red-black-green color palette (red denoting gene induction and green denoting gene repression). Does such a palette exist already? If not, is there an easy way to create one? Thanks, Jake

Dear all, I am finding difficulty in the following, I would like to create an empty matrix e.g. 10x10 of 0s and sequentially fill this matrix with randomly placed a 1s until it is saturated. Producing 100 matrices of sequentially increasing density., This process needs to be randomized 1000 times., I assume i should run this along the following lines, 1) Create 1000 matrices all zeros, 2) add

Hi you guys: I have been wondering if there is any way to change the labeling in plot.ts( ), for example , if I plot two sequences, i always got y labels as "series1", "series2", I tried to use ylab=c((expression(mu_1)),(expression(mu_2))), but it doesn't work. thanks in advance for any help best.

Dear R People: I have a few questions about multiple comparisons and contrasts for ANOVA, please. I've tried some things but with no success. Suppose I have a completely randomized design, and I want to have the contrast \mu_1 - 0.5 \mu_2 - 0.5 \mu_2 How do I set that up, please? I used the C command, and ran aov, but the results were identical to those with no contrasts. Also, is there

Dear mailing list, may some one be kind to help me solve following problem. I am trying to write a code that will combine two tables "x" and "y". The first columns of both tables are unique identification for the rows. The first column of table "X" is a sub set of the first column of "Y". I need to find the matching rows in both tables by looking on their

Hello, I am interested in performing a 2D loess smooth on microarray data, i.e. log2 ratios on a 2D grid, using different spans in the horizontal and vertical directions (the immediate reason being that replicate spots are laid out in the horizontal direction). Is it possible to do this in R? Functions like loess(stats) seem to apply the same span for all predictors, which carries over to

Hi all, is it possible to modify the way a graph obtained through image(graphics) is filled, I mean starting filling the graphical matrix by row from the upper-left rather than by the lower-left cell... In many cases, it can be usefull to have a representation of the data spatialy corresponding to a real support, as it is the case with the function image(marray) from Bioconductor packages, which

Hi, I am wirting a R script to process some data routinely and have a problem when try to get the content of object in ls()? Here is the script: > ls() [1] "a" "b" [3] "files.num" "fll92_1a_gpr" [5] "fll92_1a_gpr_norm" "fll92_1a.gpr.norm.scale" [7] "fname" "fullnames" [9] "gal" "galfile" [11]

Hello R-Users! I'm using a heatmap to visualize a matrix of values between -1 and 3. How can I set the colors so that white is zero, below zero is blue of increasing intensity towards -1 and above zero is red of increasing intensity towards red? I tried like this (using the marray and gplots packages from bioconductor): mcol <- maPalette(low="blue", mid="white",

Hi, I am trying to install Bioconductor onto R version 2.7.0 for Windows. I installed R, then followed the instructions on http://www.bioconductor.org/download, which state that you should type the following: source("http://bioconductor.org/biocLite.R") biocLite() When I do that, I get the following error: Running biocinstall version 2.2.9 with R version 2.7.0 Your version of R

The inst/doc directory of the DAAG package has 6 files coral551.spot, ... that are around 0.85 MB each. It would be useful to be able to zip then, but that as matters stand interferes with the use of the Sweave file that uses them to demonstrate input of expression array data that is in the "spot" format. They do not automatically get unzipped when required. I have checked that

Hi, How can I write the output to an excel (csv) file without printing row names (i.e without breaks). Here is my code: library( fn <- function() { q <- c(1,2,3) write.csv(q,"C:/Temp/op.xls", append = TRUE, row.names = FALSE,quote = FALSE) } # Function Call for(i in 1:3) { fn() } Present Output : x 1 2 3 x 1 2

Hello, I have some parameters from Mclust function. The parameters are in the form *parametersDf * * mu_1 mu_2 var_mc1 var_mc2 c1 c2 * *2 1.357283 2.962736 0.466154 0.1320129 0.5258975 0.4741025 * *21 8.357283 9.962736 0.466154 0.1320129 0.5258975 0.4741025 * Each row in the above data frame