similar to: gamm

Hello, Why would I get an error message with the following code for gamm? I want to fit the a gam with different variances per stratum. library(mgcv) library(nlme) Y<-rnorm(100) X<-rnorm(100,sd=2) Z<-rep(c(T,F),each=50) test<-gamm(Y~s(X),weights=varIdent(form=~1|Z)) summary(test$lme) #ok summary(test$gam) Gives an error message: Error in inherits(x, "data.frame")

Dear R-helpers, I'd like to understand how to test the statistical significance of a random effect in gamm. I am using gamm because I want to test a model with an AR(1) error structure, and it is my understanding neither gam nor gamm4 will do the latter. The data set includes nine short interrupted time series (single case designs in education, sometimes called N-of-1 trials in medicine)

Dear R-help users, I am trying to analyze some visual transect data of organisms to generate a habitat distribution model. Once organisms are sighted, they are followed as point data is collected at a given time interval. Because of the autocorrelation among these "follows," I wish to utilize a GAM-GEE approach similar to that of Pirotta et al. 2011, using packages 'yags' and

Dear Daniel and Jorge, Thank you very much and it does help. If I have a string "ABCD", how can I access the second element of the string "B"? Thanks, Miao 2012/7/27 Daniel Nordlund <djnordlund@frontier.com> > > -----Original Message----- > > From: r-help-bounces@r-project.org [mailto:r-help-bounces@r-project.org] > > On Behalf Of jpm miao

Dear all, for demonstration purposes I want to display the basis functions used by a thin plate regression spline in a gamm model. I've been searching the help files, but I can't really figure out how to get the plots of the basis functions. Anybody an idea? Some toy code : require(mgcv) require(nlme) x1 <- 1:1000 x2 <- runif(1000,10,500) fx1 <- -4*sin(x1/50) fx2 <-

Hi, I have a table which I can read via read.csv: fx1<-read.csv(file="A_FX_M.csv", header=TRUE) TIME REER NTD JPY GBP HKD 1 198001 124.26 36.030 237.96 2.263980 4.8366 2 198002 126.59 36.030 244.05 2.290426 4.8765 3 198003 128.33 36.026 248.62 2.206045 4.9960 4 198004 127.85 36.063 251.67 2.215330 4.9760 5 198005 124.40 36.050 228.35 2.302026 4.8891 6 198006

Hello - I am having a problem with the function 'integrate'. I am running R on OSX (R 2.5.1). I am trying to suppress the error message when 'integrate' attempts to integrate across a parameter set giving a non-finite function value. I'm using it in a MCMC / simulated annealing algorithm, and it is entirely possible that some parameter sets will give non-

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 R users, I have been using the pdMat class "pdIdnot" (from the mgcv package)instead of "pdIdent" to avoid overflow in GLMM fits with the MASS package function glmmPQL, of the following form: fit1 <- glmmPQL(fixed=y0~-1+xx0, random=list(gp=pdIdent(~-1+zz0)), family=binomial) # vulnerable to overflow fit2 <- glmmPQL(fixed=y0~-1+xx0,

Hi all, I am not used to using the computer to do calculus and have up to now done my differentiation "by hand" , calling on skills I learned many years ago and some standard cheat sheets. My interest at present is in getting the second derivative of a gaussian, which I did by hand and results in a somewhat messy result involving terms in sigma^5 .. I have done some spot checks

Hi all, I've been attempting to fit a logistic glmm using glmmPQL in order to estimate variance components for a score test, where the model is of the form logit(mu) = X*a+ Z1*b1 + Z2*b2. Z1 and Z2 are actually reduced rank square root matrices of the assumed covariance structure (up to a constant) of random effects c1 and c2, respectively, such that b1 ~ N(0,sig.1^2*I) and c1 ~

I've been trying to get comfy with arima() and associated functions in the ts() package. I'm thinking seriously about using this package, and R generally, in a 4th year intro time series course that I'm teaching this autumn. I have a couple of questions about arima: (1) The help file says that residuals component of the value returned by arima() consists of the

Dear all, apologies for this (perhaps recurrent) question but I did not found a question when searching mailing lists. How to write a list of a simple kind, e.g.: abc <- list(one=(1:2), two=(1:5)) # to a file? I understand that write() & co. cannot work but when I try sink("aa.txt", append=T, split=T) abc sink() # the output is indeed "split by rows" in the

Dear Ilia, this page summarize all the possible options of the "nouveau" kernel module: https://nouveau.freedesktop.org/wiki/KernelModuleParameters/ I will test the options that you have suggested as soon as possible. Try a PCI video card is a very good advice because the AGP/FX1(chipset) is on the top of the "usual suspects" list; unfortunately the c8000 has 3.3 PCI slots

The very same option list, but for the "radeon" kernel module: http://xorg.freedesktop.org/wiki/RadeonFeature/#index4h2 Simone Mannori - Italy On 4 April 2016 at 05:56, Simone Mannori <simone.mannori at gmail.com> wrote: > Dear Ilia, > > this page summarize all the possible options of the "nouveau" kernel module: > >

<!doctype html public "-//w3c//dtd html 4.0 transitional//en"> <head> <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> <meta name="Author" content="sam"> <meta name="GENERATOR" content="Mozilla/4.61 [en] (Win98; I) [Netscape]"> <title>fgfg</title>

<!doctype html public "-//w3c//dtd html 4.0 transitional//en"> <head> <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> <meta name="Author" content="sam"> <meta name="GENERATOR" content="Mozilla/4.61 [en] (Win98; I) [Netscape]"> <title>fgfg</title>

Dear "nouveau" developers, I know that many very competent guys have already spent a lot of time and efforts on this issue without success. I have started to play with "hppa" two weeks ago and, with the support of the linux-parisc mailing list people, now I have a - almost fully - working workstation (hp c8000). Everythings work perfectly BUT the video card. No matter the

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems