The read.table.ffdf function in the ff package can read in delimited files
and store them to disk as individual columns. The ffbase package provides
additional data management and analytic functionality. I have used these
packages on 15 Gb files of 18 million rows and 250 columns.
On Tuesday
I would like to incorporate multiple summary plots from the rms
package into a single device and to control the titles, and also to
open a new device when I reach a specified number of plots. Currently
I am only getting a single "plot(summary(" graph in the upper left-
hand corner of each successi
Is the labeling/naming of levels in the documentation for the
predict.glmmPQL function "backwards"? The documentation states "Level
values increase from outermost to innermost grouping, with level zero
corresponding to the population predictions". Taking the sample in
the documentation:
fit <- g
<- orig.df[,-c(1,3,5)]
> head(new.df,2)
num11.20 lc11.20 uc11.20
1 11 k K
2 12 l L
On Jan 6, 9:49 am, David Winsemius wrote:
> On Jan 6, 2012, at 10:00 AM, Mike Harwood wrote:
>
> > How does R do it, and should I ever be worried? I always
How does R do it, and should I ever be worried? I always remove
columns by index, and it works exactly as I would naively expect - but
HOW? The second illustration, which deletes non contiguous columns,
represents what I do all the time and have some trepidation about
because I don't know the mec
This issue occurs only when both the evd and ismev packages are loaded.
Please retract this posting, if possible. Thank you in advance!
Mike
-- Forwarded message --
From: Mike Harwood
Date: Mon, Dec 12, 2011 at 7:47 PM
Subject: tcplot documentation in evd package
To: r-help@r
Hello, and please advise regarding any errors/omissions on my part.
However, the documentation for R's tcplot function in the evd package
appears to contain an error. I am using evd version 2.2-4 in R 2.14.0
with Windows XP.
> data(portpirie)
> mrlplot(portpirie) ## No Error
> tlim <- c(3.6, 4.
Hello, and thank you in advance.
I would like to capture the expected survival from a coxph model for
subjects in an observational study with recurrent events, but the
survexp statement is failing due to memory. I am using R version
2.13.1 (2011-07-08) on Windows XP.
My objective is to plot the
Hello list, and thank you in advance.
I'm unable to generate predicted values when specifying newdata using
phreg and aftreg models
in the eha package, but I do not have the same problem with a
proportional hazards model from coxph. Without the newdata argument
the predicted values are returned,
I am unsure what is being returned, and what is supposed to be
returned, when using 'predict' with "type='expected'" for an aftreg
survival model. The code below first generates a weibull model, then
uses predict to create a vector of the linear predictors, then
attempts to create the 'expected'
Hello,
I would like to increase the number of iterations for running a
Buckley-James regression model in the rms package, but something is
apparently syntactically wrong. The following code initially is
exactly as it appears in the help page (which runs successfully), then
a "failure to converge"
ström
Depends: R (>= 2.2.0), survival, graphics
Maintainer: Göran Broström
Packaged: 2011-03-01 14:56:12 UTC; gb
Repository: CRAN
Date/Publication: 2011-03-01 15:50:52
Built: R 2.13.0; i386-pc-mingw32; 2011-04-15 08:22:36 UTC; windows
Mike
On Mon, May 2, 2011 at 10:38 AM, Mike Harwood wrote:
Hello,
I am apparently confused about the use of an id parameter for an event
history/survival model, and why the EHA documentation for aftreg does
not specify one. All assistance and insights are appreciated.
Attempting to specifiy an id variable with the documentation example
generates an "ove
os mingw32
system i386, mingw32
status
major 2
minor 12.2
year 2011
month 02
day25
svn rev54585
language R
version.string R version 2.12.2 (2011-02-25)
package ‘survival’ version 2.36-5
On Apr 20, 11:03
Hello,
I probably have a syntax error in trying to generate an expected
survival curve from a weighted cox model, but I can't see it. I used
the help sample code to generate a weighted model, with the addition
of a "weights=albumin" argument (I only chose albumin because it had
no missing values,
h.D.
> Statistical Data Center
> Intermountain Healthcare
> greg.s...@imail.org
> 801.408.8111
>
>
>
>
>
> > -Original Message-
> > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
> > project.org] On Behalf Of Mike Harwood
> > Se
G'morning
What does the error message "Error in x %*% coef(object) : non-
conformable arguments" indicate when calculating the response values
for
newdata with a model from bigglm (in package biglm), and how can I
debug it? I am attempting to do Monte Carlo simulations, which may
explain the loop
Hello, and thanks in advance!
What does the error message "Error in x %*% coef(object) : non-
conformable arguments" when indicate when predicting values for
newdata with a model from bigglm (in package biglm), and how can I
debug it? I am attempting to do Monte Carlo simulations, which may
expla
Hello,
Can the x and y axis limits be specified in a density plot with the
mixtools package for a finite mixture model? Uncommenting the xlim2/
ylim2 lines in the plot command below generates 'not a graphical
parameter' warnings (and does not change the axis settings), and
uncommenting the xlim/y
Is there a way to generate identical random samples using R's runif
function and SAS's ranuni function? I have assigned the same seed
values in both software packages, but the following results show
different results. Thanks!
R
===
> set.seed(6)
> random <- runif(10)
> random
[1] 0.6062683 0.9
Is there an alternative to par(mfrow=c(2,1)) to get stacked scatterplot
matrixes generated with "pairs"?
I am using version 2.11.1 on Windows XP. The logic I am using follows, and
the second "pairs" plot replaces the first plot in the current graphics
device, which is not what I expected (or desi
Marc,
Thank you for your help. I had tried "assign", but I did not enclose the
table name in quotes in my function call. I needed to see what someone else
had written before I ever would have noticed it!
On Fri, Feb 12, 2010 at 10:00 AM, Marc Schwartz wrote:
> On Feb 12, 2010, at 8:19 AM, mah
22 matches
Mail list logo