search for: fairbroth

...lap? but they do. Can anyone please explain where I've gone wrong with the following code, OR why it makes sense for the bands to overlap, despite the statistically significant beta coefficient? There may be a good statistical reason for this, but I'm not aware of it. Many thanks, Malcolm Fairbrother n <- 120 set.seed(030512) x <- rbinom(n, 1, 0.5) dat <- within(data.frame(x), ybe <- rbinom(n, 1, plogis(-0.5 + x))) mod1 <- glm(ybe ~ x, dat, family=binomial) summary(mod1) # coefficient on x is statistically significant at the 0.05 level? almost at the 0.01 level pred <- pr...

Dear all, I have a data frame f, with four variables: f <- data.frame(A=c(0,0,1,1), B=c(0,1,0,1), C=c(1,1,0,1), D=c(3,1,2,3)) f A B C D 1 0 0 1 3 2 0 1 1 1 3 1 0 0 2 4 1 1 1 3 I want to create a new variable (f$E), such that each of its elements is drawn from either f$A, f$B, or f$C, according to the value (for each row) of f$D (values of which range from 1 to 3). In the first row, D is

I have spatial data on a sphere (the Earth) for which I would like to run an gls model assuming that the errors are autcorrelated, i.e. including a corSpatial correlation in the model specification. In this case the distance metric should be calculated on the sphere, therefore metric = "euclidean" in (for example) corSpher would be incorrect. I would be grateful for help on how to

Dear all, I cannot seem to get "boot" running. Can anybody help? I've tried downloading and compiling the binary from <http://cran.r-project.org/web/packages/boot/index.html>, as well as the below. (Please see the scary-sounding error message, and the sessionInfo.) Neither works. Any help would be much appreciated. - Malcolm > install.packages("boot",