Displaying 10 results from an estimated 10 matches for "kappa1".
Did you mean:
kappa
2005 Feb 01
0
RV: problems checking a package
...ave a very strange problem. I made a package (under Windows and
Linux). The package passed the R CMD Check without problem. Then, I
installed the package and executed a function which calls to a 'dll'
mod<-frailtyPenal(Surv(time,status)~sex+age+cluster(id),
+ n.knots=8,kappa1=10000,data=kidney)
mod
Call:
frailtyPenal(formula = Surv(time, status) ~ sex + age + cluster(id),
data = kidney, n.knots = 8, kappa1 = 10000)
Shared Gamma Frailty model parameter estimates
using a Penalized Likelihood on the hazard function
coef exp(coef) se(coef) se(coef) HI...
2013 Mar 13
1
saving vector output as numeric
...al data to R.
An example of my original data frame minmaxfunc is as follows:
min max
T1 1.500000e+01 3.999954e+01
SE1 0.000000e+00 1.000000e+00
PRE 0.000000e+00 1.000000e+00
WET 0.000000e+00 5.980000e+00
BE1 4.664642e+00 5.866620e+00
Kappa1 5.500000e+03 2.000000e+04
Kappa3 1.000000e+04 2.000000e+04
Then I created a latin hypercube set using (qunif(x[,i],
minmaxfunc$min[i], minmaxfunc$max[i]). The new data frame looks as
follows:
T1 SE1 PRE
WET BE1 Kappa1 Kap...
2009 Nov 10
0
NEW release of FRAILTYPACK
...events.
For instance :
-- ADDITIVE FRAILTY MODELS for proportional hazards models with two
correlated random effects (intercept random effect with random slope):
fitAdditive<-additivePenal(Surv(t1,t2,event)~cluster(group)+var1+slope(var1),
correlation=TRUE,data=dataAdditive,n.knots=8,kappa1=10000)
-- NESTED FRAILTY MODELS for hierarchically clustered data (with 2
levels of clustering) by including two random effects:
fitnested<-frailtyPenal(Surv(t1,t2,event)~cluster(group)+subcluster(subgroup),
data=dataNested,n.knots=8,kappa1=50000)
-- JOINT FRAILTY MODELS in the cont...
2009 Nov 10
0
NEW release of FRAILTYPACK
...events.
For instance :
-- ADDITIVE FRAILTY MODELS for proportional hazards models with two
correlated random effects (intercept random effect with random slope):
fitAdditive<-additivePenal(Surv(t1,t2,event)~cluster(group)+var1+slope(var1),
correlation=TRUE,data=dataAdditive,n.knots=8,kappa1=10000)
-- NESTED FRAILTY MODELS for hierarchically clustered data (with 2
levels of clustering) by including two random effects:
fitnested<-frailtyPenal(Surv(t1,t2,event)~cluster(group)+subcluster(subgroup),
data=dataNested,n.knots=8,kappa1=50000)
-- JOINT FRAILTY MODELS in the cont...
2010 Aug 19
1
memory problem
...how can i solve it?
thanks in advance
Avi
my code
# frailtypack
library(frailtypack)
cgd.ag <- read.csv("C:/rfiles/RE/cgd.csv")
cgd.nfm <-frailtyPenal(Surv(TStart, TStop,
Status)~cluster(Center)+subcluster(ID)
Treatment,data=cgd.ag,Frailty=TRUE,n.knots=8,kappa1=50000,
cross.validation=TRUE,recurrentAG=TRUE)
--
View this message in context: http://r.789695.n4.nabble.com/memory-problem-tp2330510p2330510.html
Sent from the R help mailing list archive at Nabble.com.
2013 Mar 14
0
Error with epiR and sensitivity
...oblem with the data itself rather than
the syntax.
An original data frame of
min max
T1 1.500000e+01 3.999954e+01
SE1 0.000000e+00 1.000000e+00
PRE 0.000000e+00 1.000000e+00
WET 0.000000e+00 5.980000e+00
BE1 4.664642e+00 5.866620e+00
Kappa1 5.500000e+03 2.000000e+04
Kappa3 1.000000e+04 2.000000e+04
produces a latin hypercube sampling set called samplevalues:
T1 SE1 PRE
WET BE1 Kappa1 Kappa3
1 31.35590 0.7066388715 0.8665111432 4.965701530 5.783424 12240...
2009 Jan 26
1
Error in Surv(time, status) : Time variable is not numeric
...lt using SAS software. I also tried to change the format
using SPSS and Excell.
My (reduced) dataset has following column names:
ID entry time status family var1
I used following command:
> frailtyPenal(Surv(time, status) ~var1 + cluster(family), Frailty=TRUE
> ,n.knots=8, kappa1=1500,
+ cross.validation=FALSE)
And got this error :
Error in Surv(time, status) : Time variable is not numeric
In addition: Warning message:
In is.na(time) : is.na() applied to non-(list or vector) of type 'closure'
I think R transforms the data when importing into R, so that the
observa...
2004 Feb 03
1
Error in f(x, ...) : subscript out of bounds
...-wired in (not data dependent and not parameter dependent),
furthermore, the constant is not used in any subscripting.
Sorry I cannot provide a toy example, but the likelihood function looks
like this:
lnL <- function(theta, gsvr, gsvR) {
# theta[1]= mu, theta[2]=gamma, theta[3]=kappa1,
theta[4]=kappa2
nn <- 252
d <- 0:nn
wd <- exp((theta[3] * d + theta[4] * d^2))/(sum(exp(theta[3] *
d + theta[4] * d^2)))
sigsq <- numeric(length(gsvR$ret))
x <- numeric(length(gsvR$ret)
# Below this line can be s...
2010 Aug 09
2
recurrent events
Hello,
I have a cohort with approx 1,200 patients at the ages of 30-65 that had
their first myocardial infarction during 1992:
· They were in a follow up until 2005.
· About 400 of them died during this period of time (right censored)
· Each one of them had up to 4 mi recurrent events.
I am using the semi-parametric model in order to assess the relationship of
2006 May 12
3
Maximum likelihood estimate of bivariate vonmises-weibulldistribution
Thanks Dimitris!!! That's much clearer now. Still have a lot of work to
do this weekend to understand every bit but your code will prove very
useful.
Cheers,
Aziz
-----Original Message-----
From: Dimitrios Rizopoulos [mailto:Dimitris.Rizopoulos at med.kuleuven.be]
Sent: May 12, 2006 4:35 PM
To: Chaouch, Aziz
Subject: RE: [R] Maximum likelihood estimate of bivariate