[R] mgcv 'bam' : prediction levels for random effects

Hi R users, I'm using the 'bam' function in mgcv to examine trends in a remotely sensed vegetation index. I have one random effect variable, 'cons' (with six levels) which identifies different subjects within this analysis. My model is specified as follows: rm4<-bam(trend~factor(zone)+s(cons,bs

Re: [R] Assign name to object for each iteration in a loop.

Thanks Michael! Yeah, that dealt with one of the problems. I still get the following error message: Error in cat(list(...), file, sep, fill, labels, append) : argument 1 (type 'list') cannot be handled by 'cat' I know that this has something to do with writing the names of the output file, an

[R] Assign name to object for each iteration in a loop.

Hi R-users, I'm trying to produce decompositions of a multiple time-series, grouped by a factor (called "area"). I'm modifying the code in the STLperArea function of package ndvits, as this function only plots produces stl plots, it does not return the underlying data. I want to extract the tren

[R] RTisean executable problem in STLperArea

Hi, I’m trying to use the STLperArea function in the ndvits package. I can run the sample data (“SLPSAs_full”) without any problem. However, when I come to run the function on my own data, I get the following message. Waiting to confirm page change... Error in .checkPath(path) : no TISEAN execu

[R] Compare two time series, at each time step

Hi, I have two time series of NDVI data, one inside a protected area, and one outside. I have decomposed the time series, to remove the seasonal component and residual variation. I want to compare the long-term trends inside and outside of the protected area over time. So, at each time step (1 y

[R] Var.calc in Matching

Hi, I already have matched samples, which were matched using different software. I need to calculate Abadie-Imbens standard errors, together with the average treatment effect. I know that the Matching package enables me to calculate these after a Matching procedure, but is there any way to do it o

[R] Combining multiple MatchIt outputs.

Hi All, Sorry if this seems obvious, I'm new to R. I've run nearest-neighbour matching in MatchIt for three different treatment groups. By this I mean, I have three separate treatments, each with a different control group. I want to combine the MatchIt outputs so that I can estimate treatment

[R] Fill in missing times in a timeseries with NA

Hi, I have a irregularly spaced time series dataset, which reads in from a .csv. I need to convert this to a regularly spaced time series by filling in missing rows of data with NAs. So my data, called NtuMot, looks like this (I've removed some of the additional rows for simplicity) ELEID