similar to: reference on permutation test

Hi, I have optimized the shrinkage parameter (GCV)for ridge and got my r2 value is 70% . to check the sensitivity of the result, I did permutation test. I permuted the response vector and run for 1000 times and draw a distribution. But now, I get r2 values highest 98% and some of them more than 70 %. Is it expected from such type of test? *I was under impression that, r2 with real data set

I am perfectly aware that this question is not an R question, at least not yet, but I have not succeeded in finding what I want in other ways, so ... What I am looking for are two algorithms, preferabley in Pascal, but other languages may do. For (a) systematic (complete) permutations for grouped data with unequal group sizes, and (b) random permutations for the same kind of data. I know

Hello. I wrote the hook function for deliver. I want to add support of editheaders in the plug-in for dovecot. For this purpose I wrote the function rarules_get_stream. Remover of headrs works properly, but adding does not work. I took Timo Sirainen's advice from http://markmail.org/message/skb6arnll5gaopdr . Do I use a correct way of creation of a message? I give backtrace and a

Greetings, I have two files which contain responses to a series of multiple choice questions. One file contains responses before an "intervention" and the other contains the responses afterward. There were three possible responses to each question: D, F, T (for Don't Know, False, and True). I would like to try McNemar's test to determine if there was any significant

The other day I was looking into one of the classics in resampling, Eugene Edgington's "Randomization Tests". This type of test is simple to do in R with things like a simple correlation, the sample () function is perfect for the purpose. However, things are more complex if you have grouped data, like a one-way ANOVA. The reason is that you have to avoid the consideration of

Is there a way of repositioning the y label in scatterplot3d so that it is parallel with the y axis? Many thanks Richard

Dear R experts, I would like to explore if and to what extent bootstrapping and permutation statistics can help me for my research (functional brain imaging). I am looking for an introductory textbook, rather legible. I have statistical knowledge, but I am definitely no statistical or mathematical guru. Do you have suggestions for a useful textbook? Thanks a lot, Peter

Hi I was wondering if there is a permutation test available in R for linear models with continuous dependent covariates. I want to do a test like the one shown here. bmi<-rnorm(100,25) x<-c(rep(0,75),rep(1,25)) y<-rnorm(100)+bmi^(1/2)+rnorm(100,2)*x+bmi*x H0<-lm(y~1+x+bmi) H1<-lm(y~1+x+bmi+x*bmi) anova(H0,H1) summary(lm(y~1+x+bmi)) But I want to use permutation testing to

> De : Bill Venables <venables at acland.qld.cmis.csiro.au> > Objet : Re: [R] paste ? > Date?: mardi 25 avril 2000 08:45 > (...) > Secondly, I'm curious about the history of this kind of plot. > I've only heard it called a "Tukey mean difference" plot, (and > Trellis graphics has a function, tmd(), that does it, but no one > knows about it...).

In the past there's been a need to provide consistency between symbolic links or repository metadata during a sync. Currently, rsync renames files piecemeal. The attached patch (extremely ugly) attempts to resolve this by foregoing the rename step until the end. It adds a new option (if we didn't, ls might catch up). There are several issues to get over. The first big one in

Hey gang, I was wondering if all that documentation could or has been put into PDF format. I usually like to download documentation and read it while I''m sitting comfortably at home and I don''t want to tie up the phone line all night. Thanks, Nino p.s. If so, please feel free to attach the PDF formatted document to my e-mail ;-)

I am using the following code to produce a graphic: library(lattice) postscript("figs%03d.eps", width = 6.0, height = 6.0, ?????????? horizontal = FALSE, onefile = FALSE, paper = "special") xyplot(cases~yr|agrp*sex,data=data[tse==0 & expgrp==1,], groups=source, pch=".", type="l", ?????? main="Exposure group 1, time since entry 0") The

I am new to R and am trying to fit a survival curve with a weibull hazard function to a set of data giving the probability of survival to age x, given the year of birth, in the form: Probability of survival: Birth year 1980 1981 ... 2003 .2 0.90 0.89 ... 0.87 1 0.80 0.81 ... 0.79 age 2 0.75 0.74 ... 0.73 3 0.70 0.69 ... 0.68 5 0.50 0.49 ... 0.43 10 0.30 0.31 ... 0.26 I would like to be

I am using the following code: library(lattice) data<-read.csv("data.csv") attach(data) fig<-xyplot(S_t~month|event, key = list(text=list(lab=c("Time to first CV event - Data", "Survival post first CV event - Model", "Survival post first MIA/CA event - Data",

I saw the manual of this function but not sure what to do. I have a array contain 7 numbers and want to choose 4 to do permutation test. But using this function with parameters as onet.permutation(scores,4), it returns 0. Instead, with no parameter, onet.permutation(), it returns something but every time it's different. Please tell me how am I supposed to do it. Thanks a lot -- View this

I am looking for an algorithm (written in R (preferably) or C, but even pseudo-code in a text book maybe fine) to determine the sign of a permutation. What is that? Well, a permutation is either even or odd, the sign is +1 or -1, respectively, see, e.g., http://en.wikipedia.org/wiki/Signature_of_a_permutation which also says >> In practice, in order to determine whether a given

R users, Is it possible to do 'exact' permutation tests in R? If I run a perm.test or permutation.test.discrete on a data set that has, for example, 720 possible permutations, will that test permute ALL of those possible permutations, or will it just permute 720 random permutations (thus allowing some possibilities to repeat and some not to be permuted at all)? I'm a little confused

Hi, Is there any way to find the permutation function of the sorting and to apply the function (or its inverse) elsewhere? For example, the following permutation function from the sorting in the matrix form is c(1,2,3), c(2,1,3) > sort(c("bc","ac","dd")) [1] "ac" "bc" "dd" I try to find it in the permutations/permute

Hi, My query is regarding permutation test and reshuffling of genotype/phenotype data I have been using the haplo.stats package of R. for haplotype analysis and I would like to perform an analysis which I'm requesting your advice. I have a data set of individuals genotyped for 12 SNP and a dichotomous phenotype. At first, I have tested each of those SNP independently in order to bypass

Dear R-help-list I have two character vectors a <- c("A", "B", "C") b <- c("A", "C", "B") Then sapply(a, function(i) grep(i, b)) computes the permutation of the entries in 'b' needed to bring 'b' into the same order as 'a'. I have searched around, but haven't been able to find any existing function