search for: unpaired

...nd future code I write run more quickly/efficiently, I will be most grateful, as as this is currently written it would take what looks like days to run a thousand simulations of each possible combination of variables of interest. Best Ben Caldwell ----------------------------------------------- #unpaired verification.plots<-rnorm(30, 100, 10) # writeClipboard(as.character(verification.plots)) unpaired.test<- function(verification.plots, project.n, project.mean, project.sd, allowed.deviation, project.acres, alpha=.05){ verification.plots <-as.numeric(as.character(verification.plots)) a &l...

Please forgive a stats question. I have to sets of data (unpaired) measured on an ordinal scale. I want to test to see if the two sets are different (i.e. do they have the same location): set1: 1,3,2,2,4,3,3,2,2 set: 4,4,4,3,3,5,4,4 What is the most appropriate non-parametric test to test location? Thanks, John Confidentiality Statement: This email messa...

...eciate it. Below is what I tried and what happened (with the error highlighted in bold red font). thanks in advance, Sheila > library (samr) Loading required package: impute > > ## - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - > ## testing the two class unpaired example from the man pages ... > > set.seed(100) > x <- matrix ( rnorm(1000*20), ncol=20 ) > dd <- sample ( 1:1000, size=100 ) > u <- matrix ( 2*rnorm(100), ncol=10, nrow=100 ) > x[dd,11:20] <- x[dd,11:20] + u > y <- c(rep(1,10),rep(2,10)) > > data <- li...

Dear, I wanna to compare AUC generated by two distribution models using the same sample. The AUC for model 1 consists of two columns, column A for 0/1 and column B for probability, eahc with the same row number of 3000. The AUC for model 2 consists of two columns, column A for 0/1 and column B for probability, eahc with the same row number of 10000 rows. I am wondering what value I should put

Dear, I wanna to compare AUC generated by two distribution models using the same sample. The AUC for model 1 consists of two columns, column A for 0/1 and column B for probability, eahc with the same row number of 3000. The AUC for model 2 consists of two columns, column A for 0/1 and column B for probability, eahc with the same row number of 10000 rows. I am wondering what value I should put

...rent from the other observations (ie x[-1]). My null hypothesis is that x[1] is drawn from a Gaussian distribution of the same mean as observations x[-1], which are assumed to be iid Gaussian. The (unknown) variance of x[1] is assumed to be the same as the variance of x[-1]. This should be an unpaired t-test. But > x <- c(23,25,29,27,30,30) > t.test(x=x[1] , y=x[-1]) Error in t.test.default(x = x[1], y = x[-1]) : not enough 'x' observations > What arguments do I need to send to t.test() to test my null? -- Robin Hankin Uncertainty Analyst National Oce...

...int with 0.25 inches warning: no point with 0.25 inches Why does this happen. 2) Now I am trying to run EBAM: and when I type I get an error > find.out<-find.a0(datrmanhi,groupsnhi.cl,rand=123) Loading required package: affy Loading required package: affyio EBAM Analysis for the two class unpaired case. Warning: There are 1 genes which have variance Zero or no non-missing values. The d-value of these genes is set to NA. The following object(s) are masked _by_ .GlobalEnv : n The following object(s) are masked from mat.repeat ( position 5 ) : c...

...culated using R2.1.1. Even the example code doesn't work the same way under these two versions of R. Does anybody know what is going on? Thanks for any suggestions. library(multtest) data(golub) # golub.cl contains the class labels. golub.cl # Perform a SAM analysis for the two class unpaired case assuming # unequal variances. sam.out<-sam(golub,golub.cl,B=100,rand=123) sam.out # Obtain the Delta plots for the default set of Deltas plot(sam.out) Best Wei Zhao Postdoctoral Research Fellow Department of Human Genetics and Biostatistics University of California, Los Angele...

