similar to: GAM function with interaction

Hello R Users, I have a question regarding fitting a model with GAM{mgcv}. I have data from several predictor (X) variables I wish to use to develop a model to predict one Y variable. I am working with ecological data, so have data collected many times (about 20) over the course of two years. Plotting data independently for each date there appears to be relationships between Y (fish density)

Hello all, Thank you for the previous assistance I received from this listserve! My current question is: How can I create an appropriate matrix of values from a GAM (actually a GAMM) to make a 3-D plot? This model is fit as a tensor product spline of two predictors and I have used it to make specific predictions by calling:

Hello, I am writing regarding log transformation of data in a single matrix column, and subsequent use of these data in a glm model fit. I have a data matrix in which I am using the log function to transform the values. This transformation results in -Inf values in some places, though. I then receive an error when this matrix is used in the glm function, and would like to know this

Hello, I have a question regarding how to subset/select parts of a data frame (matrix) in order to plot data associated only with this subset. Specifically I have a large data frame in which one column contains ID values (dates), and other columns contain data I would like to plot (temperature, light, etc.). I would like to break up this large matrix so as to plot data associated

Greetings! I am looking for advice regarding the best way to compare GAMMs. I know other model outputs return enough information for R's AIC, ANOVA, etc. commands to function, but this is not the case with GAMM unless one specifies the gam or lme portion. I know these parts of the gamm contain items that will facilitate comparisons between gamms. Is it correct to simply use these values

Hello! Thanks for reading this request for assistance. I have a question regarding creating a histogram-like figure from data that are not currently in the correct format for the "hist" command. Specifically, my data have been processed and are in a matrix with columns containing the variables of interest and separate columns containing the number of times this variable was

Hi, I'm using LVM based storage pools and I'm wondering if there is a way to specify several source-dev on the command line for creating a volume group spread over several devices : one device /dev/sdc1 is ok: * virsh pool-define-as --name lvmpool --type logical --source-dev /dev/sdc1 --source-name vg --target /dev/vg I would like something like (but sadly doesn't work): ? virsh

Thank you for answering. Yes sure I can from an xml, But is it possible from the command line with define-pool-as ? Best, Matt ----- Mail original ----- > De: "Humble Chirammal" <hchiramm@redhat.com> > À: "Matthieu Simonin" <matthieu.simonin@inria.fr> > Cc: libvirt-users@redhat.com > Envoyé: Mardi 19 Novembre 2013 11:06:35 > Objet: Re:

Hello,   I am currently working on a modeling project using Fortran to run large repetitive loops (many DO loops). As part of this process I would like to use a model fit in R and currently stored as an R object. This is a rather complex model, a GAMM, and I am curious whether there is a way to call this model from Fortran. I am not sure "call" is correct terminology, but I

Hello R users, I have a basic "computer programing" question. I am a student currently taking a course that uses Fortran as the main programming language, but the instructors are open to students using any language they are familiar with. I have used R previously, and am wondering if there is any benefit to my learning Fortran, or whether I should stick with R for this class. Any

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

Hi, I'm wondering why the spectrum() phase of quadrature couple isn't purely +/-pi. But mostly, I'm looking for a recommended way to take a 2-D spectrum and convert it into a single complex array. Kindly consider: # 10 Hz sine wave 10 seconds long sampled at 50 Hz deltaT = 1/50 t = seq(0, 10, deltaT) w = 2 * pi * 10 x = ts( sin( w * t ), deltat = deltaT ) y = ts( sin(

Dear R user, I have written a function which returns max,min and variation of a power (see below) Power is a given matrix(1,n) I call the function >Variation<-VAR(p,(n-deltat)) Now the problem is when I want plot(Results[1],Results[2]). Not possible! I become the following error (in english it means: Error in as.double.default(x) :Object cannot be transformed in double) >

Hi, I am getting unexpected behaviour from aggregate.ts(). The 'ndeltat' argument is effectively being reduced by 1 in some cases, even when it is an integer, with the result that the blocks to be aggregated are not of the expected size, and also that the end() of the aggregated series is much later than the end() of the original series. rawts <- ts(rep(1:10, each = 5), start = 1) ##

Hi I am still having problems with using the stl function, when I read the csv file into R into a file called tkr and use dim(tkr) the result is 132 1 which is fine. When coerce it into a trime series using ts either: tstkr <- ts(t(tkr), deltat=1/12) or tstkr <- ts(c(tkr), deltat=1/12) and use the stl function I get the following error: Error in

I am a little bit confusing about the following help message on how to fit a GAM model with interaction between factor and smooth terms from http://rss.acs.unt.edu/Rdoc/library/mgcv/html/gam.models.html: ?Sometimes models of the form: E(y)=b0+f(x)z need to be estimated (where f is a smooth function, as usual.) The appropriate formula is: y~z+s(x,by=z) - the by argument ensures that the smooth

Greetings! I have one quick question: How do you do a Tukey test on an ANCOVA? Thanks for any tips! -Paul

Hi. I have to forecast a time series of a Internet network traffic bitrate. The data are in file http://www.forumaltavilla.it/joomla/datitesi/dati.datand the sampling time is every 0.05 seconds. Now, i want to use HoltWinters forecasting. This is my script. dt=1.58443823e-9 #0.05 seconds in years dati.ts=ts(scan("dati.dat"),start=0,deltat=dt) model=HoltWinters(dati.ts)

Hello, I am attempting to create plots using two continuous variables and it seems I should be able to use the simple "plot" command. However, my x-axis values are dates, and I believe this could be the reason I am receiving error messages reading: Error in Summary.factor(c(2L, 4L, 3L, 5L, 6L, 7L, 1L), na.rm = FALSE) : range not meaningful for factors I have attempted to

I would like a sequence of dates with a time step of 15 minutes starting: 1/1/2006 00:00:00 - 12/31/2006 23:45:00 function(x) { chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x)) } this is the piece of code I use to read in zoo objects for any help I would be grateful I have tried sequence and I can not seem to get it to work -- Let's not