search for: 7ay

..., or are these functions not yet built? Chris Rogers ----------------------------------------------------------------------- L C G Rogers, Professor of Probability tel:+44 1225 826224 Department of Mathematical Sciences fax:+44 1225 826492 University of Bath, Bath BA2 7AY, GB Email: L.C.G.Rogers at bath.ac.uk web: www.bath.ac.uk/~maslcgr/home.html ----------------------------------------------------------------------- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-...

Hi all R users ! I'm using gamm function from Simon Wood's mgcv package, to fit a spatial regression Generalized Additive Mixed Model, as covariates I have the geographical longitude and latitude locations of indexed data. I include a random effect for each district (dist) so the code is fit <- gamm(y~s(lon,lat,bs="tp", m=2)+offset(log(exp.)), random=list(dist=~1),

...een added. See ?adaptive.smooth. * The interface for adding user defined smooths has been modified and simplified. See ?smooth.construct. A fuller list of changes is at http://cran.r-project.org/web/packages/mgcv/ChangeLog -- > Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7AY UK > +44 1225 386603 www.maths.bath.ac.uk/~sw283 ------------------------------------------------------- -- > Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7AY UK > +44 1225 386603 www.maths.bath.ac.uk/~sw283 _______________________________________________ R-packag...

...een added. See ?adaptive.smooth. * The interface for adding user defined smooths has been modified and simplified. See ?smooth.construct. A fuller list of changes is at http://cran.r-project.org/web/packages/mgcv/ChangeLog -- > Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7AY UK > +44 1225 386603 www.maths.bath.ac.uk/~sw283 ------------------------------------------------------- -- > Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7AY UK > +44 1225 386603 www.maths.bath.ac.uk/~sw283 _______________________________________________ R-packag...

...h(names(intvl.periods.myrs)) { my.new.vector<-sum(intvl.periods.myrs[i] } This does not work and probably for obvious reasons. Thanks Martin -- Martin Hughes MPhil/PhD Research in Biology Rm 1.07, 4south University of Bath Department of Biology and Biochemistry Claverton Bath BA2 7AY Tel: 01225 385 437 M.Hughes at bath.ac.uk http://www.bath.ac.uk/bio-sci/biodiversity-lab/hughes.html

...here are two types of NaNs in the ieee 754 standard, QNaN (Quiet NaN) and SNaN (Signalling NaN). Does 8.29488e-311 look familiar (I didn't think so)? Regards, Simon -------------------------------------------- Simon Pickering MEng Dept. of Mechanical Engineering University of Bath Bath, BA2 7AY, UK Tel: +44 (0)1225 383314 Fax: +44 (0)1225 386928 -------------------------------------------- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or...

...y<=0 if (sum(iy)) mustart[iy] <- min(y[!iy])*.5 } else if (family$link=="inverse") { iy <- y==0 if (sum(iy)) mustart[iy] <- min(abs(y[!iy]))*.5 } }) best, Simon >- Simon Wood, Mathematical Sciences, University of Bath, Bath BA2 7AY >- +44 (0)1225 386603 www.maths.bath.ac.uk/~sw283/

...omment after it.... # << FIXME? -- reconstructing all NA object ... so I guess I'm really asking if there is any chance of fixing this soon, or whether I should just code up a work-around for predict.gam? Simon -- > Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7AY UK > +44 1225 386603 www.maths.bath.ac.uk/~sw283

...ot;\t") PA<-as.matrix(pa[,-c(1)]) rownames(PA)<-(pa[,1]) OCHAR<-c() source("DISPARITY.R") library(calibrate) ### Thanks Martin -- Martin Hughes MPhil/PhD Research in Biology Rm 1.07, 4south University of Bath Department of Biology and Biochemistry Claverton Bath BA2 7AY Tel: 01225 385 437 M.Hughes at bath.ac.uk http://www.bath.ac.uk/bio-sci/biodiversity-lab/hughes.html