What this message means is that a "{" showed up when some other bracket was unpaired. In this case, if you check your code, you'll see that "if(MatriceDist[i,j] > 0 & ((vectorID[i] > 0 | vectorID[j] > 0))" is lacking a closing ")" for the if clause. Phalaen wrote > Hi! > I am a Phd student in the university of Padua and I am tryin...

...t; 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 False Called FDR 1 0.1 0.929 292.25 293 0.927 2 0.4 0.929 43.60 56 0.724 3 0.7 0.929 12.25 20 0.569 4 1.0 0.929 7.25 14 0.481 5 1.3 0.929 2.60 7 0.345 6 1.7 0.929 1.30 5 0.242 7 2.0 0.929 1....

...tal about 20 questions. I have to say my statistic skills are very basic and very rusty, we had some lectures two years ago, where we were introduced to R. I looked through my notes, and back then we did a one sample t-test to analyse likert type questions. I believe I would need to do a 2 sample unpaired t-test. It would be great if someone could give me some feedback if this test is the most suitable one for my purpose, and maybe could explain to me what?s the easiest way to do this in R? You would help me loads!! Many thanks in advance Mona -- View this message in context: http://n4.nabble.c...

Some of my students were interested in a plot where the sizes of the symbols were proportional to the number of repeats in the data (I have them plotting discrete data, so there are often overlaps). Here's the solution I came up with, which seems clever (to me) but not necessarily the best way. Can anyone suggest improvements? sizeplot <- function(x,y,maxsize=0.2,

...int with 0.25 inches warning: no point with 0.25 inches Why does this happen. 2) Now I am trying to run EBAM: and when I type I get an error > find.out<-find.a0(datrmanhi,groupsnhi.cl,rand=123) Loading required package: affy Loading required package: affyio EBAM Analysis for the two class unpaired case. Warning: There are 1 genes which have variance Zero or no non-missing values. The d-value of these genes is set to NA. The following object(s) are masked _by_ .GlobalEnv : n The following object(s) are masked from mat.repeat ( position 5 ) : c...

I'm running samr (Two class unpaired), but keep getting the following error: perm= 1 Error in if (logged2) { : argument is of length zero <code> library (impute) library (samr) data = list (x=dat, y=y, geneid = matrix(twoUnpaired.data[,1],ncol=1), genenames = matrix(twoUnpaired.data[,2], ncol=1)) samr.obj <- samr (dat...

Hi, I am not sure which kind of test is applied to the data if you use cor.test(x, y) ? Is it an unpaired t-Test? Regards -- View this message in context: http://www.nabble.com/cor.test%28x%2Cy%29-tp22492993p22492993.html Sent from the R help mailing list archive at Nabble.com.

Does anyone know of an R-function or method to compare two C-statistics (Harrells's C - rcorr.cens) obtained from 2 different models in partially paired datasets (i.e. some similar and some different cases), with one continuous independent variable in each separate model? (in a survival analysis context)? I have noticed that the rcorrp.cens function can be used for paired data. Thanks

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" &...

Hi! I'm attempting to start setup.exe from office97 package. The installer starts but then bails out with the dialog Wrong parameters on the command line: " (see the unpaired double quote) and then pops another dialog with its usage, then exits. I'm calling it just with wine setup.exe, no quotes nor other characters there. I was trying to compensate this character by putting one more from my command line but I was unable to do it (for example, calling wine setup.e...

...the data, so I'd rather use distribution-free methods. I would like to have statistical justification for the statement that the snow depth in the second period is less than in the first half, and that the variability decreased as well. For the difference in central measures, I am using the (unpaired) wilcox.test, but I really have no idea how to address the question of changes in variability using nonparametric means. Any ideas? Thanks, Tony

Hello list, I am trying to apply the paired t.test between diseased and not diseased patients to identify genes that are more expressed in the one situation under the other. In order to retrieve the genes that are more expressed in the positive disease state I do: p.values<-c() for(i in 1:length(Significant[,1])){ p.values[i]<-try(t.test(positive[i,],negative[i,],alternative