similar to: Beginner's Question on Linear Regression Models

Hi, does anyone out there have a recommendation for multilevel / random effects and longitudinal analysis? My dream book would be something that's both accessible to a non-statistician but rigorous (because I seem to be slowly turning into a statistician) and ideally would use R. Peter

[Using R 2.7.2 on Windows XP] After re-building our heplots package, I've begun to get the following error from sessionInfo(), even though it passes R CMD check and builds without errors: > sessionInfo() Error in x$Priority : $ operator is invalid for atomic vectors In addition: Warning message: In FUN(c("MASS", "heplots", "car", "rgl",

Hello, I'm writng this message to see if anyone knows how to analyze a data with geometrically distributed outcome. My situation now is : I have a outcome which is geometrically distributed. And I've also self-derived the link function for geometric distribution. I want to use generalized linear models to analyze it, but there're only models for binary data, poisson distributed

How about this (I'm showing it as a pipe because it's easier to read that way): library(magrittr) "f 147/1315/587 2820/1320/587 3624/1321/587 1852/1322/587" %>% ? strsplit(' ') %>% ? unlist %>% ? sub('^[^/]*/*','',.) %>% ? sub('^[^/]*/*','',.) %>% ? paste(collapse = ' ') Georges Monette -- Georges Monette,

I have a file containing "words" like a a/b a/b/c where there may be multiple words on a line (separated by spaces).? The a, b, and c strings can contain non-space, non-slash characters. I'd like to use gsub() to extract the c strings (which should be empty if there are none). A real example is "f 147/1315/587 2820/1320/587 3624/1321/587 1852/1322/587" which

Hi, I couldn't resist these two suggestions: strings <- c("ASk/20005-01-45/90", "Alldatk/25-17-4567/990") x <- as.numeric(gsub("^[^-]*-|-.*$","",strings)) or x <- as.numeric(sub("^[^-]*-([0-9]+)-.*$","\\1",strings)) Best, Georges --------------------- Georges Monette, York University, Toronto On 08/04/2016 10:53 PM,

Dear R list members, I'd like to produce rug plots at the bottom of panels in a trellis display (using the lattice package), but par("usr") doesn't return user coordinates for panels, and consequently rug fails, as the following example (suggested to me by Georges Monette) illustrates: > x <- rnorm(50) > y <- rnorm(50) > f <-

Hi, I'm using glm function to do logistic regression and now I want to know if it exists a kind of R-squared with this function in order to check the model. Thank you.

Thomas, may I also suggest, from the Documentation>Contributed section of CRAN, "Econometrics in R" by Grant Farnsworth http://cran.at.r-project.org/doc/contrib/Farnsworth-EconometricsInR.pdf (see the chapter on Time series) and, in case you can read Italian, "Analisi delle serie storiche con R" by Vito Ricci http://cran.at.r-project.org/doc/contrib/Ricci-ts-italian.pdf

We are pleased to announce the following book: Title: Beginner's Guide to Zero-Inflated Models with R Authors: Zuur, Ieno Book website: http://www.highstat.com/BGZIM.htm Paperback or EBook can be order (exclusively) from: http://www.highstat.com/bookorder.htm TOC: http://www.highstat.com/BGS/ZIM/pdfs/TOCOnly.pdf Keywords: 430 pages. Zero inflated count data. Zero inflated continuous data.

Sorry if this is obvious, but my S skills aren't great and I haven't been able to find it documented anywhere. I want to write a new function for use with lme objects; the function will simply calculate an ICC (aka "rho") for each level of a mixed-effects model. What I need for this is pretty simple: (c(var1..varn, residual)) / sum(c(var1..varn, residual)) where var1..varn

After upgrading to 1.7.0 under debian linux, I can't get e1071 working properly. The first problem I had was that g++-3.0 was the standard compiler but wasn't installed, so I installed it. e1071 then installed correctly, but I get the following: aperrin at perrin:~/afshome/papers/authoritarian/R$ R R : Copyright 2003, The R Development Core Team Version 1.7.0 (2003-04-16) R is free

Greetings- An odd situation has developed. I use the following code to create .eps files of two very similar graphs: postscript(file='resources.bygt.eps', onefile=FALSE, horizontal=TRUE) barplot(resources.bygt.matrix, beside = TRUE, legend.text=c('narrative','doubt'),

We would like to announce the following statistics course: Course: Introduction to Linear Mixed Effects Models and GLMM with R Where:? Lisbon, Portugal When:?? 19-23 February 2018 Course website: http://highstat.com/index.php/courses Course flyer: http://highstat.com/Courses/Flyers/2018/Flyer2018_02LisbonV2.pdf Kind regards, Alain Zuur -- Dr. Alain F. Zuur Highland Statistics Ltd. 9 St

Dear R users, I have a couple of linear model related questions. 1) How do I produce a fixed effect linear model using lme? I saw somewhere (this may be Splus documentation since I use Splus and R interchangeably) that using lme(...,random= ~ -1 | groups,...) works, but it gives the same as lme(...,random= ~ 1 | groups,...), ie. fits a random effect intercept term. The reason why I want to do

Please forgive my ignorance on grapics. I'm trying to make a relatively simple plot with two line plots, same axes, by mean over a series of dates. I can make the plot well like this: plot(sort(tapply(first.anti.auth.sum,date,mean), partial=1), type="l", col="yellow",ylim=c(0,2.0)) par(new=TRUE) plot(sort(tapply(first.pro.auth.sum,date,mean), partial=1),

Dear Val, Your question isn't entirely clear (to me), but this is what I think you want to do: ------------------ snip ---------------- > strings <- c("ASk/20005-01-45/90", "Alldatk/25-17-4567/990") > location <- regexpr("-[0-9]*", strings) > x [1] "01" "17" > x <- substring(strings, location + 1, location +

Hi can anyone recommend a perl module that I can use to run R? Stephen

Hi, Can R calculate an intraclass correlation coefficient for clustered data, when the outcome variable is dichotomous? By now I calculate it by hand, estimating between- and intracluster variance by one-way ANOVA - however I don't feel very comfortable about this, since the distributional assumptions are not really met.... Maybe anyone can help me? Best regards and many many thanks,

Returning after the long weekend, I get the following: aperrin at perrin:~/afshome/papers/microcultures/R$ R R : Copyright 2003, The R Development Core Team Version 1.7.0 (2003-04-16) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type `license()' or `licence()' for distribution details. R is a collaborative