similar to: Branching on 'grep' returns...

I'm trying to use boot.stepAIC for feature selection; I need to be able to specify the name of the dependent variable programmatically, but this appear to fail: In R-Studio with MS R Open 3.4: library(bootStepAIC) #Fake data n<-200 x1 <- runif(n, -3, 3) x2 <- runif(n, -3, 3) x3 <- runif(n, -3, 3) x4 <- runif(n, -3, 3) x5 <- runif(n, -3, 3) x6 <- runif(n, -3, 3) x7

Failed? What was the error message? Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Aug 22, 2017 at 8:17 AM, Stephen O'hagan <SOhagan at manchester.ac.uk> wrote: > I'm trying to use boot.stepAIC for

The error is "the model fit failed in 50 bootstrap samples Error: non-character argument" Cheers, SOH. On 22/08/2017 17:52, Bert Gunter wrote: > Failed? What was the error message? > > Cheers, > > Bert > > > Bert Gunter > > "The trouble with having an open mind is that people keep coming along > and sticking things into it." > -- Opus (aka

It seems that if you build the formula as a character string, and postpone the "as.formula" into the lm call, it works. instead of frm1 <- as.formula(paste(trg,"~1")) use frm1a <- paste(trg,"~1") and then strt <- lm(as.formula(frm1a),dat) regards, Heinz Stephen O'hagan wrote/hat geschrieben on/am 23.08.2017 12:07: > Until I get a fix that works, a

SImplify your call to lm using the "." argument instead of manipulating formulas. > strt <- lm(y1 ~ ., data = dat) and you do not need to explicitly specify the "1+" on the rhs for lm, so > frm2<-as.formula(paste(trg," ~ ", paste(xvars,collapse = "+"))) works fine, too. Anyway, doing this gives (but see end of output)" bst <-

Until I get a fix that works, a work-around would be to rename the 'y1' column, used a fixed formula, and rename it back afterwards. Thanks for your help. SGO. -----Original Message----- From: Bert Gunter [mailto:bgunter.4567 at gmail.com] Sent: 22 August 2017 20:38 To: Stephen O'hagan <SOhagan at manchester.ac.uk> Cc: r-help at r-project.org Subject: Re: [R] boot.stepAIC

OK, here's the problem. Continuing with your example: strt1 <- lm(y1 ~1, dat) strt2 <- lm(frm1,dat) > strt1 Call: lm(formula = y1 ~ 1, data = dat) Coefficients: (Intercept) 41.73 > strt2 Call: lm(formula = frm1, data = dat) Coefficients: (Intercept) 41.73 Note that the formula objects of the lm object are different: strt2 does not evaluate the formula. So

https://bugzilla.samba.org/show_bug.cgi?id=12569 Bug ID: 12569 Summary: Missing directory errors not ignored Product: rsync Version: 3.1.2 Hardware: All OS: Linux Status: NEW Severity: normal Priority: P5 Component: core Assignee: wayned at samba.org Reporter: axkibe

Hi, I have a PRI euroisdn link between an Alcatel PBX and Asterisk. I'm having some trouble with overlap dialing. Suppose I dial '874053' from an Alcatel extension ('7034') where '87' is an Alcatel prefix of type "ARS Prof.Trg Grp Seiz.with overlap". I'm expecting Asterisk to receive '1004053' (where '100' is a prefix which always shows

Hey, I'm using rsync to backup my server but there is a problem because rsync is using very huge amounts traffic. But first to the setup. The server I backup has 4GB of data and I use the following command to backup this data. /usr/bin/rsync -aze 'ssh -i /root/.ssh/backup.key -l backupuser' --rsync-path='sudo rsync' --delete --exclude-from=ex.list $SRC $TRG The problem is

I'm still learning how to program with R and I was hoping someone could take the time to show me how I can rewrite this code? Many thanks Tom data.intersects<-data.frame( x=c(0.230,0.411,0.477,0.241,0.552,0.230), y=c(0.119,0.515,0.261,0.431,0.304,0.389), angle=vector(length=6), length=vector(length=6),

Dear listers, I have two questions: (1) Is there a way in R to change the base-n of the calculations. I wnat to run some calculations either in binary (base-2) or base-4. Is there a way to specify that in R - to chnage from the decimal? (2) I also want to extract the digits from a larger number and store them as a vector. I could do it through converting top string, parsing the string and

hello...everybody... i have a question about mount point ... if i want to separate the streaming by category, let say i want to play an english category song. Then what should i do is just go to http://putera-trg.kustem/english and the song plays. if i want to play rock category songs, then just go to http://putera-trg.kustem/rock and the songs play. This is about a mount point but i still not

Hi, Im' using deliver with postfix (mailbox_command = /usr/lib/dovecot/deliver -n -m "$EXTENSION"). In dovecot.conf I've set mail_location to maildir:/var/spool/imap/users/%u/Maildir It works, but not for ONE "special" email. If I send mail to e.g. th at trg-oha.de it gets correctly delivered (see noerror.txt), but I have one mail with an attachment (i cannot

Hi! I am a relatively new user of R and I use it to prepare my dissertation. I have come to some very usefull and specific libraries written for S-PLUS 4 and would like to use them in R. Is that possible? I just found out that one of these libraries has already been transfered to R, while 3 others have not. For the matter of beeing more exact I''m interested in the dealing with missing

Dear listers, I must first say that all the suggestions and help I got from the list so far was of great help to my work. I approach you with another question. Is it possible to generate an automatic specification of the formula for lme (or other models in R)? Let me clarify. I'm using models in a simulation run and can have variable number of variables in the data. I would like to generate

Greetings. At the moment, I'm applying R to some AIX 'nmon' output, trying to get a handle on some disk performance metrics. In case anyone's interested: http://docs.osg.ufl.edu/tsm/pdf/ some of them are more edifying than others. (ahem) I'm trying to develop a somewhat general framework for plotting these measures, in the hopes that it's of some use to people other

Dear listers, First let me appologize if the same mail arrives multiple times. Recently I had some probelms sending my e-mails to the list. I encountered a problem when running glmmPQL procuedure doing multilevel modeling with a dichotomous outcome. Those are the two error messages I usually get: Error in logLik.reStruct(object, conLin) : NA/NaN/Inf in foreign function call (arg 3)

I'm working to automate the building of RPM packages for CRAN &c. In the process, I'm trying to get a sense of the correct dependencies. It's my sense that R CMD CHECK is the Right Way to check to see if a package is built properly. But according to Writing R Extensions, I could reasonably expect anything in Suggests or Imports to be used by CHECK. In RPM parlance, this turns

Dear R-devel, Is there a reason that lapply(NULL, ...) returns the empty list, rather than NULL? It seems intuitive to expect the latter, and rather counterintuitive that lapply(list(), ... ) returns the same value as lapply(NULL, ...). > lapply(list(), function(x) 1) list() > lapply(NULL, function(x) 1) list() > version _ platform i386-pc-mingw32 arch