similar to: random effects

Dear R users, I am trying to draw error bars in a bar plot, I use this code (tried many others which did not work): install.packages() library(gplots) y <-c(39.02, 46.42) se <- c(7.57,7.35) plot <- barplot(y, beside=TRUE, col=0, ylim=c(0,47), axis.lty=1, main="far-red", xlab="latitude", names.arg=c("56N", "68N")) superpose.eb(plot, y, se,

Hello I have a data set with outlier and it is not normally distributed. I would instead like to use a more robust distribution like t-distribution. My question is if the coefficients of the regression are different from zero, but assuming a t-distribution. Could someone hint me what package to use or.... Thanks in advance Rosario

Hello I am trying to draw a basic black and white map of two European countries. After searching some key words in google and reading many pages I arrived to the conclusion that persp() could be used to draw that map. I have prepared three small example files, which are supposed to be the files required for running that function. xvector is a vector with the longitudes yvector is a vector with

Hi all- I fit a zero-inflated Poisson model to model bycatch rates using an offset term for effort. I need to apply the fitted model to a datasets of varying levels of effort to predict the associated levels of bycatch. I am seeking assistance as to the correct way to code this. Thanks in advance! Laura [[alternative HTML version deleted]]

Hello I have a data set summarized like this: File name= Height Group Ind Age Trait 1 1 1 20 1 1 2 21 1 2 1 22 1 2 2 21 1 3 1 24 1 3 2 45 1 4 1 23 1 4 2 26 2 1 1 45 2 1 2 12 2 2 1 25 2 2 2 26 2 3 1 45 2 3 2 43 2 4 1 23 2 4 2 47 . . . I would like to plot Trait ~ Age but a different plot for each Group. I tried: > plot(Height$Trait ~ Height$Age | Group) But does not work. Any

Hello I am using box.cox() and I get this error message: Warning message: 'box.cox' is deprecated. Use 'bcPower' instead. See help("Deprecated") and help("car-deprecated"). I went to help but I did not understand the explanation, I am still wondering what is really happening. Thanks /R

Hello I have a problem on keeping the format when I export a matrix file with the write.table() function. When I import the data volcano from rgl package it looks like this in R: > data[1:5,] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 100 100 101 101 101 101 101 100 100 100 101 101 102 102 [2,] 101 101 102 102 102 102 102

I am trying to model survival data with a Weibull distribution using survreg. Units are clustered two apiece, sometimes receiving the same treatment and sometimes opposing treatment.

Hello I am trying to get out of an lm model the fstatistics, however after I run the model I write > names(Model) and the fstatistic does not appear only these. names(Model) [1] "coefficients" "residuals" "effects" "rank" "fitted.values" [6] "assign" "qr" "df.residual"

R users doing data analysis may be interested in the following paper: http://methodsblog.wordpress.com/2009/11/13/first-paper-now-online/?utm_source=feedburner&utm_medium=feed&utm_campaign=Feed%3A+wordpress%2Fmethodsblog+(methods.blog) All data and R code is available. Alain ----- -------------------------------------------------------------------- Dr. Alain F. Zuur First author of:

Hello I have a model with several hundred Y variables, and also several 1000 X variables. The model is linear lm(Y ~ X). My questions are: 1.- how to avoid writing all Xs variables? is list() the right function? 2.- about the multiple Ys with dependence among some of them, how to incorporate that information in the linear model? Thank you Rosario

Hello I have a annoying problem with R (which I am running in MAC). Every time I enter a function and I press enter to run it, I get many additional "<?s", it is really bothering me. I have tried to google to find out how to fix the problem but I did not succeed. Anyone with similar experience that could help me to change the set ups of the program or something. Thanks in advance

Hello I seem to find only two types of rotation for the factanal function in R, the Varimax and Promax, but is it possible to run a orthogonal and oblique rotations in R? Thanks in advance Rosario

Hello I am trying to compare two models using anova(), however I get a message error (see below). In the net I only found some information on certain library(car) for which one should use anova with A capital letter (Anova instead of anova), but I could not find car library as it says it does not exist. > Model <- lm(interceptG ~ SW + TSC + FSC + PF + SlopeG + K, data=AllTrait) >

Dear All, the lm() function has the possibility to create a subset of the possible explaining variables that you have. However, in the help there is no example how to use this subset option. I tried the following: model<-lm(dependent.data$MPFD~.,data=explaining.data,subset=c(1,0,0,0,0,0,0,0,1,1,0,0)) MPFD is the dependent variable stored in the data frame dependent.data, and all 12

Hi - I read through dse package manual a bit. I'm not quite certain how I can use it to estimate a time varying coefficient regression model? I might pick up an inappropriate package. Any suggestion would be greatly appreciated. Thank you. rh

Hello, I am using R to analyze a large multilevel data set, using lmer() to model my data, and using anova() to compare the fit of various models. When I run two models, the output of each model is generated correctly as far as I can tell (e.g. summary(f1) and summary(f2) for the multilevel model output look perfectly reasonable), and in this case (see below) predictor.1 explains vastly more

Let me start with an example: par(mfrow=c(2,3)) for (i in 1:5){ x = rnorm(100) y = .5*x + rnorm(100, 0, sqrt(1-.5^2)) plot(x,y) } Note that there's five plots and six spaces for those plots via mfrow, leaving one row empty. Is there a way to have the bottom two plots centered? I think it looks weird to have them left-justified. Thanks in advance for the help! -- Dustin Fife PhD

Hi, I am trying to add multiple variance structures such as the first example below: vf1 <- varComb(varIdent(form = ~1|Sex), varPower()) However my code below will not work can anybody please advise me? VFcomb<-varComb(varExp(form=~depcptwithextybf),varFixed(form=~FebNAO)) also if you have two variables with the same weights function would you write that as:

Hello everyone, I'm quite new to R and am trying to run a logistic model to look at how various measures of boldness in individual animals influences probability of capture, however I also want to include random effects and I'm not sure how to construct a model that incorporates both of these things. Data was collected from 6 different groups of 6 individuals with 10 replicates for