similar to: samr error

I am trying to use a new (to me) package (samr) and even when I try to run a very simple example, I get this "cannot open the connection" error. The reason I am writing to r-help rather than to the authors of samr is I think this may be a more general R problem rather than a samr-specific problem. Perhaps something with my installation and write access to some particular place ? I am

Hi BioC user, I have a problem extracting the gene set I would like to work with. Here is I work with my data: normData <- read.delim("normalizedData.txt",sep ="\t") ######### two class unpaired comparison # y must take values 1,2 classes <- c(-1,-2,1,2) #prepere the data for the samr analysis data.x <-as.matrix(normData[,8:11]) d=list(x=data.x,y=classes,

Hello Everybody, I am trying to perform SAM with the samr package. I am using the following code: sink ("R005") library(siggenes) library(samr) library(nnet) A <- as.matrix(read.table("D:\samrgenes1000.txt")) B <- as.matrix(read.table("D:\genenames1000.txt")) y1 <- c(rep(1,20),rep(2,6)) #there are 20 chips of one kind and 6 of the other kind. datasam =

Hi R-help, I am trying to use 'samr' for 10 pre and post paired samples to test whether post is different from pre (i.e., the location shift for the delta of (post-pre)). However, I got an error message saying > samr.obj<-samr(d, resp.type="Two class paired", nperms=100, random.seed=100) perm= 1 Error in !logged2 : invalid argument type Does anyone know what this

Hi, I was working on a classification problem using the pamr package. I used the pamr.adaptthresh() function to find the optimal accuracy of the classifier. I must not be doing it right, since it doesn't return the threshold values for optimum classification. For example,if I run it on a dataset, I get the following result using pamr.adaptthresh(): predicted true

Hi, I was working on a classification problem using the pamr package. I used the pamr.adaptthresh() function to find the optimal accuracy of the classifier. I must not be doing it right, since it doesn't return the threshold values for optimum classification. For example,if I run it on a dataset, I get the following result using pamr.adaptthresh(): predicted true (1)

Hi, Sorry about the earlier formatting errors... I was working on a classification problem using the pamr package. I used the pamr.adaptthresh() function to find the optimal accuracy of the classifier. I must not be doing it right, since it doesn't return the threshold values for optimum classification. For example,if I run it on a dataset, I get the following result using

Hello all! I am currently trying to sort a data frame in a particular way, but I am having some difficulties with this. Specifically I want to sort the below dataset in such a way that there is only one line per ProteinID and if there are multiple GeneID or GeneName entries for a single proteinID, that they be concatenated with a comma separating them. The way I have done it earlier worked fine

Dear list, I have gene expression measurements obtained by PCR on 11 genes, tabulated as a data matrix. I'm attempting to use GSA package to distinguish any significant changes in these genes as a pathway. My response variable is binary, 0=no disease, 1=disease. I have read the PCR data into R as follows: data <-

Dear all, I am running R 2.4.1. > library(siggenes); > library(multtest); > cl<-rep(c(0,1),c(3,3)); > sub<-exprs(AffyExpData[,c(1:3,7:9)]); > gn<-geneNames(AffyRAwData); > sam.out<-sam(sub,cl,rand=123,gene.names=gn); We're doing 20 complete permutations > sam.out SAM Analysis for the Two-Class Unpaired Case Assuming Unequal Variances Delta p0

Dear R-Helpers In R 2.5.1, the command library(Rgraphviz) fails on my Windows (XP SP2) system with error popup "The procedure entry point Rf_allocString could not be located in the dynamic link library R.dll". Thanks in advance for any suggestion in solving the error. My D. Coyne Imagination is more important than knowledge... (Albert Einstein) mcoyne@boninc.com

Hi All, I'm new with R; this is a basic question. I was given a matrix I of (nrow, ncol), I would like to create another matrix A with some data in the matrix I, say [1,4] (row 1, column 4) to [271,19000] (row 271, column 19000). How do I do this? Please help. Thank you very much. --mc [[alternative HTML version deleted]]

Hello, I received the following warning when running chi-square; n Is there a way to catch the 'error' code of 'warning' after run chisq.test(x)? n What does this error mean? Thank you for your help. [[alternative HTML version deleted]]

I don't know if this is the correct forum to ask the following question; however, when I search the aroma.affymetrix discussion group, it suggested that I should posted the question to r-help. Here it goes. I followed the instructions on aroma.affymetrix trying to install the packages; following are the steps: > install.packages(c("R.oo", "R.utils"),

Hi, here's my siggenes.table$genes.up snippet. Two class unpaired SAMR analysis. "Row" "Gene ID" "Gene Name" "Score(d)" "Numerator(r)" "Denominator(s+s0)" "Fold Change" "q-value(%)" "1" "25" "RPL15P22" "RPL15P22" "-1.44115338424578" "-18"

Hi, Can you please tell me what does the the function logged2 in R do im list or..? As I have > ?logged2 No documentation for 'logged2' in specified packages and libraries: you could try '??logged2' > ??logged2 No help files found with alias or concept or title matching ‘logged2’ using fuzzy matching. Thank you in advance, Amor [[alternative HTML version

Hello, Newbie question and hope you can help . I have two vector V1 and V2, where length(V2) = length of (V1) * 2; length(V1) ~ 16,000. For each member in V1, I need to compare 2 element of V2 for equality i.e. for (I in 1:length (V1)) { if ( v2[i] == v1[i] & v2[i+1]==v1[i] ){ statement_1 statement_2 . } } This for-loop is too slow

Hi, I place a t.test in a loop and would like to continue to process the loop even when t.test encounter error. How do I do that? For example, in one iteration, the data is completely constant and t.test gives error, the entire program terminates. I would like to write the information out to a file, and the loop should continue. Thanks My D. Coyne [[alternative HTML version

R-newbie question I have 2 matrices (a) P1 has only one column of 32K rows (b) PC has 2 column {P, C} of 3200 rows Every values in P1 matches with a value in PC[,p] (column p). I would like to use Which to search for all value in P1 that matchex PC[,p] and get the PC[,c]. However because P1 and PC does not have the same number of rows, I got lots of 'NA'. Thanks for your

Running R2.4.0 on Apple Mac OS X 10.4.8, in Emacs ESS mode, and also R.app. In an attempt to learn a bit more about a particular method (geneNames in package affy) I invoked getMethods("geneNames") which produced geneNames methods, but not the one in affy (output below). I had to know the signature (AffyBatch) in order to find the method > getMethod("geneNames",