Displaying 20 results from an estimated 7000 matches similar to: "datadist options error, DESIGN library"
2008 Sep 10
2
relsurv package
Dear R-users,
I have a couple of questions about the relsurv package:
1) when I try to run the example:
fit <- rsmul(Surv(time,cens)~sex+as.factor(agegr)+ratetable(age=age*365.24,sex=sex,year=year),ratetable=slopop,data=rdata)
with the datasets in the package (rdata and slopop) it gives me an error:
Error in nrow(x) : object "x" not found
2)
If I have a date format
2009 Jul 09
2
datadist() in Design library
Hi I got an error message using datadist() from Design package:
> library(Design,T)
> dd <- datadist(beta.final)
> options(datadist="dd")
> lrm(Disease ~ gsct+apcct+rarct, x=TRUE, y=TRUE)
Error in eval(expr, envir, enclos) : object "Disease" not found
All variables inclduing response variable "Disease" are in the data frame
2008 Apr 17
1
Error in Design package: dataset not found for options(datadist)
Hi,
Design isn't strictly an R base package, but maybe someone can explain
the following.
When lrm is called within a function, it can't find the dataset dd:
> library(Design)
> age <- rnorm(30, 50, 10)
> cholesterol <- rnorm(30, 200, 25)
> ch <- cut2(cholesterol, g=5, levels.mean=TRUE)
> fit <- function(ch, age)
+ {
+ d <- data.frame(ch, age)
+
2007 Jul 25
2
Subscript out of bounds when using datadist() from Design library
I am running R version 2.4.1 on Windows XP. I have a question regarding the datadist() function from the Design library. I have a data.frame (call it my.data) with 4 columns. When I submit the code
datadist(data=my.data)
I get the following error message:
Error in X[[1]] : subscript out of bounds
I suspect there may be something wrong with my data.frame (I'm certain there is nothing
2011 Apr 12
1
Datadist error
Dear all,
I have performed a simple logistic regression using the lrm function from
the Design library. Now I want to plot the summary, or make a nomogram. I
keep getting a datadist error: options(datadist= m.full ) not created with
datadist.
I have tried to specify datadist beforhand (although I don't know why it
should be done):
ddist<-datadist(d) ##where d is my dataset
2007 Aug 27
2
validate (package Design): error message "subscript out of bounds"
Dear R users
I use Windows XP, R2.5.1 (I have read the posting guide, I have
contacted the package maintainer first, it is not homework).
In a research project on renal cell carcinoma we want to compute
Harrell's c index, with optimism correction, for a multivariate
Cox regression and also for some univariate Cox models.
For some of these univariate models I have encountered an error
2011 Aug 25
1
survplot() for cph(): Design vs rms
Hi, in Design package, a plot of survival probability vs. a covariate can be generated by survplot() on a cph object using the folliowing code:
n <- 1000
set.seed(731)
age <- 50 + 12*rnorm(n)
label(age) <- "Age"
sex <- factor(sample(c('male','female'), n, TRUE))
cens <- 15*runif(n)
h <- .02*exp(.04*(age-50)+.8*(sex=='Female'))
dt <-
2007 Oct 19
1
plot.Design
Dear R-users:
I am trying to use the following code to reproduce the figures on page 140
of Prof. Frank Harrell's book 'Regression Modeling Strategies':
rm(list=ls())
options(width=128)
library(Hmisc)
library(Design)
getHdata(counties)
counties$older <- counties$age6574 + counties$age75
label(counties$older) <- '% age >= 65, 1990'
counties$pdensity <-
2008 Mar 18
0
Trouble fitting models with nonNULL datadist in Design Library
Having trouble with datadist in the Design library. When I set the
datadist option, all of the fitting functions give error messages
similar to "Error in exists(x, envir, mode, inherits) : invalid first
argument". Code below was run in R version 2.5.0 under windows, but get
similar results with R version 6.2, except with additional problem that
it will not fit model with
2005 Jul 11
1
validation, calibration and Design
Hi R experts,
I am trying to do a prognostic model validation study, using cancer
survival data. There are 2 data sets - 1500 cases used to develop a
nomogram, and another of 800 cases used as an independent validation
cohort. I have validated the nomogram in the original data (easy with
the Design tools), and then want to show that it also has good results
with the independent data using 60
2009 Apr 29
1
Error with Design.Function(fit)
Hi all,
I'm reposting this with a more appropriate subject.
Do I need to define limits as the error message seems to suggest? If so, how? The error message, my code, the output and the first few lines of my data are all below.
Thank you!
"Error in Getlim(at, allow.null = TRUE, need.all = TRUE) :
variable dmodel.df does not have limits defined in fit or with datadist"
My code:
2008 May 29
2
Troubles plotting lrm output in Design Library
Dear R-helpers,
I'm having a problem in using plot.design in Design Library. Tho
following example code produce the error:
> n <- 1000 # define sample size
> set.seed(17) # so can reproduce the results
> age <- rnorm(n, 50, 10)
> blood.pressure <- rnorm(n, 120, 15)
> cholesterol <- rnorm(n, 200, 25)
> sex <-
2006 Apr 11
1
cph() in Design package
Hi there, I encountered a weird problem using cph()
with Design package:
I have 2 datasets, say "dat1" and "dat2", both data
frames with 3 columns "time","status" and "scores",
all numeric
If I run the following:
dd<-datadist(dat1)
options(datadist='dd')
dd
time status scores
Low:effect 37.0 0 -6.018
2004 Dec 15
2
using Hmisc and Design library
Hi, I encountered a weird problem when using the
Design and Hmisc problem. I have 2 data frame called
"a" and "b", both have 3 columns: "time", "status" and
"scores", a sample of the data frame is like:
data frame "a":
time status scores
1 21 1 99.61
2 38 0 101.11
3 51 0 100.62
4 48 0 87.52
5 78 0
2010 Oct 04
2
i have aproblem --thank you
dear professor:
thank you for your help,witn your help i develop the nomogram successfully.
after that i want to do the internal validation to the model.i ues the bootpred to do it,and then i encounter problem again,just like that.(´íÎóÓÚerror to :complete.cases(x, y, wt) : ²»ÊÇËùÓеIJÎÊý¶¼Ò»Ñù³¤(the length of the augment was different))
i hope you tell me where is the mistake,and maybe i have
2007 Oct 02
1
Design package: plot summary
Hi everybody,
I am a new user of R, design package.
I am trying to plot the estimated hazard ratio's of my cox regression
model with the confidence intervals. But I keep getting the
message:Error in `contrasts<-`(`*tmp*`, value = "contr.treatment") :
contrasts can be applied only to factors with 2 or more levels
I have dichotomous, categorical as well as continuous
2009 Apr 14
1
Function call error in cph/survest (package Design)
Dear UseR,
I do not know if this a problem with me, my data or cph/survest in package
design. The example below works with a standard data set, but not with my
data, but I cannot locate the problem.
Note that I am using an older package of survival to avoid a problem with
the newly renamed function in survival meeting Design.
Dieter
# First, check standard example to make sure
library(Design)
2012 Nov 07
1
error message no lines available in input
Dear all,
a program that worked well for weeks today gave me consistently the error
message
no lines available in input
referring to the lines
for (i in (0:(timeintervals-1)))
{ j=subjectquantity+6+i*(subjectquantity+7)
print (j)
results<-read.table(file, header=F, skip=j, nrows=subjectquantity)
timeintervals have been specified as argument as 12
subjectquantity as 9
the file it should
2008 Apr 03
1
Design package lrm summary and factors
Hello, I have question regarding the lrm function and estimating the odds
ratio between different levels of a factored variable.
The following code example illustrates the problem I am having. I have a
data set with an outcome variable (0,1) and an input variable (A,B,C). I
would like to estimate the effect of C vs B, but when I perform the summary
I only get A vs B and A vs C, even though I
2007 Sep 27
1
R: anova.Design
Dear All:
I tried to replicate a case study described by Prof. Harrell in Chapter 7 of
his Regression Modeling Strategies book, but failed on using anova.Design to
reproduce his table 7.1, Following is the code:
rm(list=ls())
library(Hmisc)
library(Design)
getHdata(counties)
counties$older <- counties$age6574 + counties$age75
label(counties$older) <- '% age >= 65, 1990'