search for: grueber

...f place here. Best wishes, Giovanni Giovanni Millo, PhD Research Dept., Assicurazioni Generali SpA Via Machiavelli 4, 34132 Trieste (Italy) tel. +39 040 671184 fax +39 040 671160 ----------------- original message ------------ Message: 8 Date: Wed, 2 May 2012 05:45:47 -0700 (PDT) From: Philipp Grueber <philipp.grueber at ebs.edu> To: r-help at r-project.org Subject: Re: [R] error in La.svd Lapack routine 'dgesdd' Message-ID: <1335962747113-4603097.post at n4.nabble.com> Content-Type: text/plain; charset=UTF-8 Dear R Users, I have an unbalanced panel with (on average) approx...

Dear all, I am working on a large dataset and need to use biglm() to perform OLS regressions. I have detected significant ARCH effects which I try to account for using the Newey-West correction. So far, I have worked with NeweyWest() in the sandwich package. NeweyWest() however seems to be unable to handle an object of class "biglm". Looking into the code, I figured out that

...:Propfarm+ (1|Socialgroup)+(1|Year)+(1|Tattoo)),REML=FALSE, data=propspec) pss<-standardize(ps,standardize.y = FALSE) psdrg<-dredge(pss) summary(model.avg(get.models(psdrg,subset=delta<2))) #-------------------------------------------------------------------------------------# REf -Grueber, C.E., Nakagawa, S., Laws, R.J. & Jamieson, I.G. (2011) Multimodel inference in ecology and evolution: challenges and solutions. Journal of evolutionary biology, 24, 699-711. Any help would be much appreciated Regards Andrew Robertson PhD student Centre for Ecology and Conservation Universit...

Hello R users I'm using MuMIn but for some reason I'm not getting the adjusted confidence interval and uncoditional SE whe I use model.avg(). I took into consideration the steps provided by Grueber et al (2011) Multimodel inference in ecology and evolution: challenges and solutions in JEB. I created a global model to see if malaria prevalence (binomial distribution) is related to any life history traits of 14 different birds species, while controling for Family and genus in a GLMM: global.m...

...Lik(logLik=log_Lik,start=c(coef(res1),1,1),method="nr") _______________________________________ Am I on the right track? Is there an easier way to do this? Did I miss something important? Any help is appreciated, thanks a lot in advance! Best, Philipp ______________________ Philipp Grueber EBS Universitaet fuer Wirtschaft und Recht Wiesbaden, Germany ----- ____________________________________ EBS Universitaet fuer Wirtschaft und Recht FARE Department Wiesbaden/ Germany http://www.ebs.edu/index.php?id=finacc&L=0 -- View this message in context: http://r.789695.n4.nabble.com/M...

...n though being close). I guess I am missing something. Every comment pointing me to the right solution would be of great value to me. Also, comments on how to increase the efficiency of my function would help! Please find an example based on the Grunfeld data below. Thank you very much! Philipp Grueber ########################## library(MASS) getQ<-function(object,t.index,i.index){ t.ind<-object[,t.index] i.ind<-object[,i.index] nam<-unique(i.ind) tim<-unique(t.ind) n<-nrow(object) N<-length(nam) T<-length(tim) I<-matrix(0,n,n) I[row(I)==col(I)] <- 1 I_...