similar to: rcspline.problem

Dear all, As I am new to the R community - although eager to advance- I would like to pose a question to the community. I have an SPSS file which I have imported it in R (with the read.spss command) which conists of scale (continuous) variable "adiponectin" and the corresponding categorical value "death" (0=No, 1=Yes). In all there are 60 observations (among which

Hello All, Although I have downloaded some libraries such as multivariate data analysis library (multiv) and ade4, their functions such as pca or reconst are not recognised. Should I install any thing else or use any instruction so that R could find the location of libraries? Thanks for your quick response, Haleh Yasrebi

Hello, Like QR decomposition, I am looking for decomposition to get U and V matrix of SVD decomposition (according to LINPACK, X = UDV'). Do you know if there is a function which could calculate this decomposition? Look forward to your reply, Haleh

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

Hello, Why do I get doubt (NA) in the factor of test classification even if I fix l (minimum vote)? By setting l, no doubt should be occurred. Look forward to your reply Haleh

Hello, I wanted to install MergeMaid package in v 2.2.1. I could install it but couldn't use without its dependant, Biobase. at biobase installation, I got the following error message In method for function "split": expanding the signature to include omitted arguments in definition: drop = "missing" Error in .MakeSignature(new("signature"), def, signature) :

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

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),

Hello, It might be an easy question but if you have many variables to fit in the lm function, how do you take all without specifying var1+var2+...+var2100 in the terms parameter in response ~ terms? Cheers, Carol

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

Hi folks, Using R Version 1.6.2 (2003-01-10) on SuSE Linux 7.2, I just installed Hmisc_1.5-3.tar.gz and Design_1.1-5.tar.gz These were taken from http://hesweb1.med.virginia.edu/biostat/s/library/r Checked the dependencies: Hmisc: grid, lattice, mva, acepack -- all already installed Design: Hmisc, survival -- survival already installed, so installed Hmisc first All seems to go

I have been attempting to use cph models to get better calibration of my models for which I had originally used logistic regression. I tried running with 40 repetitions and got an error. I then tried 500 repetitions (thinking that the NaNs in the output below might be caused by that choice) and then let my computer crunch for several hours and got only the same error message and

Please forgive me if you already received this. I had an e-mail sending glitch this morning. http://cran.r-project.org/src/contrib/checkSummary.html reported an error in Design.trans.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

Hello All, I would like to know what fit criterion (abstol arg) is in nnet. Is it the threshold for the difference btw the max output and target values? Is the value at each iteration also the difference btw max of output and target values over all output units (case of multiple classes)? How could value displayed at each iteration be related to SSE and abstol be related to threshold SSE,