similar to: Significance test for discrete data

Dear everyone, I am trying to do a Wilcoxon one-sided test on my gene expression data. These are the data i have in R: data.matrix (matrix, numeric) containing all gene expression data (42 rows=genes, 42 columns=tumors), no column header or row names data.cl (vector, numeric) consisting of 42 0's and 1's to indicate class 0 or class 1 for each column in data.matrix I want to do a

Hello everyone, I am employing the waba method from the multilevel package for obtaining a within-group correlation (Description: http://bg9.imslab.co.jp/Rhelp/R-2.4.0/src/library/multilevel/man/waba.html). Does anybody know a way or a calculation for obtaining a significance value for that correlation? And another question: Does anybody know whether it is possible to save individual

Hello, we are trying to use samba 3.6 over a big network with multiple xdsl link (512Kbit to 10 Mbits not symetric). We have a main site with two domain controlers (Windows 2k8 R2) and about one hundred middle sites with samba member servers. We have a big issue with Active Directory (AD) ACL which disapear when member server (powered by samba and winbind) can't discuss (for network

Hello all, I have the square symetric matrix A: 2 1 1 1 2 1 1 1 2 My first question is what is the easiest way to enter this matriz in R? Second, matrix A has an eigenvalue with multiplicity 2, in this case, how could I find the two related ortogonal eigenvectors given below by R, without the help of R, I mean, I want to know how R calculate this eigenvectors related to the same eigenvalue.

Hello, I am inputing a 17 x 17 symetric matrix to sammon. The matrix is a co-occurance matrix with no missing data. If this is at all relevant, running hclust on this matrix works. > samx <- sammon(q23axproduct) I receive the following error: Error in sammon(q23axproduct) : initial configuration must be complete In addition: Warning messages: 1: some of the first 2 eigenvalues are

Hi! I come from Poland, and my English is not as good as it should be, so I have to apologize for my mistakes. I would like to ask you for a package, where I can find a test, which is named The Test Of Series (in Poland). In econometric models we use this test to check if the residuals are completely random ( I mean with destiny). And one more thing - I would like to know where I can find a test

Hi: Recently, I want to perform a transformation on my data to make it more normal, meanwhile the order statistics is unchanged. So I decided to use a box-cox transformation. below is the qq-plot of the original data http://n4.nabble.com/file/n1011015/start%2Bvalue%2Bproblem%2B02.jpeg Note that the min of my data is -1099, so I add a fix value 1200 to the original sample. I choose the

Hi, I would need some help to fix polybench/symm: void kernel_symm(int ni, int nj, DATA_TYPE alpha, DATA_TYPE beta, DATA_TYPE POLYBENCH_2D(C,NI,NJ,ni,nj), DATA_TYPE POLYBENCH_2D(A,NJ,NJ,nj,nj), DATA_TYPE POLYBENCH_2D(B,NI,NJ,ni,nj)) { int i, j, k; DATA_TYPE acc; /* C := alpha*A*B + beta*C, A is symetric */ for (i = 0; i < _PB_NI; i++) for (j = 0; j < _PB_NJ; j++) {

Hi there! I installed the BudgetTone (GrandStream) on my LAN without any problems. Then, I moved it to another location using a D-Link NAT. I opened 5060 (SIP) and 5000 to 5008 for RTP. I also fixed the IP address of the BudgetTone. When I receive a call on my Asterisk, it would ring my FXS as before. However, after I pick up, it hangs within a few seconds (Hungup Zap1-1 in the log). The

I frequently want to test for differences between animal size frequency distributions. The obvious test (I think) to use is the Kolmogorov-Smirnov two sample test (provided in R as the function ks.test in package ctest). The KS test is for continuous variables and this obviously includes length, weight etc. However, limitations in measuring (e.g length to the nearest cm/mm, weight to the nearest

Hello, I am running R under Linux/x11. I would like to call R from a shell script and have it display a series of graphics. The graphics should remain visible until the user clicks or presses a key. I first tried R BATCH, but it does not load the x11 module, making it impossible to open x11 or png devices. Then, I tried to call R with a 'here' script: R --vanilla --quiet --args

Hello, Gurus: I tried to test if the sample mean of a dataset is zero. The data has 1500 numbers with a lot of zeros and some small positive numbers. The data range on [0,1] but the distribution is unknown. It is zero inflated anyway. I tried to use the Wilcoxon Signed Ranks test. But I read from this website that it does assume the population pdf is symmetric.

Hi, I am trying to write a function with the following codes and I would like it to return the values for "alpha beta para parab " seperately. Then I would like to use this funstion for "variable" with factor "a" and "b". But the result turns out to be a matrix with element like "Numeric,2" ... I guess they are just the values for

Dear All, I am doing some exercise in statistics textbook on comparison of two experimental means. Is it possible to use t.test() do t-test when I have only two means, sample size, two standard deviations ? (no raw data). Thanks. Pramote

Is there a conventional way to test for nested factors? I.e., if 'a' and 'b' are lists of same-length factors, does each level specified by 'a' correspond to exactly one level specified by 'b'? The function below seems to suffice, but I'd be happy to know of a more succinct solution, if it already exists. Thanks, Tim. --- "%nested.in%" <-

I'm trying to learn how to do a repeated measures ANOVA in R using lme(). A data set that comes from the book Design and Analysis has the following structure: Measurements (DV) were taken on 8 subjects (SUB) with two experimental levels (GROUP) at four times (TRIAL). In SAS, I use the code: PROC MIXED DATA=[data set below]; CLASS sub group trial; MODEL dv = group trial group*trial;

I have 2 files containing data analysed by 2 different methods. I would like to find out which genes appear in both analyses. Can someone show me how to do this? _________________________________________________________________ [[trailing spam removed]] [[alternative HTML version deleted]]

I have a bunch of benchmark measurements that look something like this: sample.1 0.0000066660 0.0000062500 0.0000058330 0.0000058330 0.0000058330 sample.2 0.0000058330 0.0000058330 0.0000058330 0.0000058330 0.0000058330 sample.3 0.0000062500 0.0000062500 0.0000070830 0.0000062500 0.0000066660 i.e each measurement take on one of a set of values. The set values isn't fixed, but

Hi, all I have a rather un-ideal dataset that I am trying to work with, and would appreciate any advice you have on the matter. I have 4 years worth of data taken at 3 depth-zones from which samples have been taken at random. I am looking at the abundance of organism A between depth zones and across years, and am interested in the possible interaction of organism A distributions shifting

Hello, I routinely use aov and and the Error term to perform analyses of variance of experiments with 'within-subject' factors. I wonder whether a notion like 'multistratum models' exists for glm models when performing a logit analysis (without being 100% sure whether this would make sense). I have data of an experiment where the outcome is a categorical variable: 20