search for: seeto

...p[, i]]) : # NA/NaN/Inf in foreign function call (arg 1) Gls(y ~ x, data=d, correlation = corARMA(p=1)) #This works I would rather use Gls than gls so that I can represent a variable with a spline using rcs. I'm using version 3-5.0 of rms in R 2.15.0. Thanks for any help you can give. Mark Seeto

...s(y ~ x1 + I(x1^2), data=d) Error in if (!length(fname) || !any(fname == zname)) { : missing value where TRUE/FALSE needed ols(y ~ I(x1^2), data=d) # No error message, but lacks linear term Is there a way to do these things without first creating new variables in the data frame? Thanks, Mark Seeto National Acoustic Laboratories

...from an ols object (using rms package)? library(rms) x <- rnorm(10) y <- x + rnorm(10) ols1 <- ols(y ~ x) Typing "ols1" displays adjusted R^2 among other things, but how can I assign it to a variable? I tried str(ols1) but couldn't see where to go from there. Thanks, Mark Seeto

...39;t use Nabble), with cc to the mailing list? If I choose the option to reply to the person by email, I don't see an option to cc to the mailing list. If I reply to the list, there's an option to email the post to someone, but I don't know the person's email address. Thanks, Mark Seeto -- View this message in context: http://r.789695.n4.nabble.com/Replying-on-Nabble-tp3691214p3691214.html Sent from the R help mailing list archive at Nabble.com.

...quot;a.-5" d x a.-5 1 0 0 2 1 1 d <- data.frame(d, y = c(0,1)) d x a..5 y 1 0 0 0 2 1 1 1 names(d)[2] <- "a.-5" d x a.-5 y 1 0 0 0 2 1 1 1 Why does the "a.-5" column name change to "a..5" when another column is added? Thanks, Mark Seeto

...to what Frank Harrell says should happen. I am either misunderstanding something or making a mistake in the code (I'm almost 100% certain I'm not mixing up the axes). I would be most appreciative if anyone could explain where I'm going wrong. Thanks for any help you can provide. Mark Seeto

...error happens if x.knots is simply defined as a vector like c(-1, 0, 1) (i.e. not using quantile). Is this the intended behaviour? The requirement that x.knots be in the workspace seems strange, given that the knot locations are stored in ols1$Design$parms. Thanks for any help you can give. Mark Seeto National Acoustic Laboratories, Australia

...nts: (Intercept) x3 x4 -0.1087 0.2830 0.1024 How can I specify the formula in the line marked *** so that the output will show "formula = y ~ x3 + x4" instead of "formula = paste..."? Thanks for any help you can give. Regards, Mark -- Mark Seeto Statistician National Acoustic Laboratories A Division of Australian Hearing

...live with this R^2 optimism? It can be decreased by taking a bigger penalty, but then the corrected R^2 is reduced. Also, a penalty of 9 gives a corrected slope of about 1.17 (corrected slope of 1 is achieved with a penalty of about 1 or 2). Thanks for any help/advice you can give. Mark -- Mark Seeto Statistician National Acoustic Laboratories A Division of Australian Hearing

...4 0.047 0.899 I don't understand the warning message, but for these examples, the ordinary t interval appears to be better than the bootstrap BCA interval. I would really appreciate any recommendations anyone can give on when bootstrap confidence intervals should be used. Thanks, Mark -- Mark Seeto National Acoustic Laboratories, Australian Hearing

...nd predict(fmi,d) to be the same except for the final NA, but they are not. Instead, fmi$fitted is the same as ols.imp2$fitted. Also, anova(ols.imp2) anova(fmi) unexpectedly give the same result in the ERROR line. I would be most grateful if anyone could explain this to me. Thanks, Mark -- Mark Seeto Statistician National Acoustic Laboratories <http://www.nal.gov.au/> A Division of Australian Hearing

...ion. A related question: does ols allow "y ~ ." notation? lm(y ~ ., data = d[, 2:4]) # This works ols(y ~ ., data = d[, 2:4]) # Gives error Error in terms.formula(formula) : '.' in formula and no 'data' argument Thanks for any help you can give. Regards, Mark -- Mark Seeto Statistician National Acoustic Laboratories <http://www.nal.gov.au/> A Division of Australian Hearing 126 Greville Street Chatswood NSW 2067 Australia P +61 2 9412 6800 F +61 2 9411 8273

...4 bit versions of R 2.14.2, and it also happens with ggplot2 plots. I'm using RGui on Windows 7. It did not happen with R 2.13.1. It's not a major problem, because the plots still appear to be produced correctly, but if anyone can tell me how to fix it, I'd appreciate it. Thanks, Mark Seeto

Hello, I'm fairly new to R and I can't work out how to produce a double inequality like (LaTeX) $0 \leq x \leq 1$ in the legend of a graph. If I try > legend(50, 0.1, legend = c(expression(0 <= x <= 1), c(2 <= x <= 3)), pch = c(1,1), col = c(2, 3)) then I get an error message "unexpected '<=' in ...". I've checked the help files for plotmath and