similar to: What is the test statistics in perm.test

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 I am working on a Nested one-way ANOVA. I don't know how to implement R code to test the significance of the random factor My R code so far can only test the fixed factor : anova(lme(PCB~Area,random=~1|Sites, data = PCBdata)) numDF denDF F-value p-value (Intercept) 1 12 1841.7845 <.0001 Area 1 4 4.9846 0.0894 Here is my data and my hand

Hi! I am doing a two sample permutation test and trying to find confidence interval. perm.test() seems can do that but when the vector contains negative numbers it occurs an error. And for perm.test(x,y), elements in x can't be more than in y. I am wondering how I can fix these two problems? Thanks a lot. -- View this message in context:

Hi All: OK, this is probably a stupid thing to do but....... Say I have 706 controls and 6 treatments, and want to do a basic two-sample permutation test. I could do a monte carlo approximation but instead want to use perm.test. Now it happens that my data is counts with > unique(x) [1] 1 2 > unique(y) [1] 1 2 5 3 10 9 7 4 6 11 8 23 18 30 16 29 42 so there aren't that many

Hi, Please suggest a method to answer below questions: Factory_ID Factory_Location Factory_Size Total_Sample Good_Sample Fair_Sample Bad_Sample ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ 1 City_A

Hello list, I wish to perform a randomization test on the F-statistics of a 2 way ANOVA but have not been able to find out how to do so - is there a package / function that can perform this that I am unaware of? FactorA has 6 levels (0,1,2,3,4,5) whereas FactorB has 3 (1,2,3). A sample: Resp. FactorA FactorB 2 0 2 3 1 2 1 2 2 0 3 2 0 4 2 0 5 2 4 0 1 6 1 1 1 2 1 0 3 1 1 4 1 0 5 1 2 0 2 3 1 2 1

Hi all, Suppose we have the following matrix m <- matrix(c(1,2,3,2,1,3,3,1,2), ncol = 3, byrow=T) where in each row each number occurs only once. I'd like to define a permutation, e.g. 1 -> 2, 2 -> 1, 3 -> 3 and apply it to the matrix. Thus, the following matrix should result: m.perm <- matrix(c(2,1,3,1,2,3,3,2,1), ncol = 3, byrow=T) i.e. each 1 should map to 2 and vice

Hi everyone, I am somewhat new to R and I am trying to write a permutation function such that it inputs a character vector and from an arbitrary length "n" which is the length of the combinations for the character vector. I know there are R packages for permutation but this is for an assignment. So far this is what I have: alphabet =

Dear R Programmers, I am trying to run a Poisson regression on all pairs of variables in a data set and obtain the permutation distribution. The number of pairs is around 100000. It seems my code will take weeks to run, unless I try something else. Could you give me any suggestions on how to improve the speed of the code below, or any general suggestions on how I may accomplish this task. Thanks

Hi, Are you guys interested in implementing if-conversion as a GSoC 2015 project? Last year, I did a literature review about approaches of if-conversion and the if-conversion in LLVM. This was the undergraduate thesis of my bachelor degree. It seems that, the if-conversion used in LLVM is a very simple approach instead of following the literature. So I want to implement the approaches in the

Hi, I'm trying to run a permutation test on paired samples. First I tried the package "exactRankTests": require("exactRankTests") x <- c(1.83,0.50,1.62,2.48,1.68,1.88,1.55,3.06,1.30) y <- c(0.878,0.647,0.598,2.05,1.06,1.29,1.06,3.14,1.29) wilcox.test(x,y,paired = TRUE,alternative = "greater") perm.test(y,x,paired = TRUE,exact = TRUE,alternative =

OK, Let me describe. There is nothing wrong with if-conversion in LLVM. The algorithm implemented in LLVM can handle the if(???){do something} and if(???){do something}else{do something else} case very well. But it can handle complicated case like when there are a lot of gotos in the program. The more systematic way to do if-conversion is based on Hyperblock [Scott A. Mahlke et al 1992]

Hi! I am doing a permutation test on 7 data randomly selected into two groups, one has 4 and the other 3. The purpose is to list all the 35 possible outcomes and calculate p-value. I googled the function perm.test() and not sure if my usage is correct. So I set a=c(7 data), b=(1,1,1,1,2,2,2) which represents the group. Then I use the function pt=perm.test(a,b...). There's a p-value popped

Hi, I have a little problem that puzzles me about contradictory results returned by a bootstraped t-test (as in MASS 3rd ed p. 146) versus a permutation t-test (as in MASS 3rd ed, p 147). Data are measurements done on 100 cells in 9 slides of fish blood. With one method, cells are randomly sampled, and with the other method, the operator selects cells arbitrarily (in a way it is done usually

Dear all, I did a RDA and when I looked to the signification of the test with permutest, the output was non-significant. But when I used the envfit function, some of the vectors are significant. All the test's conditions are respected. What it means? Is it an error in the script? Commands and output: > permutest(rda.ind, perm=999, first=TRUE) Permutation test for rda  Call:

Hello, I am using logistic discriminant analysis to check whether a known classification Yobs can be predicted by few continuous variables X. What I do is to predict class probabilities with multinom() in nnet(), obtaining a predicted classification Ypred and then compute the percentage P(obs) of objects classified the same in Yobs and Ypred. My problem now is to figure out whether P(obs) is

Hi *, is there a way to obtain some kind of p-value for a model fitted with RMA using the lmodel2 package? I know that p-values are discussed and criticized a lot and as you can image from my question I'm not very much of a statistican (only writing my bachelor thesis). As fare as I understood the confidence interval statistic correctly, a coefficient is regarded as statistically significant

Hi, Urgent help- I have not been using R and statistics in my research for a long time, but still remember some concept. I would like to have a sample code for Manova analysis of Split-plot experiment. Could someone please post a sample code and a short input sample as well? Thank you so much! [[alternative HTML version deleted]]

Hello dear subscribed Users, this is my first post, so please forgive me for any inconveniences. The following problem: I have a dataframe containing a factor column. For each column i would like to compare means as parted by the factor. Using the normal t.test function I have already achieved my goal. But i think, that for my original data, a permutational t.test would be better suited.

Hi Al, On 2023/3/10 11:15, Al Viro wrote: > On Thu, Mar 09, 2023 at 11:21:23PM +0800, Yangtao Li wrote: >> Use i_blockmask() to simplify code. > > Umm... What's the branchpoint for that series? Not the mainline - > there we have i_blocksize() open-coded... Actually Yue Hu sent out a clean-up patch and I applied to -next for almost a week and will be upstreamed for