similar to: Time-varying correlation calculation

Dear experts, I am asking for help with a question regarding to stratified bootstrap. My dataset is a longitudinal dataset (3 measurements per person at year 1, 4 and 7) composed of multiple clinic centers and multiple participants within each clinic. It has missing values. I want to do a bootstrap to find the standard errors and confidence intervals for my variance components. My model is a

There is now a tutorial on bootstrapping and other resampling methods at: http://www.burns-stat.com/pages/Tutor/bootstrap_resampling.html Corrections and other suggestions are welcome. The project started because a novice asked me about bootstrapping. My response was, "How dare you bug me while I'm playing with my cats, just google for it." My correspondent was not very impressed

Hello all, I am trying to run a two sample t-test on a matrix which is a 196002*22 matrix. I want to run the t-test, row-wise, with the first 11 columns being a part of the first group and columns 12-22 being a part of the second group. I tried running something like (temp.matrix being my 196002*22 matrix) t.test(temp.matrix[,1:11],temp.matrix[,12:22],paired=TRUE) or somthing like

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

Hello, I learnt there's an advanced R/Splus course in Boston this january. Anyone got the announcement? please kindly forward it to me. Best, Eugene

Sometimes one might like to obtain pointwise bootstrap bias-corrected, accelerated (BCA) confidence intervals for a large number of statistics computed from a single dataset. For instance, one might like to get (so as to plot graphically) bootstrap confidence bands for the fitted values in a regression model. (Example: Chiu S et al., Early Acceleration of Head Circumference in Children with

Dear R experts, I have defined a class "myclass" and would like the slots to be extractable not only by "@" but also by "$". I now try to write a method for "$" that simply executes the request object at slotname, whenever someone calls object$slotname for any object of class "myclass". I don't manage to find out how I can provide this

Hi there, I want to generate different samples using the followindg code: g<-sample(LETTERS[1:2], 24, replace=T) How can I specify that I need 12 "A"s and 12 "B"s? Thank you, Judith ____________________________________________________________________________________ Be a better friend, newshound, and

I apologize for this post. I am new to R (two days) and I have tried and tried to calculated confidence intervals for medians. Can someone help me? Here is my data: institution1 0.21 0.16 0.32 0.69 1.15 0.9 0.87 0.87 0.73 The first four observations compose group 1 and observations 5 through 9 compose group 2. I would like to create a bootstrapped 90% confidence interval on the difference of

If my sample size is small is there a particular switch option that I need to use with t.test so that it calculates the t ratio correctly? Here is a dummy example? á =0.05 Mean pain reduction for A =27; B =31 and SD are SDA=9 SDB=12 drgA.p<-rnorm(5,27,9); drgB.p<-rnorm(5,31,12) t.test(drgA.p,drgB.p) # what do I need to give as additional parameter here? I can do it manually but

----- Quick summary: In the File:Open dialog, please change "S files (*.q)" to "S files (*.q, *.ssc, *.S)" and show the corresponding files (including .SSC and .s files). ----- Background This is motivated by the following query to R-help: >Date: Thu, 16 Mar 2006 22:44:11 -0600 >From: "xpRt.wannabe" <xprt.wannabe at gmail.com> >Subject: [R] Is

I'm trying: R CMD Rdconv file.Rd --type=Ssgm If file.Rd contains \code{x} then this is currently translated as <s-expression>x</s-expression> I suggest instead translating to <code>x</code> (provided that R CMD Sd2Rd is changed to support the <code> tag; I just submitted that bug separately). Note that this is an enhancement, not a bug. This change would

Hi dear all, Can someone help me about detection of outliers using jackknife after bootstrap algorithm? -- View this message in context: http://r.789695.n4.nabble.com/jackknife-after-bootstrap-tp3041634p3041634.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]]

is() does not catch parent S3 classes: > library(splines) > temp <- bs(1:99, df=5) > class(temp) [1] "bs" "basis" > is(temp, "basis") [1] FALSE In contrast, is() does catch parent S4 classes: > library(copula) > norm.cop <- ellipCopula("normal", param = c(0.5, 0.6, 0.7), + dim = 3, dispstr = "un")

On Mon, 30 Apr 2007 bill at insightful.com wrote: > On Tue, 10 Apr 2007 timh at insightful.com wrote: > > > I've created a .Rd file (below), then converted that to .sgml using > > R CMD Rdconv --type=Ssgm combn.Rd > combn.sgml > > The output (shown below) is missing some of the sections: > > arguments > > seealso > > examples > > If

Insightful is seeking a pre-doctoral student and an undergraduate student for two internship positions. The primary responsibilities are to assist in the development of software for high-dimensional regression and machine learning applications using least angle regression (LARS). The pre-doctoral candidate should have a background and interest in statistical methodology, algorithms, data

On Tue, 10 Apr 2007 timh at insightful.com wrote: > I've created a .Rd file (below), then converted that to .sgml using > R CMD Rdconv --type=Ssgm combn.Rd > combn.sgml > The output (shown below) is missing some of the sections: > arguments > seealso > examples > If instead I convert to .d (below), the same sections are missing, > and the "note"

Doing help("R_LIBS") brings up a help file (the same one as help(library)), but the help file doesn't mention R_LIBS. It does have a link to .libPaths, which does document R_LIBS. The quickest fix would be for help("R_LIBS") to bring up the .libPaths help file. --please do not edit the information below-- Version: platform = i386-pc-mingw32 arch = i386 os = mingw32

Hi! When I apply the lars (least-angle-regression) method to my data (3655 features, only 355 data points, no I did not mistype), I observe a strange behaviour: 1) The beta values tend to grow into real high values quite fast up to a point where they overflow and get negative. The overflow is not a problem, I don't need the last part of the analysis anyway, but why do they just

library(glmpath) data(heart.data) # heart.data is a list, $y a vector, $x a matrix data <- data.frame(x=I(heart.data$x), y = heart.data$y) > data[1:2,] x.1 x.2 x.3 x.4 x.5 x.6 x.7 x.8 x.9 y 1 160 12 5.73 23.11 1 49 25.3 97.2 52 1 2 144 0.01 4.41 28.61 0 55 28.87 2.06 63 1 > dimnames(heart.data$x)[[2]] [1] "sbp"