search for: rcs

Hi, I have fit a series of ols() models, by group, in this manner: l <- ols(y ~ rcs(x, 4)) ... where the series of 'x' values in each group is the same, however knots are not always identical between groups. The result is a table of 'coefs' derived from the ols objects, by group: group Intercept top top' top'' 1 6.864 0.01 2...

I was shocked to find no RCS on SXCE 107. I needed it to update some RCS archives I had copied over. No problem - go to sunfreeware and copy it over. No OpenSolaris branch? The Solaris 10 package should work. Only the package install appears to complete and does not - as documented only in a log file, not on my screen. OK...

I am encountering an error I do not know how to debug. The error arises when I try to add an interaction term involving two continuous variables (defined using rcs functions) to an existing (and working) model. The new model reads: model5 <- lrm( B_fainting ~ gender+ rcs(exactage, 7) + rcs(DW_nadler_bv, 7) + rcs(drawtimefrom8am, 7)+ DW_firsttime+ DW_race_eth + rcs(DW_BPS, 5)+ rcs(DW_BPD, 5)+ rcs(pulse3, 5)+ locat2+ (rcs(exactag...

...>= 65, 1990' pdensity <- logb(counties$pop.density+1, 10) label(pdensity) <- 'log 10 of 1992 pop per 1990 miles^2' counties <- cbind(counties, older, pdensity) # add 2 vars. not in data frame dd <- datadist(counties) options(datadist='dd') f <- ols(democrat ~ rcs(pdensity,4) + rcs(pop.change,3) + rcs(older,3) + crime + rcs(college,5) + rcs(income,4) + rcs(college,5) %ia% rcs(income,4) + rcs(farm,3) + rcs(white,5) + rcs(turnout,3), data=counties) incomes <- seq(22900, 32800, length=4) show.pts <- function(college.pts, income...

...$age6574 + counties$age75 label(counties$older) <- '% age >= 65, 1990' counties$pdensity <- log10(counties$pop.density+1) label(counties$pdensity) <- 'log 10 of 1992 pop per 1990 miles^2' dd <- datadist(counties) options(datadist='dd') f <- ols(democrat ~ rcs(pdensity,4) + rcs(pop.change,3) + rcs(older,3) + crime + rcs(college,5) + rcs(income,4) + rcs(college,5) %ia% rcs(income,4) + rcs(farm,3) + rcs(white,5) + rcs(turnout,3), data=counties) f an <- anova(f) and the error message reads as: Error in solvet(cov[idx, idx], coef[idx...

...$age6574 + counties$age75 label(counties$older) <- '% age >= 65, 1990' counties$pdensity <- log10(counties$pop.density+1) label(counties$pdensity) <- 'log 10 of 1992 pop per 1990 miles^2' dd <- datadist(counties) options(datadist='dd') f <- ols(democrat ~ rcs(pdensity,4) + rcs(pop.change,3) + rcs(older,3) + crime + rcs(college,5) + rcs(income,4) + rcs(college,5) %ia% rcs(income,4) + rcs(farm,3) + rcs(white,5) + rcs(turnout,3), data=counties) f r <- resid(f) windows(width=10.67, height=6.60) par(mfrow=c(3,3)) par(mar=c(4, 4, 2, 2))...

Dear, When fitting the following model knots <- 5 lrm.NDWI <- lrm(m.arson ~ rcs(NDWI,knots) I obtain the following result: Logistic Regression Model lrm(formula = m.arson ~ rcs(NDWI, knots)) Frequencies of Responses 0 1 666 35 Obs Max Deriv Model L.R. d.f. P C Dxy Gamma Tau-a R2 Brier 701 5e-07...

Dear all, Is the restricted cubic spline function working properly in the glm model? We used glm(y~rcs(x,5), family=binomial) but it seems that for some theoretical reasons the rcs, restricted cubic spline function can not be fitted by a glm function. Is this correct? Regards, Jan ((Originally, we used lrm(y~ rcs(x,5)) but we couldn't find how to derive the AIC value of the fitted model. Is t...

...859 -0.04884 -0.00589 -0.05392 ... ..- attr(*, "names")= chr [1:886393] "1" "2" "3" "4" ... $ means : num [1:16] 43.472 4.162 0.563 52.815 8.503 ... $ terms :Classes 'terms', 'formula' length 3 srv900 ~ rcs(Age, c(35, 50, 65)) + Sex + rcs(BL_HDL.A, 3) * rcs(BL_CHOLEST.A, 3) + rcs(BL_GGT.A, 3) + BL_TRIGLYC.A + rcs(BL_GGT.A, 3) + ... .. ..- attr(*, "variables")= language list(srv900, rcs(Age, c(35, 50, 65)), Sex, rcs(BL_HDL.A, 3), rcs(BL_CHOLEST.A, 3), rcs(BL_GGT.A, 3), BL...

...formulas for the ols function in the rms package. I want to have two variables represented as restricted cubic splines, and also include an interaction as a product of linear terms, but I get an error message. library(rms) d <- data.frame(x1 = rnorm(50), x2 = rnorm(50), y = rnorm(50)) ols(y ~ rcs(x1,3) + rcs(x2,3) + x1*x2, data=d) Error in if (!length(fname) || !any(fname == zname)) { : missing value where TRUE/FALSE needed ols(y ~ rcs(x1,3) + rcs(x2,3) + I(x1*x2), data=d) Error in if (!length(fname) || !any(fname == zname)) { : missing value where TRUE/FALSE needed I get the same err...

...ntity' link arguments, and finally, the error message when testing just an iter.max argument. This was run in R 2.13.0 with rms version 3.3-0 on Ubuntu Linux 11.04. Thank you in advance, and all insights and criticisms are appreciated. Mike library(rms) > f <- bj(Surv(ftime, stroke) ~ rcs(age,5) + hospital,x=TRUE, y=TRUE) > f <- bj(Surv(ftime, stroke) ~ rcs(age,5) + hospital, link='identity',x=TRUE, y=TRUE) No convergence in 50 steps Failure in bj.fit > f <- bj(Surv(ftime, stroke) ~ rcs(age,5) + hospital, link='identity', iter.max=200, x=TRUE, y=TRUE) Er...

Hi, I read data from a file. I'm trying to understand how to use Design.rcs by using simple test data first. I use 1000 integer values (1,...,1000) for x (the predictor) with some noise (x+.02*x) and I set the response variable y=x. Then, I try rcs and ols as follows: m = ( sqrt(y1) ~ ( rcs(x1,3) ) ); #I tried without sqrt also f = ols(m, data=data_train.df); print(f); [...

...Rd * checking Rd files ... ERROR Rd files with syntax errors: /var/mnt/hda3/R.check/r-devel/PKGS/Design/man/Design.trans.Rd: unterminated section 'alias' The .Rd file is attached. It begins with \name{Design.trans} \alias{Design.trans} \alias{asis} \alias{pol} \alias{lsp} \alias{rcs} \alias{catg} \alias{scored} \alias{strat} \alias{matrx} \alias{\%ia\%} \title{ Design Special Transformation Functions } \description{ This is a series of functions (\code{asis}, \code{pol}, \code{lsp}, .... No error is detected when I run CMD CHECK using the latest R-devel locally, using the 5...

Hi Apologies if this is a FAQ, I've tried Googling... I'm trying to set up CD-RCS on Windows boxes talking to our Linux Server (Samba 3.0.10-1.fc1.1.legacy). Apparently, the way this works is that it needs to be able to delete files that are read-only and not owned by you; they recommend that you set "delete readonly = yes" in the config file. The problem is that...

We just upgraded our server to run Mandrake 10.1. It is using kernel 2.6.8 and we are now using Samba 3.0.10. I am using the RCS commands on a Windows XP Professional computer. Before the upgrade, I had no problem checking files out with the "co -l" command no matter who owned the file. Now I get a "Permission denied" error unless I'm the owner of the RCS control file. The "delete readonly&q...

Hi, I am using the rcs() function in the Design library to model non-linearity that is not well characterized by an otherwise mechanistic function. I am able to make the model 'available' to others through the excellent nomogram() function and the set of tables that it can create. However, I would like to present...

...Network Theory Ltd -- Publishing Free Software Manuals 15 Royal Park Bristol BS8 3AL United Kingdom Tel: +44 (0)117 3179309 Fax: +44 (0)117 9048108 Web: http://www.network-theory.co.uk/ Index: src/library/base/man/Arithmetic.Rd =================================================================== RCS file: /cvs/R/src/library/base/man/Arithmetic.Rd,v retrieving revision 1.11 diff -u -r1.11 Arithmetic.Rd --- src/library/base/man/Arithmetic.Rd 2003/11/29 10:56:07 1.11 +++ src/library/base/man/Arithmetic.Rd 2003/12/12 15:44:47 @@ -30,7 +30,7 @@ by element operations. The elements of shorter vec...

I shouldn't try to do this on Monday morning. Can anyone suggest why the make for R 0.63 is failing for me under Solaris (SunOS 5.6). Paul Gilbert ... creating src/scripts/html2dos creating tests/Makefile creating tests/Examples/Makefile creating src/include/Platform.h R is now configured for sparc-sun-solaris2.6 Source directory: . Installation directory: /usr/local C

Hello again. Question as in topic. Thanks in advance.

Dear guRus, I am doing a logistic regression using restricted cubic splines via rcs(). However, the fitted probabilities should be nondecreasing with increasing predictor. Example: predictor <- seq(1,20) y <- c(rep(0,9),rep(1,10),0) model <- glm(y~rcs(predictor,n.knots=3),family="binomial") print(1/(1+exp(-predict(model)))) The last expression should be a non...