OOPPS...my apologies, I don't know why I jumped straight to
correlation...for the coefficients of linear regression (and again based
only on my limited knowledge) it would be like this:

x <- runif(15,1,10)            # make play data
# don't worry, it's linear now...
y <- x + rnorm(15,0,0.5)

plot(x,y)

# now begin the jacknife...
# loop through and drop one
jack.cor <- matrix(NA,nrow=length(x),ncol=2)
for(i in 1:length(x)) {
    jack.cor[i,] <- lm(y[-i]~x[-i])$coeff
}

# some different ways you could look at the SLOPE:
plot(density(jack.cor[,2]))    # pdf of test statistic
plot(x,jack.cor[,2])        # which x values are "outliers"
plot(y,jack.cor[,2])        # which y values are "outliers"

# to get the intercept, same thing, different column
# example:
plot(density(jack.cor[,1]))

SORRY for any confusion!
Andy


On Tue, Jun 1, 2010 at 6:41 AM, Carmen Burghelea
<carmen_burghelea@yahoo.es>wrote:

>
> Hi there,
>
> I would very much appreciate if somebody can help me with a syntax example
> of
> jackknife to obtain estimates of correlation coefficients in R.
>
> Many thanks in advance.
>
> Carmen
>
>
>
>        [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-sig-ecology mailing list
> R-sig-ecology@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
>
>

	[[alternative HTML version deleted]]

