similar to: Unequal variance ANOVA using gls function in nlme

Hi I want to extract the random slopes from a lmer (I am doing a random regression), but are the answers obtained from ranef or coef? My model is: mod1<-lmer(B~ A +(A|bird), family=quasibinomial) And I want to obtain a slope for each individual bird but am not sure which output I need and can't find the answer anywhere. Thanks Sam Dr Samantha Patrick EU INTERREG Post Doc Davy 618

Hello, I've assembled a Multivariate ctv (with a lot of help from Achim Zeleis, who has now posted the view on CRAN). I'd be grateful for comments regarding missing packages / functions. Opinions on the organisation of the view would also be appreciated, as well as having any errors pointed out. I've adopted a rather broad and vague definition of "multivariate", which may

I'am trying to develop some code if R, which would correspond to what I did in SAS. The data look like: Treatment Replicate group1 GSI Control A 1 0.81301 Control B 1 1.06061 Control C 1 1.26350 Control D 1 0.93284 Low A 2 0.79359 Low B

I found this paper on ANOVA on unequal error variance. Has this be incorporated to any R package? Is there any textbook that discuss the problem of ANOVA on unequal error variance in general? http://www.jstor.org/stable/2532947?cookieSet=1

Hello, I would be extremely grateful if anyone is able to provide any (rather obscure) advice on using R with Condor. I think I'm following Xianhong Xie's instructions (R News 5(2) 13-15) correctly, but my job just stays held in the queue (for days / months). I've checked condor_status to make sure there are plenty of machines available, but can't see any way to attack the

Hello, (apologies, I'm not entirely sure whether this question is about R or my limitations with PowerPoint). I've submitted a paper (which has been accepted) but the journal now require me to submit graphs that are "editable in PowerPoint". I would be grateful for suggestions as to how I should do this. The best route seems to be to copy-and-paste the figures from the

Hi all, right now i am simply comparing means. obviously this can be done by the simple t.test respectively the welch test, if var.equal is set to FALSE. just like this t.test( Y ~ group) t.test( Y ~ group, var.equal = FALSE) now that i need to compare weighted means i am using the lm function as an adjusted t-test: like lmtest <- ( Y ~ group ) anova(lmtest)

Dear Douglas and R-help, Does lme assume normal distribution AND equal variance among groups like anova() does? If it does, is there any method like unequal variance T-test (Welch T) in lme when each group has unequal variance in my data? Thanks, Shirley

Dear list, Inspecting residuals of my linear models, I detected spatial autocorrelation. In order to take this into account, I decided to use the GLS method with the correlation = corGaus ( ~ X + Y). Then, I can sort my GLS models based on their AIC. But ... how to know the proportion of the variance explained by the best one (it can be best of the worst models) ? R-squared value has not the

Dear R-help members, I have a question regarding how to use varComb function to specify a variance function for the "weights" in the gls. I need to fit a linear model with heteroscedasticity. The variance function is exp(c0+nu0*W +nu1*W^2) where W is a covariate. Initially I want to use varFunc to define my own variance function following the instruction in the Pinheiro and Bates

Hello everybody, I have a dataset where each row has number of subjects and that gives me natural weights for the variance function. Additionally I see that variance increases with Age, which is a regressor. So using gls I have weights=varFixed(~1/n) but don't know how to include the extra effect of the regressor. Fitted values show a quadratic curve vs. age, not sure if that helps.

Hi all, I have probably a basic question, but I can't seem to find the answer in the literature or in the R-archives. I would like to do a robust ANCOVA (using either rlm or lmRob of the MASS and robust packages) - my response variable deviates slightly from normal and I have some "outliers". The data consist of 2 factor variables and 3-5 covariates (fdepending on the model).

Dear R-users I am relatively new to R, i hope my many novice questions are welcome. I have problems accessing some objects (specifically the random effects, correlation structure and variance function) from an object of class gls and lme. I used the following models: yah <- gls (outcome~ -1 + as.factor(Trial):as.factor(endpoint)+

A couple of questions: How to be sure that gls allowes errors to be correlated and/or have unequal variances? (is this on auto or is there a switch?) How to calculate confidence limits for a linear regresssion? -------------- next part -------------- A non-text attachment was scrubbed... Name: dthompson.vcf Type: text/x-vcard Size: 303 bytes Desc: Card for David Thompson Url :

On Tue, 2018-12-04 at 08:19 +0000, J Martin Rushton via CentOS wrote: > The '/etc/named.conf.rpmnew' file supplied is a bare minimum to > "configure the ... server as a caching only nameserver (as a localhost > DNS resolver only)". As soon as you start adding any structure to it > things change, not just are added to. See >

How does one implement a likelihood-ratio test, to test whether the variances of the random effects differ between two groups of subjects? Suppose your data consist of repeated measures on subjects belonging to two groups, say boys and girls, and you are fitting a linear mixed-effects model for the response as a function of time. The within-subject errors (residuals) have the same variance in

Hello, I can't figure out how to use the upper.panel and lower.panel options in pairs to alter the label options for either panel independently of the other. I would like to be able to show the pairwise scatter plots for the data as they are (a vanilla pairs plot?) but separately to be able to label the points according to a factor level. It is easy enough to do this independently, but I

Recently after upgrading to 1.2 Beta3 from 1.1.13, I started to defer mail. After looking into it deeper, it looks like Dovecot/sieve now asks for a tmp folder in the imap (maildir) folder. Once I added the tmp folder as the logs indicants, the mail was processed normally. I have split my Inbox and imap store to different folder. Looking at what it's trying to do, it looks to me that

Hello, In trying to resolve a problem with a cron job, we can see when the job starts by looking in the /var/log/cron log file. However, I was asked if when the job ends could also be logged. (It seems to be something that crops up every so often over the years.) I found on the 'net this article https://serverfault.com/questions/248915/crontab-is-there-any-log-with-begin-and-end-time which

I am trying to analyse a data with gls/lm using the following set of models prcn.0.lm <- lm( log10(Y)~(cond-1)+(cond-1):t ,prcn) prcn.1.gls <- gls( log10(Y)~(cond-1)+(cond-1):t ,prcn,cor=corAR1()) prcn.0.gls <- gls( log10(Y)~(cond-1)+(cond-1):t ,prcn) prcn.1m.gls <- gls( log10(Y)~(cond-1)+(cond-1):t ,prcn,cor=corAR1(),method="ML") I get the following AICs for these models: