similar to: GAM without intercept

Hi everybody, I am trying to generate two random walks with an specific correlation, for example, two random walks of 200 time steps with a correlation 0.7. I built the random walks with: x<-cumsum(rnorm(200, mean=0,sd=1)) y<-cumsum(rnorm(200, mean=0,sd=1)) but I don't know how to fix the correlation between them. With white noise is easy to fix the correlation using the function

Dear List, I've just tried to specify a GAM without an intercept -- I've got one of the (rare) cases where it is appropriate for E(y) -> 0 as X ->0. Naively running a GAM with the "-1" appended to the formula and the calling "predict.gam", I see that the model isn't behaving as expected. I don't understand why this would be. Google turns up this old

Hello, I've used this command to analyse changes in brain volume: mod1<-aov(Volume~Sex*Lobe*Tissue+Error(Subject/(Lobe*Tissue)),data.vslt) I'm comparing males/females. For every subject I have 8 volume measurements (4 different brain lobes and 2 different tissues (grey/white matter)). As aov() provides only type I anovas, I would like to use lmer() with type II, however, I have

Hi all, I constructed a GAM model with a linear term and two smooth terms, all of them statistically significant but the intercept was not significant. The adjusted r2 of this model is 0.572 and the deviance 65.3. I decided to run the model again without intercept, so I used in R the following instruction: regression= gam(dependent~ +linear_independent +s(smooth_independent_1)

Hey all, I'm using the gam() function inside the mgcv package to fit a penalised spline to some data. However, I don't quite understand what exactly the intercept it includes by default is / how to interpret it. Ideally I'd like to understand what the intercept is in terms of the B-Spline and/or truncated power series basis representation. Thanks!

On 10/29/2014 01:07 PM, Thomas Stein wrote: > About the --quiesce option. Do i need to do something inside the vm? The most > commonly would probably be a sql server running inside the vm. Do i need to > tell the sql server something about the --quiesce option i use? I read this > article here which suggests such a procedure. Okay, it's vmware, but... Is > that right? For

Dear All: I apologize for asking such an elementary question, but I could not find an adequate response on line. I am hoping to receive some help with the interpretation of the Intercept coefficient in the gam model below. I1 through I3 are dummy coded "Item difficulty" parameters in a data set that includes 4 items. If the Intercept is the value of Y when all other terms are 0, am I

Dear useRs, I am using mgcv version 1.7-16. When I create a model with a few non-linear terms and a random intercept for (in my case) country using s(Country,bs="re"), the representative line in my model (i.e. approximate significance of smooth terms) for the random intercept reads: edf Ref.df F p-value s(Country) 36.127 58.551 0.644

A customer is asking for a "manual". He's not talking about a How-To. He's talking about a PDF/DOC that shows what files do what and what parameters can be used and the syntaxis. I think Asterisk Business Edition has one that comes with the box. For the recent *not ABE* stable release will be nice. Comments? -- ------------------------------------------- Erick Perez Linux

Hi, Using version 1.0.rc15 I'm trying to let Evolution connect with IMAP over SSH. When run from the command line, this is what I get: ssh host dovecot --exec-mail imap * BAD [ALERT] imap binary must not be started from inetd, use imap-login instead. ssh host MAIL=/root/Maildir/ /usr/local/libexec/dovecot/imap * BAD [ALERT] imap binary must

Dear colleagues I found duplicate coordinates trying to convert a dataframe object (see the attached file) into a geodata one ('geoR' package is needed). My dataframe object consists of 2 columns with UTM geographical coordinates and a 3rd column with fish densities . To overcome the handicap I must add a tiny value generated by rnorm() to each geographical coordinate. And it must be in

Dear R-users, While making a prediction using the randomForest function (package randomForest) I'm getting the following error message: "Error in predict.randomForest(model, newdata = CV) : No forest component in the object" Here's my complete code. For reproducing this task, please find my 2 data sets attached ( http://www.nabble.com/file/p17855119/data.rar data.rar ).

Hi again! Please try something like this: > zone1.df zone sex len 1 1 male 15 2 1 fema 20 3 1 fema 17 4 2 fema 19 5 2 male 18 > str(zone1.df) 'data.frame': 5 obs. of 3 variables: $ zone: Factor w/ 2 levels "1","2": 1 1 1 2 2 $ sex : Factor w/ 2 levels "fema","male": 2 1 1 1 2 $ len : int 15 20 17 19 18 >

Dear all, I've recently switched from EViews to R with RMetrics/fSeries (newest version of july 10) for my analysis because of the much bigger flexibility it offers. So far my experiences had been great -prior I had already worked extensively with S-Plus so was already kind of familiar with the language- until I got to the fSeries package. My problem with the documentation of fSeries is that

Hi all, when I have c code like --- c code ------------- struct stest { /* deklariert den Strukturtyp person */ int age; float weight; } foo={44,67.2}; /* deklariert Variable des Typs person */ int main() { callAFunction(foo.weight); ------------------------ The generated c++API code to me seems to be too chatty in the sense that the foo.weight's data

Hello, I'm using aov() to analyse changes in brain volume between males and females. For every subject (there are 331 in total) I have 8 volume measurements (4 different brain lobes and 2 different tissues (grey/white matter)). The data looks like this: Subject Sex Lobe Tissue Volume subect1 1 F g 262374 subect1 1 F w 173758 subect1 1 O g 67155 subect1 1 O w 30067 subect1 1 P g 117981

Hi i have to table with IDs in each one. I want to make a join (as in sql) by the ID. Is any way to use the RODBC package (or other) in local tables (not a access, mysql, sql, etc. ) and made the join? Thanks in advance

Dear list, I need to transform the DateTime of my GPS data from: "666.1751" into "yyyy/mm/dd hh:mm:ss" I have the following code: d$Date <- ISOdatetime(2009, 1, 1, 0, 0, 0, tz = "GMT")+d$Date*(24*3600) This gives me: 2010-10-29 04:12:09, which is wrong. It should be 2010-10-29 06:12:09 Another example: 418.3219 corresponds to: 2010-02-23 07:43:30, but it

Hi, I'm trying to set up an iscsi 12.5 TB storage for some data backup. Doing so, I had some difficulties to find the right tool, maybe it's also a question of the system settings... The server is a 32Bit CentOS 5.3 with the recent updates. Ths iscsi connection can be establised. fdisk and parted fail to create any information on the device or fail completely. using the lvm tools

Dear all I have configure asterisk with 100 SIP PHONE ( SNOM ) but now thing is that my boss need phonebook feature find extention number by Pbook so i have read about it there is a feature in asterisk but it is with voicemail now i have IP SIP phone of SNOM so how to fine phone number by SIP phone ?? how to asterisk directory work ? Rgd satish patel