similar to: rcspline.plot query

Dear Dr ; I am a PhD student at Epidemiology department of National University of Singapore. I used R command (rcspline.plot) for plotting restricted cubic spline ??? the model is based on Cox. I managed to get a plot without adjustment for other covariates, but I have a problem regarding to adjusting the confounders. I applied below command to generate the matrix for

I would like to impute missing data in a set of correlated variables (columns of a matrix). It looks like transcan() from Hmisc is roughly what I want. It says, "transcan automatically transforms continuous and categorical variables to have maximum correlation with the best linear combination of the other variables." And, "By default, transcan imputes NAs with "best

I have an attached data.frame with a binary composite variable (comp) and a predictor (egfr). I am trying to plot a restricted cubic spline using: require (Hmisc) require (Design) rcspline.plot(x=egfr,y=comp,model="logistic",nk=5,knots=NULL, xlab="label1", ylab="Label2", plim=c(0,1), plotcl=TRUE, showknots=TRUE,

Hi all, I have a vector of proportions (post_op_prw) such that >summary(amb$post_op_prw) Min. 1st Qu. Median Mean 3rd Qu. Max. NA's 0.0000 0.0000 0.0000 0.3985 0.9134 0.9962 1.0000 > summary(cut2(amb$post_op_prw,0.0001)) [0.0000,0.0001) [0.0001,0.9962] NA's 1904 1672 1

Hi, I'm trying to get the 'xts' library from CRAN packaged for Fedora 12, and one of its dependencies is 'Hmisc'. I have created a spec file for that using R2spec, but building it fails with something in building the manual pages: Warning: ./man/escapeRegex.Rd:22: unknown macro '\backslash' [similar warnings snipped] Warning: newline within quoted string at

http://r.789695.n4.nabble.com/file/n3058505/file.csv file.csv Hi, Rusers I have a problem in making a rcspline.plot with a Hmisc package. My data is in the upload attachment. My programme as follows: library(Hmisc) A<-read.csv("file.csv",header=TRUE) attach(A)

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.

Dear R-help, I don't understand why Gls gives me an error when trying to fit a model with AR(2) errors, while gls (from nlme) does not. For example: library(nlme) library(rms) set.seed(1) d <- data.frame(x = rnorm(50), y = rnorm(50)) gls(y ~ x, data=d, correlation = corARMA(p=2)) #This works Gls(y ~ x, data=d, correlation = corARMA(p=2)) # Gives error # Error in

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 the model in an 'expanded' format-- probably what rcspline.restate()

Apologies if this is this too obscure for R-help. In package splines, ns(x,,knots,intercept=TRUE) produces an n by K+2 matrix N, the values of K+2 basis functions for the natural splines with K (internal) knots, evaluated at x. It does this by first generating an n by K+4 matrix B of unconstrained splines, then postmultiplying B by H, a K+4 by K+2 representation of the nullspace of C (2 by K+4),

Dear list, I'm currently analyzing some count data using a hurdle model. I've used the rcspline.eval function in the Hmisc-library to contruct the spline terms for the regression model, and what I want in the end is the ability to compute coefficients and confidence intervals for different changes in the smooth function as well as plotting the smooth function along with the

Dear R-users, I'm currently having trouble with the implementation of a groupedData object in the lme() function. Executing the following function > applyScalingSimp <- function(input.population) > { > ## GA is a time value > varInOrder <- c("GA","weight","grouping","sex") > modelVar <-

Hello, I am trying to reformat some data so that it is organized by group in the columns. The data currently looks like this: group X3.Hydroxybutyrate X3.Hydroxyisovalerate ADP 347 4 4e-04 3e-04 5e-04 353 3 5e-04 3e-04 6e-04 359 4 4e-04 3e-04

Am 22.07.2015 um 17:41 schrieb Windsor Dave (AdP/TEF7) <Dave.Windsor at us.bosch.com>: > Sorry for the top post, Outlook defaults strike again..... Outlook forces you to write above ? :-) -- LF

Hi, Using package gclus in R, I have created some graphs that show the trends within subgroups of data and correlations among 9 variables (v1-v9). Being interested for more details on these data I have produced also the var-covar matrices. Question: From a pair of two subsets of data (with 9 variables each, I have two var-covar matrices for each subgroup, that differ for a treatment on one

I was looking at the manpage for rsyslog.conf, primarily because I need to filter my manager's new fedora 22 logs coming to our loghost, because of the bug that I forwarded (if it gets through). At any rate, I am surprised: under selectors, I see that " The keywords error, warn and panic are deprecated and should not be used anymore." Huh? If I only want warn or more severe,

Hi I have attempted to install Samba on an ageing Motorola running SVR3. Unfortunately it has a limitation of 14 characters for filenames, which means that when I untar the package, many of the files are not created. Does anyone have a solution to this problem, or am I stuck? Many thanks -- Pete ---------------------------------------------------------------- morph@morhp.dircon.co.uk

Dear all, I am using function glmer from package lme4 to fit a generalized linear mixed effect model. My model is as follows: model1 <- glmer(fruitset ~ Dist*wire + (1|Site), data, binomial) summary(model1) Generalized linear mixed model fit by the Laplace approximation Formula: fruitset ~ Dist * wire + (1 | Site) Data: data AIC BIC logLik deviance 68.23 70.65 -29.11 58.23 Random

Just thought I'd pass along an error I received this weekend after upgrading two servers to 5.6. Both are HP Proliant DL380 servers running 64-bit, but one is a G6 model and one is a G7. After the upgrade and the reboot to the 2.6.18-238.9.1.el5 kernel, both servers displayed the following error at boot: "pci_mmcfg_init marking 256MB space uncacheable." Some Googling

Hi all, How can I set open files that contain values of Zero =0? These are valid values for the parameters I need to evaluate. I have tried CSV and tab formats. Trying XL Connect and/or XLConnectJars dies not seem to work to open Excel files so I am at a loss on how to get the data into a DF. Sample of data with 0 values: Filename Dur TBC Fmax Fmin Fmean Fc S1 Sc Pmc g8221843.13#