similar to: Case weights in nlme models

Dear list, I was trying to create a VarClass for nlme to work with Fay-Herriot (FH) models. The idea was to create a modification of VarComb that instead of multiplying the variance functions made their sum (I called it varSum). After some fails etc... I found that the I was not getting the expected results because I needed to make sigma fixed. Trying to find how to make sigma fixed I run into

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 :

Works perfectly. Thank you. -Hugh Rand -----Original Message----- From: Spencer Graves [mailto:spencer.graves at pdf.com] Sent: Sunday, January 15, 2006 6:41 PM To: Rand, Hugh Cc: 'r-help at lists.R-project.org' Subject: Re: [R] trouble with extraction/interpretation of variance structure para meters from a model built using gnls and varConstPower How about this: >

Hello - I've been using R for years, but have always been able to find what I need already available. Now I find that I would like to write a new varFunc class for the lme() or nlme() packages. There is some guidance for this in Problem 4 Chapter 5 of Bates and Pinheiro Mixed Effects Models in S and S-Plus. However, I find that I am unable to even get started and so have just purchased

Hello. I am trying my hand at modifying the varFunc class varExp, but I must be missing a step. All I want to do right now is make a working copy of varExp, call it varExp2, and then later change it. coef.varExp2, coef<-.varExp2, and Initialize.varExp2 all seem to work properly after I construct them. I can successfully use the commands: v2 <- varExp2(form = ~age|Sex,fixed =

some times ago i was playing with coding, shannon theoremes and other stuff, i have tried without success to compress audio wave, and i have notice that simply oversampling audio material enacnhe a lot compression ratio i only take awav file, oversampled it by 20 tiimes and then compressi it using pkzip or rar. i don0't remember if i also do a CONSTANT PREDICTION, iony know that pkzip and

Hi, I am running R2.8.1 under Linux, and I am having trouble using the variance functions in nlme My basic model was something like: model0 <- lme( log(growth) ~ light * species.group , data=data, random=~light|species ) # with 20 odd species divided in 2 groups Following the methods in Pinheiro&Bates I tried to put a variance function in the model: model1 <- update(model0,

It seems like every time I try to do something a little different, I cannot get output saved just right. This is on RedHat 7.2 with R 1.4.1. The png output looks fine, but the eps output has the problem that the bounding box on the legend cuts the legend in half. I put a copy of a bad one here: http://lark.cc.ukans.edu/~pauljohn/ResearchPapers/meanProtest-box.eps When I asked about these

some times ago i was playing with coding, shannon theoremes and other stuff, i have tried without success to compress audio wave, and i have notice that simply oversampling audio material enacnhe a lot compression ratio i only take awav file, oversampled it by 20 tiimes and then compressi it using pkzip or rar. i don0't remember if i also do a CONSTANT PREDICTION, iony know that pkzip

Dear listers, glmmPQL (package MASS) is given to work by repeated call to lme. In the classical outputs glmmPQL the Variance Structure is given as " fixed weights, Formula: ~invwt". The script shows that the function varFixed() is used, though the place where 'invwt' is defined remains unclear to me. I wonder if there is an easy way to specify another variance

Hello everybody. I am trying to upgrade from R-1.9.0 to R-1.9.1 on a RedHat linux 2.4.18 system. I get the following error after "tar -xvzf R-1.9.1.tgz ; cd ./R-1.9.1/ ; ./configure" and "make" : [make works for 10 minutes ... snip ...] varExp text html latex varFixed text html latex varFunc

Hi, I would like to use a weighted lm model to reduce heteroscendasticity. I am wondering if the only way to generate the weights in R is through the laborious process of trial and error by hand. Does anyone know if R has a function that would automatically generate the weights need for lm? Thanks, -- Tom [[alternative HTML version deleted]]

Dear all, First, I would like to thank you for your immense work. My question is about a frequent topic which I am not able to solve - even after hours of search in the mailing lisy. I would like to analyse random-effects (and fixed-effects)models of longitudinal / panel data with sampling weights. I have an unbalanced panel of different individuals in 5 years and income data as well as their

I have been using gnls with the weights argument (and varConstPower) to specify a variance structure for curve fits. In attempting to extract the parameters for the variance model I am seeing results I don't understand. When I simply display the model (or use "summary" on the model), I get what seem like reasonable values for both "power" and "const". When I

In species distribution modeling where one uses a large sample of background points to capture background variation in presence\pseudo-absence or use\available models (0\1 response) it is frequently recommended that one weight the data so the sum of the absence weights is equal to the sum of presence weights so that the model isn?t swamped by an overwhelming and arbitrary number of background

Hi all, I am trying to predict a target variable that takes values 0 or 1 using the rpart command. In my initial dataset I have few positive observations of the target variable; therefore I have oversampled the rare event by a multiple of 6 (i.e. from 762 to 4572). However, in my results, I end up with a number of positives in one of the terminal nodes that is not divisible by 6. As I have the

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 [To: flac-dev@lists.sourceforge.net] On Mon, 11 Nov 2002, Josh Coalson wrote: > > Very strange. I wounder if my two problems are related. The > > bitbuffer > > gets data moved around when read it called, so maybe they are > > related. > > It really seems like there is some memory craziness going on. > > Gets

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 [To: flac-dev@lists.sourceforge.net] I've made my very first release of MMIOFLAC at <http://math.berkeley.edu/~roconnor/MMIOMP3/>. It is an IOProc DLL that plugs into OS/2 multimedia system. The DLL allows playback of FLAC files. Seeking and encoding will be added. I just haven't gotten around to it yet. Thanks for the easy to use

I would like to be able to extract the significance levels (or standard errors) of the variance-covariance parameters in nlme objects (i.e. from the varClasses and/or corClasses). Is there a way to do it, or do I have to reverse the intervals function calculations ? Thanks. Rob _________________________________________________________________ [[alternative HTML version deleted]]

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 [To: flac-dev@lists.sourceforge.net] I'm working on making a plugin to OS/2 Multimedia subsystem for FLAC. I can't seem to find example FLAC files anywhere on FLAC's website. Where can I get some test files (perferably with WAV equiavlents). - -- Russell O'Connor <http://www.math.berkeley.edu/~roconnor/> ``[Law