Displaying 20 results from an estimated 7000 matches similar to: "R-devel Digest, Vol 98, Issue 19"
2014 Jun 16
1
model.frame and parent environment
Someone has reported a problem with predict.coxph that I can't seem to solve. The
underlying issue is with model.frame.coxph; the same issue is also found in lm so I'll use
that for the example.
--------------------------
> test <- data.frame(y = 1:10 + runif(10), x=1:10)
> myfun <- function(formula, nd) {
fit <- lm(formula, data=nd, model=FALSE)
2014 Jun 17
0
model.frame question
Brian,
So let me ask an r-core opinion. Should I change the default to model=TRUE? Survival
is heavily used and there is something to be said for consistency within the central
packages. Sometimes old habits die hard, and there is a "save memory" part of me that
hates to save a large object that likely won't be used. Not nearly as relevant today as
when I started my
2011 Apr 18
2
help with eval()
I've narrowed my scope problems with predict.coxph further.
Here is a condensed example:
fcall3 <- as.formula("time ~ age")
dfun3 <- function(dcall) {
fit <- lm(dcall, data=lung, model=FALSE)
model.frame(fit)
}
dfun3(fcall3)
The final call fails: it can't find 'dcall'.
The relevant code in model.frame.lm is:
env <- environment(formula$terms)
2006 Sep 18
0
Permission denied
Hii All,
I'm really a linux newbie, I managed to run into a problem with
rsync. I want to conduct backups of various files and directories over
about miles of Internet, so I chose ssh and rsync for the job.
Earlier it used to work fine, but suddenly donno what
happened it started to show the following errors
sync: mkstemp
2002 Jun 27
1
Building from a source-code library under windows
Dear All,
I have a pair of .cpp and .def file can be compiled using VC++ and works
perfectly well in S-PLUS.
I wanted to do the same for R; so I followed the guidline given in "Building
from a source-code library under Windows" as much as possible and manage to
compile them using VC++ and call it from R. But it gives different answer
from the one called from S-Plus.
I know that I did
2004 Mar 01
0
se.contrast ....too hard??? .... Too easy????? .....too trivial???? ...... Too boring.....too????????
Hi all,
Regular and avid readers of this column will know that Don Driscoll and
I have recently posted two messages requesting assistance concerning an
apparent failure of "se.contrast" to produce an se for a contrast. So
far, an ominous silence rings in our ears, but read on Gentle Reader,
and see if even the machinations of "debug" doesn't stimulate you to
respond with a
2004 Feb 26
0
se.contrast ???????????
Hi all,
Just to follow up Don Driscoll's earlier post, can anyone please explain
why "se.contrast" fails here??
> shp<-factor(rep(c("reserve","strip"),each=96))
>
>
>
site<-factor(rep(c("1g","1p","1t","2g","2p","2t","3g","3p","3t","4g","4p
2006 Jun 26
0
[klibc 21/43] alpha support for klibc
The parts of klibc specific to the alpha architecture.
Signed-off-by: H. Peter Anvin <hpa at zytor.com>
---
commit 5e5ce29210ac33a0b3704eb9ab5e5d5b55375575
tree 2ec24df596e13c21b68da4d905f546770d36fdad
parent 8529b52550ba78984998d3a9cc9deb467217fa3e
author H. Peter Anvin <hpa at zytor.com> Sun, 25 Jun 2006 16:58:14 -0700
committer H. Peter Anvin <hpa at zytor.com> Sun, 25 Jun
2009 May 21
3
[LLVMdev] Passing a pointer to a function
I recently began hacking around with my first LLVM pass. The big
picture is that I would like to insert function calls for each
instruction type, and pass some parameters based on the instruction
type. Then I will link the output to some C file that implements those
functions.
Things were going well until I started trying to make function calls
with a pointer as a parameter. For example, I would
2011 Jan 06
1
Help spruce up a ggplot graph
Given the data structure below and the call to ggplot2, how can I increase the size of the axis scale points, the line weight, and the size of the legend?
ddata <-structure(list(year = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("2003", "2007"), class = "factor"),
area = structure(c(7L, 6L, 1L, 2L, 3L, 4L, 5L, 7L, 6L, 1L,
2009 May 08
2
Malformed HTML
I''m using Mechanize to parse an extraordinarily malformed html page.
After submitting a form like so:
page = mech.submit(dform)
The result I get back is truncated. I suspect that it''s because the
source HTML looks like this:
<html>
<head> yadda yadda</head>
<p>some text</p>
<html>
<table yadda yadda>
My
2004 Feb 24
0
(no subject)
G'day,
I'm fitting a simple one-way nested anova and would like to obtain standard
errors or confidence limits. I'm using se.contrast to try to get se for
the contrast between the two levels of the main effect, but I get an error
message (see below). What is going on?
Don
> shp<-factor(rep(c("reserve","strip"),each=96))
>
2004 Feb 24
0
se.contrast
G'day,
I'm fitting a simple one-way nested anova and would like to obtain standard
errors or confidence limits. I'm using se.contrast to try to get se for
the contrast between the two levels of the main effect, but I get an error
message (see below). What is going on?
Don
> shp<-factor(rep(c("reserve","strip"),each=96))
>
2011 Mar 15
1
Problem with nls.lm function of minpack.lm package.
Dear R useRs,
I have a problem with nls.lm function of minpackl.lm package.
I need to fit the Van Genuchten Model to a set of data of Theta and hydraulic conductivity with nls.lm function of minpack.lm package.
For the first fit, the parameter estimates keep changing even after 1000 iterations (Th)
and
I have a following error message for fit of hydraulic conductivity (k);
Reason for
2001 Oct 16
0
plot function
Hola!
It is somewhat inconvenient to use plot.function, when add=TRUE.
The following (miniscule) change makes it behave better:
plot.function <-
function (fn, from , to, xlim = NULL, ...)
{
if (!is.null(xlim)) {
if (missing(from))
from <- xlim[1]
if (missing(to))
to <- xlim[2]
}
curve(fn, from, to, xlim = xlim, ...)
}
The only
2009 May 24
1
filling area under a function line
Hi R collective,
I quite like the "curve" function because you can chuck a R function
into it, and see the graph in one line of R.
I had a google and found some threads for filling under the line;
http://tolstoy.newcastle.edu.au/R/e2/help/07/09/25457.html
However they seem to miss the point of the simplicity of going, "I
wonder what that looks like, and can I have some colour
2001 Apr 30
2
plotting an expression
I am sure it is just me not understanding how R works, but could somebody
explain why
curve(cos(x))
works and
curve(expression(cos(x))
does not?
I have done some investigating and here is what I found. If I comment out
the line of curve indicated below, both calls work fine.
function (expr, from, to, n = 101, add = FALSE, type = "l", ylab = NULL,
log = NULL, xlim =
2009 Dec 15
2
split.data.frame
Hello,
I very much enjoy "with" and "subset" semantics for data frames and was
wondering if we could have something similar with split, basically by
evaluating the second argument "with" the data frame :
split.data.frame
function(x, f, drop = FALSE, ...){
call <- match.call( )
fcall <- call( "with", data = call[["x"]],
2013 Oct 09
1
frailtypack
I can't comment on frailtypack issues, but would like to mention that coxme will handle
nested models, contrary to the statement below that "frailtypack is perhaps the only ....
for nested survival data".
To reprise the original post's model
cgd.nfm <- coxme(Surv(Tstart, Tstop, Status) ~ Treatment + (1 | Center/ID), data=cgd.ag)
And a note to the poster-- you should
2005 Nov 09
1
script/console
hi all,
on win xp (fr) when i start a console , i cant use any ctrl caracteres.
so i cant do nothing, as a = [ "a","b"] became a = "a" , "b"
even if I use cut/paste.
what the problem ?
"Ce message et toutes les pices jointes sont des informations
strictement confidentielles et rserves au(x) destinataire(s). Ce
courriel n''a pas de valeur