Displaying 4 results from an estimated 4 matches for "fmendolia".
Did you mean:
mendolia
2011 May 30
2
Group by multiple variables
Hello,
I would like to create a group variable that is based on the values of three variables:
For example,
> dat <- data.frame(A=c(1,1,1,1,1,2,2,2,2,2),
B=c(1,1,1,5,5,5,9,9,9,9),
C=c(1,1,1,1,1,2,2,7,7,7))
> dat
A B C
1 1 1 1
2 1 1 1
3 1 1 1
4 1 5 1
5 1 5 1
6 2 5 2
7 2 9 2
8 2 9 7
9 2 9 7
10 2 9 7
All rows that that are
2010 Aug 11
4
Arbitrary number of covariates in a formula
Hello!
I have something like this:
test1 <- data.frame(intx=c(4,3,1,1,2,2,3),
status=c(1,1,1,0,1,1,0),
x1=c(0,2,1,1,1,0,0),
x2=c(1,1,0,0,2,2,0),
sex=c(0,0,0,0,1,1,1))
and I can easily fit a cox model:
library(survival)
coxph(Surv(intx,status) ~ x1 + x2 + strata(sex),test1)
However, I want to
2011 Sep 28
1
Multiplying a list of matrices with a vector
Hi all!
I have a list of matrices and I want to multiply the ith element of the list with the ith element of a another vector. That is,
> LL <- list(A=diag(3),B=diag(3),C=diag(3))
> vec <- 1:3
> for(i in 1:3)
+ {
+ LL[[i]] <- LL[[i]]*vec[i]
+ }
> LL
$A
[,1] [,2] [,3]
[1,] 1 0 0
[2,] 0 1 0
[3,] 0 0 1
$B
[,1] [,2] [,3]
[1,] 2 0
2010 Aug 31
1
Q about package Icens: change the color of the shading in plot function
Hello!
I want to use the Icens package for analyzing interval-censored data. This code from the manual gives me what I want.
library(Icens)
data(cosmesis)
csub1 <- subset(cosmesis, subset=Trt==0, select=c(L,R))
e1 <- VEM(csub1)
plot(e1)
However, I would like to change the color of the shading from green to something less green, say gray. Any ideas how I could do that? I looked at par, but