similar to: AR(1) models with gls

Hi there, I am using gls from the nlme library to fit an AR(1) regression model. I am wondering if (and how) I can separate the auto-correlated and random components of the residuals? Id like to be able to plot the fitted values + the autocorrelated error (i.e. phi * resid(t-1)), to compare with the observed values. I am also wondering how I might go about calculating confidence (or

Hi, I was wondering if someone could please help me. I am doing a logistic regression to compare size at maturity between 3 seasons. My model is: fit <- glm(Mature ~ Season * Size - 1, family = binomial, data=dat) where Mature is a binary response, 0 for immature, 1 for mature. There are 3 Seasons. The Season * Size interaction is significant. I would like to compare the size at 50%

Dear r-help, Sleep-deprivation from having 2 youngsters under 2 around the house is fuzzing my brain, so please be gentle if the answer to this query is obvious! In the example below, I'm trying to use grid.rect to add grey rectangles to the panels of a lattice plot to indicate which months spawning occurred of a (very cute) native Tasmanian fish. The fish in the two lakes spawned at

Does anyone know how to do hermite or Bezier splines in R? I can find Matlab routines, but really need to implement them in R. Failing that I'd be interested in other conformal splines. I need to smooth and interpolate animal tracking data. Antarctic Wildlife Research Unit School of Zoology University of Tasmania PO Box 252-05 Hobart TAS, 7001 Australia Phone: (0)3 6226 2645 Fax:

On Thu, 4 Aug 2016 at 17:53 Michael Sumner <mdsumner at gmail.com> wrote: > On Thu, 4 Aug 2016 at 11:17 Paul Murrell <paul at stat.auckland.ac.nz> wrote: > >> Hi >> >> Just to clarify, I think this IS a problem with grid.path() as well as >> polypath(). >> >> > Hi, oh dear - sorry about that > > I appreciate the deeper explanation, I

On Thu, 4 Aug 2016 at 11:17 Paul Murrell <paul at stat.auckland.ac.nz> wrote: > Hi > > Just to clarify, I think this IS a problem with grid.path() as well as > polypath(). > > Hi, oh dear - sorry about that I appreciate the deeper explanation, I knew about the id aspect in grid, but just forgot in my haste. I'll be more careful with examples if I find any more clues.

Hello, it's probably worth adding that this is not a problem with pathGrob, only polypath. This code is sufficient to demonstrate the problem in Windows. ## overlapping, both clock-wise x <- cbind(c(.1, .1, .6, .6, NA, .4, .4, .9, .9), c(.1, .6, .6, .1, NA, .4, .9, .9, .4)) ## only a problem on Windows windows() and png() plot(x);polypath(x, rule = "winding", col =

This is now fixed in development on RForge, you can try it out by installing from there, or from the Github mirror with devtools::install_github("rforge/raster/pkg/raster"). (To get fixes into raster email the maintainer directly - you might not get a response but it'll be addressed). Cheers, Mike. On Thu, 24 May 2018 at 20:08 Michael Sumner <mdsumner at gmail.com> wrote:

Hi Just to clarify, I think this IS a problem with grid.path() as well as polypath(). For the example you give, grid.path() diverts to drawing a polygon (because there is no 'id' specified), and the NAs in 'x' generate two separate polygons, which get drawn one on top of the other. The correct analogy to the polypath() example is ... x2 <- matrix(x[!is.na(x)], ncol=2)

Hi, Is it possible to have non-authenticating shares on an server with security=ADS ? I have a RHEL server, with Centrify Express, and joined to a domain, but I would like to have a samba share that doesn't request a username/password for machines not on the domain. When I have a plain windows XP machine (not on the domain) attempt to connect, I get asked for a username/password. Is

Try text(0.5, 0.5, label = "?text") On Sat, 10 Feb 2018, 16:22 Mohammad Tanvir Ahamed via R-help, < r-help at r-project.org> wrote: > Hi, > I want to label a polygon (circle or polygon) inside. > As for example code > > library(spatstat) > x <- runif(20) > y <- runif(20) > X <- ppp(x, y, window=disc(0.7)) > plot(X) > > Now I want to

Dear All, I am fitting scallop count data to negative binomial GAMs. I have two significant parameters that explain 43%of the deviance. The adjusted r square is 0.25. The gam.check function gives me the figure attached. In the graph of linear predictor vs. residuals there seems to be more negative residual values than positive. Is that telling me that the fit is underestimating the response? Can

Hi, I see different results in png() and pdf() for polypath() on Windows when using the "winding" rule ## overlapping, both clock-wise x <- cbind(c(.1, .1, .6, .6, NA, .4, .4, .9, .9), c(.1, .6, .6, .1, NA, .4, .9, .9, .4)) pfun <- function() { plot(x) polypath(x * 0.8 + 0.2, rule = "winding", col = "#BEBEBE80") polypath(x, rule =

Hi, Is it possible to have non-authenticating shares on an server with security=ADS ? I have a RHEL server, with Centrify Express, and joined to a domain, but I would like to have a samba share that doesn't request a username/password for machines not on the domain. When I have a plain windows XP machine (not on the domain) attempt to connect, I get asked for a username/password. Is

On 20 April 2017 at 09:33, Michael Sumner wrote: | Perhaps worth pointing out "yesno" as a candidate: | | https://CRAN.R-project.org/package=yesno | <https://cran.r-project.org/package=yesno> Not really as i) base R will not have dependencies outside of it ii) more importantly, Duncan asked for a good example of a "consistent" (as in: predictable) question

Hi All I'm fairly new at R so there might be on obvious solution to my problem with calling the lattice function xyplot( ) from my own function. The function script is given below. The problem is that the xyplot() call does not do anything! I just get a blank graphics window and no error messages. The simple plot() function call works fine when not commented out so I know the data is

Hi all, Using R v2.7.1, platform i386-pc-mingw32 Can someone please shed some light on the behaviour of ifelse() for me? My intent is to calc relative proportions of z$b, at the same time subsetting z$b based on z$a. I could attack the problem other ways (suggestions welcome) but I am also intrigued by the _order_ in which ifelse seems to assign values, and how recycling works. For instance,

Hi, I have fitted a model using a glm() approach and using a gls() approach (but without correcting for spatially autocorrelated errors). I have noticed that although these models are the same (as they should be), the AIC value differs between glm() and gls(). Can anyone tell me why they differ? Thanks, Geertje ~~~~ Geertje van der Heijden PhD student Tropical Ecology School of Geography

Perhaps worth pointing out "yesno" as a candidate: https://CRAN.R-project.org/package=yesno <https://cran.r-project.org/package=yesno> Cheers, Mike On Wed, 19 Apr 2017 at 21:55 Duncan Murdoch <murdoch.duncan at gmail.com> wrote: As described in <https://bugs.r-project.org/bugzilla3/show_bug.cgi?id=17242>, R base functions are currently inconsistent in asking

Hi I am struggling with the correct code syntax to generate xyplots that are different from the default graphics option. This is simplified version of my data: site size dnaconc prey A 44 22.6 1 A 47 18.2 1 B 38 25.6 0 B 47 20.5 1 C 45 30.5 0 C 32 18.5 0 D 45 15.0 1 D 22 20.1 0 So I am using xyplot(dnaconc~size|site,data=bug,groups=prey) to generate 4 panels (one for each site). My