similar to: randomisation

Hello, While doing some enrichment tests using chisq.test() with simulated p-values, I noticed some strange behaviour. The computed p-value was extremely small, so I decided to dig a little deeper and debug chisq.test(). I noticed then that the simulated statistics returned by the following call tmp <- .Call(C_chisq_sim, sr, sc, B, E) were all the same, very small numbers. This, at first,

It is not about "really arge total number of observations", but: set.seed(4711);tabs <- r2dtable(1e6, c(2, 2), c(2, 2)); A11 <- vapply(tabs, function(x) x[1, 1], numeric(1));table(A11) A11 0 1 2 166483 666853 166664 There are three possible matrices, and these come out in proportions 1:4:1, the one with all cells filled with ones being most common. Cheers, Jari

Hello, In an attempt to randomise the basic blocks in a function, is it possible that I can randomise the entry block as well? And maybe insert some instructions in the pass to call entry block while running the program ? Is it feasible? What does entry block consist of ? Thanks Tanu --------------------------------- Do you Yahoo!? Yahoo! Search presents - Jib Jab's 'Second

Tanu Sharma wrote: > Hello, > > In an attempt to randomise the basic blocks in a function, is it > possible that I can randomise the entry block as well? And maybe insert > some instructions in the pass to call entry block while running the > program ? > > Is it feasible? > > What does entry block consist of ? The entry block, by definition, is the first basic

I would like to be able to randomise presence-absence (i.e. binary) matrices whilst keeping both the row and column totals constant. Is there a function in R that would allow me to do this? I'm working with vegetation presence-absence matrices based on field observations. The matrices are formatted to have sites as rows and species as columns. The presence of a species on a site is

I've tried and give up. I have a matrix of say 200 columns and 400 rows. For each odd ( or even i suppose if i wanted to )column, I want to know the number of rows in which the value is greater than zero. So, I did sapply(tempMatrix,2,function(x) sum( x > 0 )) this almost works but i don't know how to tell it to only do the odd columns. my guess is , like everything else in R, this

Hi, I want to ue a randomisation test, to compare two sample means, for a small set of data, I've looked at towt.permutation, perm.test and permtest, as well as permTS(), which is the best one to use for such a simple purpose? I know perm.test package isn't used anymore and coins recommended, but reading the manual I couldnt find the function most like it, ore more likely I

Hi to all, I am looking for a randomisation procedure for a single matrix, including a possibility to set the number of randomisations and the to set the number of row and columns . Knut

Hello everybody, I have a question. I would like to get a correlation between constitutive and induced plant defence which I messured on 30 plant species. So I have table with Species, Induced defence (ID), and constitutive defence (CD). Since Induced and constitutive defence are not independant (so called spurious correlation) I should do a randomisation test. I have a syntax of my

I'd like to carry out a Monte Carlo simulation test where given data is a contingency table. I think this is something to do with using rmultinonom(), but I'm not sure how to code this, to simulate contingency tables. Could anyone please help with how to use R to simulate contingency tables like this? -- View this message in context:

Hi I want to use order() to get the order of a vector. But I would need a different behavior when ties occur: similar to the parameter ties.method = "random" in the rank() function, I would need to randomise the ties. Is this possible? Example: x <- rep(1:10, 2) order(x) [1] 1 11 2 12 3 13 4 14 5 15 6 16 7 17 8 18 9 19 10 20 order(x) [1] 1 11 2 12 3 13 4 14 5 15

Dear R-people, I hope asking this is not too cheeky, but I do have a R Problem. I hope that some of you like to play around with R and can help me. Its like this. I have several plant species (A,B,C) and 10 replicates per species. 5 plants per species are damaged, 5 not. I let a caterpillar feed on each plant and measured the growth of the caterpillars on control plants (CR) and on damaged

chisq.test with simulate.p.value=TRUE uses the Patefield algorithm, this is not documented, and the original reference is not given, as it ought to be. The reference is: Patefield,W. M. (1981) An efficient method of generating r * c tables with given row and column totals (algorithm AS 159). Applied Statistics 30, 91-97. Kjetil Halvorsen

Dear R-users, We expect to create N homogenic groups of n features from an experimentation including N*n mesures. The aim of this is to prevent from group effects. How to do that with R functionalities. Does anyone know any methodes enabling this ? Best regards. Alexandre MENICACCI Bioinformatics - FOURNIER PHARMA 50, rue de Dijon - 21121 Daix - FRANCE a.menicacci at fr.fournierpharma.com t??l

Dear R-users, We expect to form N homogeneous groups of n features from an experimentation including N*n data. The aim is to prevent group effects. How to do that with R functionalitites ? Does anyone know any way enabling this ? Example : 100 patients are observed. 3 biochemical parameters are mesured for each one (Red and white globules ans glycemia). Patient RG RW

Hello I'm very much a beginner on meta-analysis, so apologies if this is a trivial posting. I've been sent a set data from separate experimental studies, Treatment and Control, but no measure of the variance of effect sizes, numbers of replicates etc. Instead, for each study, all I have is the mean value for the treatment and control (but not the SD). As far as I can tell, this forces

Hi I have FC6 system in the office running SVN-trunk-r63567 It is behind a NAT router which I have configured to do port forwarding etc. Asterisk connects and registers correctly to my SIP service (Sipgate.co.uk) and I can make and receive calls from any SIP phone on the office LAN. The problem comes when I try to use a SIP phone at home (also behind a NAT router). The phone registers correctly

I''ve spent some time converting my website from php to RoR, and so far, I''m loving every bit of Rails. However, I''ve stumbled across a problem. On my old page, I''ve had a random quote display at a specific location, and I''d like to replicate this with RoR. Basically, I thought I would be able to add all quotes into a specific table in the database,

Dear all, I have a question about the chisq.test command. As an option one can chose the computation of p-values by Monte-Carlo simulation (simulate.p.value=T). Is there any documentation available how this calculations are done and how this simulation based test behaves in small samples? Thanks Klaus Abberger University of Konstanz, Germany [[alternate HTML version deleted]]

On Mon, 14 Feb 2005, Reid Spencer wrote: > isn't necessarily tied to LLVM. Anyway, lets cross that bridge when we > get there. Sounds good. -Chris -- http://nondot.org/sabre/ http://llvm.cs.uiuc.edu/