Displaying 5 results from an estimated 5 matches for "dtemp".
Did you mean:
temp
2012 Sep 12
1
Plotting US States: Issue with filled colors cycling
...outh Dakota",
"Tennessee","Texas","Utah","Vermont",
"Virginia","Washington","West Virginia","Wisconsin","Wyoming")
dState = data.frame( col=rbinom(51,size=2,prob=.6)+1, state=state )
for( i in 1:51 ){
dTemp = dState[1:i,]
map("state", regions = dTemp$state,
# boundary = TRUE,
# lty = 1, lwd =1,
col= dTemp$col,
fill=T)
readline( paste( "Current i is", i ) )
}
[[alternative HTML version deleted]]
2007 Oct 11
0
Cox with time varying effect
...sed with schoenfeld residuals Chi2 test and graph) so I first want to split time into 2 separate intervals : t<6months and t>=6months, to estimate one hazard ratio (hr) for each interval.
I am analysing Overall survival according to 3 prognostic factors (age,deep,ldh).
my dataframe (dtemp) looks like this:
> dtemp[1:3,]
id time status age_diag deep ldh
1 1 1.544148 1 65.42368 present elevated
2 2 17.051335 1 46.92676 present elevated
4 4 84.829569 0 65.86448 present normal
my first model without time varying hr is:
> mod<-coxph(...
1998 Nov 16
1
PB Mandeville can't be reached
# Peter B. Mandeville kindly offered to send me code for Hotelling's T^2
# Test. Unfortunately there seems to be no route to his machine.
# So i'm trying to reach him via the Mailing List.
------------------------------------------------------------------------
Sir,
this morning i recieved your message about the availability of the code
for Hotelling's Test. I hurried to find out
2010 Feb 11
1
Code find exact distribution for runs test?
...roximation for distribution of the total number
of runs in a sample of n1 observations of one type and n2 observations
of another type.
I've been unable to find R code to generate the exact distribution and
would like to see how this could be done (not homework).
For example, given the data:
dtemp <- c(12, 13, 12, 11, 5, 2, -1, 2, -1, 3, 2, -6, -7, -7, -12, -9,
6, 7, 10, 6, 1, 1, 3, 7, -2, -6, -6, -5, -2, -1)
The Monte Carlo permutation approach seems to get me part way.
# calculate the number of runs in the data vector
nruns <- function(x) {
signs <- sign(x)
runs <- r...
2009 Jun 01
1
survreg.distributions() error
...her function derfun() whose first few lines
look like this:
derfun <- function(y, eta, sigma, density, parms) {
ny <- ncol(y)
status <- y[, ny]
z <- (y[, 1] - eta)/sigma
dmat <- density(z, parms) ## line A: THIS COMMAND DOES NOT
GIVE AN ERROR
dtemp <- dmat[, 3] * dmat[, 4]
if (any(status == 3)) {
z2 <- (y[, 2] - eta)/sigma
dmat2 <- density(z2) ## line B: THIS IS THE COMMAND
THAT GIVES THE ERROR
}
else {
dmat2 <- matrix(0, 1, 5)
z2 <- 0
}
>Fro...