similar to: ANOVA in quantreg - faulty test for 'nesting'?

I am helping a student with some logistic regression analyses and we are getting some strange inconsistencies regarding a warning about factor levels being dropped when running predict.glm(, newdata = ournewdata) on the logistic regression model object. We have checked multiple times that the factor levels have been defined similarly on both data sets (one used to estimate model and the newdata)

I know I must be missing something obvious, but checking help and googling a bit did not turn up a useable answer. When I've estimated a glm() model object (my example is with just identity link with gaussian family so I could have used lm() instead), one of the terms returned in the model object is listed as $effects. What are these quantities? I have not been able to relate them to the

A guess (treat accordingly): Different BLAS versions are in use on the two different machines/versions. In one, near singularities are handled, and in the other they are not, percolating up to warnings at the R level. You can check this by seeing whether the estimated fit is the same on the 2 machines. If so, ignore the above. -- Bert Bert Gunter "The trouble with having an open mind

List: I'm helping a colleague with some Poisson regression modeling. He uses SAS proc GENMOD and I'm using glm() in R. Note on the SAS and R output below that our estimates, standard errors, and deviances are identical but what we get for likelihoods differs considerably. I'm assuming that these must differ just by some constant but it would be nice to have some confirmation

Hello, I have used nlme to fit a model, the R syntax is like fmla0<-as.formula(paste("~",paste(colnames(ldata[,9:13]),collapse="+"),"-1")) > fmla1<-as.formula(paste("~",paste(colnames(ldata[,14:18]),collapse="+"),"-1")) >

Dear moderator, If possible I would like to send in the following question for R-help: I am analyzing a small data set using PGLS with phylogenetic uncertainty taken into account and thereby including 100 potential phylogenetic tree scenarios. I've managed to run models on all of the different trees and performed model averaging to get parameter estimates for the intercept and most of the

I have a recurring graphics issue that I've not been able to resolve with R. If I make a series of regression estimates and then plot the estimated function for the regression lines over a scatter plot of the data, e.g., using a sequence of plot( ) and lines ( ) similar to those below

Hi, I am currently trying to do a GLMM on a dataset with percent cover of seagrass (dep. var) and a suite of explanatory variables including algal (AC) and epiphyte cover (EC), rainfall, temperature and sunshine hours. M2=glmer(SG~AC+EC+TP+SS+RF+(1|Location/fSi/fTr), family=binomial,data=data,nAGQ=1) As the dependent variable is percent cover, I used a binomial error structure. I also have a

I think I must be missing something obvious, but I'm having trouble getting a data transformation to work on groupings of data within a data frame (csss3) as defined by 2 factors (population, locid). The data are sorted by year within locid within population and I want to lag another variable (dbc), i.e, shift them down by 1 row replacing the first row with NA, within groups defined by

Just noticed that I get a similar error about object 'kronecker' in "Matrix" package when trying to load "lme4". So this is a more pervasive problem. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: brian_cade@usgs.gov tel: 970 226-9326 ----- Forwarded by Brian S

Hi, A following function was kindly provided by GGally?s maintainer, Barret Schloerke. function(data, mapping, ...) { p <- ggplot(data = data, mapping = mapping) + geom_point(color = I("blue")) + geom_smooth(method = "lm", color = I("black"), ...) + theme_blank() + theme(panel.border=element_rect(fill=NA, linetype =

Hi Everyone, I am uncertain that I am writing the contrast statements correctly. Basically, I'm unsure when to use a -1 and a 1 when writing the contrasts. Specifically I am interested in comparing the slopes between different temperature regimes. Temperature is therefore a factor. Time and percent are numerical. Using the gmodels package I made the following model:

Having finally updated from R 1.91 to R 2.2.0 with my installation of a new computer, I discovered that something has changed drastically about the way code for contributed packages is stored when installed in a local version of R. In the 1.* versions it was easy for me to go in and modify some of the code for a contributed package by using a text editor to change the script files (these

I've looked for this information in all the R help sources I could find and found nothing. Is it possible to use some function key to stop the execution of some R command without ending the R session (Windows, R 1.91)? I've several times started functions that for various reasons are not executing properly and it would be nice to stop them without killing the R session. I've been

Your understanding isn't similar to mine. Mine says robust/resistant methods are for data with heavy tails, not heteroscedasticity. The common ways to approach heteroscedasticity are transformation and weighting. The first is easy and usually quite effective for dose-response data. The second is not much harder. Both can be done in R with nls(). Andy From: Quin Wills > > I am

Hi I am using the quantreg package for median regression for a large series of subsets of data. It works fabulously for all but one subset. When it reaches this subset, R takes the command and never responds. I end up having to kill R and restart it. It appears to be something with the particular data subset, but I can't pinpoint the problem. Here are some details Operating system:

I recently lost the partitions on my hard drive (second time in 6 months) so I had to have our IT folks image all my files over to a new drive. I completely reinstalled R (now 2.15.2) and all my libraries to my computer (Dell Latitude running Windows 7). A few of my previous workspaces (created with R 2.14.1) can't be restored, reporting an error similar to the one I get when I try to

I have a set of values and their corresponding weights. I can use the function weighted.mean to calculate the weighted mean, I would like to be able to similarly calculate the weighted median and quantiles? Is there a function in R that can do this? thanks, Spencer [[alternative HTML version deleted]]

Dear colleagues, has anyone an idea how to carry out a nonparametric manova for comparing K groups? has anyone a good package of non-parametric stats? Thank you for your help. -- José David Gómez Químico Farmacéutico Universidad Nacional [[alternative HTML version deleted]]

Hello, I'm looking for a R function proceeding MRPP (Multi-Response Permutation Procedures). Is it available? Thanking you in anticipation, Jeanne Vallet, PhD student [[alternative HTML version deleted]]