similar to: Names of variables needed in newdata for predict.glm

all.vars works fine, EXCEPT, it give a bit too much. I only want the regression variables, but in the following example I also get "k" the variable holding the chosen knots. Any machinery to find only "real" regression variables? cheers, Bendix library( splines ) y <- rnorm(100) x <- rnorm(100) k <- -1:1 ml <- lm( y ~ bs(x,knots=k) ) mg <- glm( y ~

Hi, Some try: > names(mi$xlevels) [1] "f" > all.vars(mi$formula) [1] "D" "x" "f" "Y" > names(mx$xlevels) [1] "f" > all.vars(mx$formula) [1] "D" "x" "f" When offset is indicated out of the formula, it does not work... Marc Le 07/03/2018 ? 06:20, Bendix Carstensen a ?crit?: > I would like

Hey I am trying to setup a dovecot cluster with 2 servers using replication /dsync. In front of it I got a Fortinet ADC (Load balance) and I think that I messing up the dsync. I see mails duplicated in the sync progress. If I disable one of the servers in the ADC, it seems to work and the sync if working without a problem. But if I use both servers with a round robin on the ADC, I see mailed

The example here puzzles me. It seems like the < operator doesn't work as expected. > l <- 0.6 > u <- seq(0.4, 0.7, 0.1) > u [1] 0.4 0.5 0.6 0.7 > mygrid <- expand.grid("l" = l, "u" = u) > mygrid l u 1 0.6 0.4 2 0.6 0.5 3 0.6 0.6 4 0.6 0.7 > mygridcollapsed <- mygrid[mygrid$l < mygrid$u, ] > mygridcollapsed l u 3 0.6 0.6 4

Dear all, I use Sweave to create my reports. my pdf output contains only NA's whenever my text contains special Danish characters like ø, ø and å. I have upgraded R from 2.13 to 2.15 and the scripts worked fine under R 2.13. Sweave("Test.Rnw", encoding = "utf8") texi2dvi("Test.tex", pdf = TRUE) %"Test.Rnw" \documentclass{article} \usepackage

Hello, Has anybody successfully called the Rmath library from C using the MS Visual Studio compiler (I am using Visual Studio 6.0)? I have compiled the Rmath library using gcc, and the 'test.c' program (which makes a call to qnorm) works fine when compiled with gcc. However, I get a fatal memory error when I run it after compiling it with Visual C. Would this memory problem be related

Email eller vedhæftet fil afsendt fra din adresse (eller med din adresse som afsender) er blevet afvist fra Allerød Kommune. Spam og virus bliver typisk sendt under dække af andre afsendere og den blokerede email behøver derfor ikke oprinde direkte fra dig. (Husk dog altid at have et opdateret antivirusprogram på din computer.) Du kan evt. scanne din computer med det gratis' værktøj

When I am using geom_bar I get these white slashes through the legend keys. I cannot figure out how to remove them. ggplot(diamonds, aes(clarity, fill=cut)) + geom_bar() I have tried using opts(legend.key = theme_blank()) but with no luck. Any suggestions would be much appreciated. I am using R vers. 2.15.0 and ggplot 0.9.1, win xp Best wishes Jonas Hal

Dear all, I would like to know how can I concatenate 2 data.frames into a single one. Both data frames have the same number of columns and the same class type in each correspondent column. So what I want is to have a new data.frame where I have first the values from one data.frame and then the values from a second data.frame would came after in this new data.frame. Thanks in advance. Med

Dear Sirs, I need to make a backup of email accounts (Maildir) a Linux server to a Windows server with NTFS, but get the following error: rsync: rename "/mnt/rsync/cur/137D94_29.ml110c.qnet.com.pe,S=32828:2,S.M92y5w" -> "cur/137D94_29.ml110c.qnet.com.pe,S=32828:2,S": Invalid argument (22) I have installed Rsync v3.0.9 The command is executed: rsync -rltDvu

If I have a contourplot (in the lattice package) and I want to add straight lines to it, how do I do this? I see that there are llines() and lsegement() functions for lattice plots, but they don't seem to do anything in this case: library(lattice) library(KernSmooth) x=rnorm(10000) y=x+rnorm(x,0,.5) a=bkde2D(cbind(x,y),.7) z=as.vector(a$fhat) grid=expand.grid(x=a$x1,y=a$x2) grid$z=z

Good day everyone, Can anyone suggest an effective method to solve the following problem: I have 2 dataframes D1 and D2 as follows: D1: dates ws wc pwc 2005-10-19:12:00 10.8 80 81 2005-10-20:12:00 12.3 5 15 2005-10-21:15:00 12.3 3 15 2005-10-22:15:00 11.3 13 95 2005-10-23:12:00 12.3 13 2 2005-10-24:15:00 10.3 2 95 2005-10-25:15:00 10.3 2 2 D2:

Dear List-Members, since 3 weeks I have been heavily working on reproducing the results of an economic paper. The method there uses the numerical solution of an integral within nonlinear least squares. Within the integrand there is also some parameter to estimate. Is that in the end possible to implement in R [Originally it was done in GAUSS]? I'm nearly into giving up. I constucted an

Hello rsync folks I'm trying to use rsync for backing up our servers. This mostly works extremely well, with no problems. However, 1 server is giving me a lot of trouble. It has a directory with (currently) 734088 files in it, and every time I try to backup this dir, rsync hangs after transferring roughly 2000 files. Sometimes it's around 1800, sometimes it's over 2100 (I think), but

Below is a SIMEX object that was generated with the "simex" function from the "simex" package applied to a logistic regression fit. From this mountain of information I would like to extract all of the values summarized in this line: .. ..$ variance.jackknife: num [1:5, 1:4] 1.684 1.144 0.85 0.624 0.519 ... Can someone suggest how to go about doing this? I can extract the

I'm trying to develop some graphic methods for glm objects, but they only apply for models where all predictors are discrete factors. How can I test for this in a function, given the glm model object? That is, I want something that will serve as an equivalent of is.discrete.glm() in the following context: myplot.glm <- function(model, ...) { if (!inherits(model,"glm"))

Hi, I've performed an lda and obtained a classification table for some of my data: > efa.dfa<-lda(groups~.,efa.scores.8,CV=T) > str(efa.dfa) List of 5 $ class : Factor w/ 2 levels "1","2": 1 2 1 2 1 1 2 2 1 2 ... $ posterior: num [1:160, 1:2] 0.99083 0.00852 0.93983 0.23186 0.85931 ... ..- attr(*, "dimnames")=List of 2 .. ..$ : chr [1:160]

library(effects) mod <- lm(Measurement ~ Age + Sex, data=d) e <-effect("Sex",mod) The effect is evaluated at the mean age. > e Sex effect Sex F M 43.33083 44.48531 > > e$model.matrix (Intercept) Age SexM 1 1 130.5859 0 23 1 130.5859 1 To evaluate the effect at Age=120 I tried: e

when using an nt 4.0 client, i can get non-guest access to all shares as i want. however, i cannot get access to any printers unless the printer is setup as a guest printer. nt 4.0 always seems to try to use the printer as guest user. nt 4 never gives me the chance to give a passwd for access to the server. (i do have SP3 on nt). if i use win 95 client, i can get access to all shares and

Greetings R experts, I'm having some difficulty recovering lda objects that I've saved within sublists using the save.image() function. I am running a script that exports a variety of different information as a list, included within that list is an lda object. I then take that list and create a list of that with all the different replications I've run. Unfortunately I've been