HELLO ALL, I AM GETTING AN ERROR MESSAGE WHEN TRYING TO RUN A GAMM MODEL LIKE THE ONE BELOW. I AM USING R VERSION 2.14.1 (2011-12-22) AND MGCV 1.7-12. M1 <-gamm(DepVar ~ Treatment + s(Year, by =Treatment), random=list(Block=~1), na.action=na.omit, data = mydata, correlation = corARMA(form =~ Year|Treatment, p = 1, q = 0)) THIS IS THE ERROR MESSAGE Error in `*tmp*`[[k]] : attempt to

Can anyone tell me the trick for obtaining the smoother matrix from smooth.spline when there are non-unique values for x. I have the following code but, of course, it only works when all values of x are unique. ## get the smoother matrix (x having unique values smooth.matrix = function(x, df){ n = length(x); A = matrix(0, n, n); for(i in 1:n){ y = rep(0, n); y[i]=1; yi =

...an angle Can anyone tell me of a way I can have both - ie text that is a specified distance from the margin AND tilted at and angle Thanks Martin -- Martin Hughes MPhil/PhD Research in Biology Rm 1.07, 4south University of Bath Department of Biology and Biochemistry Claverton Bath BA2 7AY Tel: 01225 385 437 M.Hughes at bath.ac.uk http://www.bath.ac.uk/bio-sci/biodiversity-lab/hughes.html

Dear useRs, This is the first time I post to this list and I would appreciate any help available. I've used the excellent mgcv package for a while now to investigate geographical patterns of language variation, and it has has always worked without any problems for me. The problem below occurs using R 2.14.0 (both 32 and 64 bit versions in Windows and the 64 bit version in Unix) and mgcv (both

hi Using a ts or tprs basis, I expected gcv to decrease when increasing the basis dimension, as I thought this would minimise gcv over a larger subspace. But gcv increased. Here's an example. thanks for any comments. greg #simulate some data set.seed(0) x1<-runif(500) x2<-rnorm(500) x3<-rpois(500,3) d<-runif(500) linp<--1+x1+0.5*x2+0.3*exp(-2*d)*sin(10*d)*x3

Hello, I'm running mixed models in GAMM4 with 2 (non-nested) random intercepts and I want to include a spline term for one of my exposure variables. However, when I include a spline term, I always get reported degrees of freedom of less than 1, even when I know that my spline is using more than 1 degree of freedom. For example, here is the code for my model: >

Dear All Hopefully someone out there can point out what I am missing! I have a (large, several hundred) dataset of gardens in which over two years the presence/absence of a particular bird species is noted each week. I have good reason to believe there is a difference between the two years in the weekly proportion of gardens and would like to assess this, before going on to look in more detail at

...oothness selection are now available. * A Tweedie family has been added. * `gam.method' has been replaced (see arguments `method' and `optimizer' for `gam') For other changes see the changeLog. Simon -- > Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7AY UK > +44 1225 386603 www.maths.bath.ac.uk/~sw283 _______________________________________________ R-packages mailing list R-packages at r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages

...one know a better way to do it? Thanks Chris Rogers ----------------------------------------------------------------------- L C G Rogers, Professor of Probability tel:+44 1225 826224 Department of Mathematical Sciences fax:+44 1225 826492 University of Bath, Bath BA2 7AY, GB Email: L.C.G.Rogers at bath.ac.uk web: www.bath.ac.uk/~maslcgr/home.html ----------------------------------------------------------------------- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-...

...roblem to me and how to test whether it is occurring. I have just compiled R for the Zaurus 5500, and want to see whether it works or requires tweaking. Regards, Simon -------------------------------------------- Simon Pickering MEng Dept. of Mechanical Engineering University of Bath Bath, BA2 7AY, UK Tel: +44 (0)1225 383314 Fax: +44 (0)1225 386928 -------------------------------------------- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or...

Hello all, I was wondering if anyone is aware of formal approaches and tools for comparing partial response curves produced in GAM? My interest is in determining if two partial response curves are "statistically" different. I recognize that point-wise standard error estimates can be produced using the GAM package but Im not certain how to translate this into a statistical test for the