search for: lambda0

...situation predict.coxph gives an error message. Namely: stratified data, predict='expected', new data, se=TRUE. I think I found the error but I'll leave that to you to decide. Thanks, Chris ######## CODE library(survival) set.seed(20121221) nn <- 10 # sample size in each group lambda0 <- 0.1 # event rate in group 0 lambda1 <- 0.2 # event rate in group 1 t0 <- 10 # time to estimate expected numbers of events # 'correct' number of events at time t0 = rate of events (lambda) times time (t0) t0*lambda0 t0*lambda1 # generate uncensored data T0 <- rexp(nn, rate=l...

Dear all I want to simulate a stochastic jump variance process where N is Bernoulli with intensity lambda0 + lambda1*Vt. lambda0 is constant and lambda1 can be interpreted as a regression coefficient on the current variance level Vt. J is a scaling factor How can I rewrite this avoiding the loop structure which is very time-consuming for long simulations? for (i in 1:N){ ... N <- rbinom(n=1, size=1...

...nd(rep(1,108),Orthodont$age) #DesignMatrix of RandomEffects Z<-matrix(data=c(rep(1,4),rep(0,108)),nrow=108,ncol=27) #Corr(RanEf)^0.5 = 27 x 27 Identity, since RandomIntercepts are independent sqrt.Sigma<-diag(27) K<-27 #number of subjects/ random intercepts n<-nrow(X) p<-ncol(X) lambda0 <- 2.182382 #actually not a sensible choice as Null-Hypothesis, but that doesn't pertain to the problem #Projection-Matrix for Fixed-Effects-Model: Y -> errors P0=diag(n)-X%*%solve((t(X)%*%X))%*%t(X) mu<-eigen(sqrt.Sigma%*%t(Z)%*%P0%*%Z%*%sqrt.Sigma)$values # mu # [1] 4.00000e+00 4....

Dear all, I am trying to find a suitable R-function for parametric proportional hazard regressions. The package survival contains the coxph() function which performs a Cox regression which leaves the base hazard unspecified, i.e. it is a semi-parametric method. The package Design contains the function pphsm() which is good for parametric proportional hazard regressions when the underlying base

Dear all, 1. How can I compute the intersecting area between 2 polygons ? 2. I have polygons with coordinates in decimal degrees (i.e. 13 deg 30 min = 13.5 decimal degrees). I want to compute their area and get the results in square meters or square kiometers. Can anyone give me a conversion coefficient or a pointer where I can find this information (sorry for this off topic question) ? Thanks

....1,1.1,2.1,3.1,100.1,1.1,2.1,3.1,100.1,1.1,2.1,3.1,100.1,1.1,2.1,3.1,100.1) b = c(1.1,1.1,1.1,100.1,1.1,1.1,1.1,100.1,1.1,1.1,1.,100.1,1.1,1.1,1.1,100.1,1.1,1.1,1.1,100.1,1.1,1.1,1.,100.1) K <-ts(cbind(a,b),start=c(1.1,1.1), names=c("a","b")) fit.bvar <- szbvar(K, p = 1, lambda0=0.6, lambda=0.1, lambda3=2, lambda4=0.25,lambda5=0,mu5=0, mu6=0, prior=0) y <- matrix(c(rep(0,8),c(1.10,1.10,1.10,100.1,1.10,1.10,1.0,100.10)),ncol=2) h <- hc.forecast(fit.bvar,y,nsteps=8,burnin=300, gibbs=500,exog=NULL)

...::DWARFCompileUnit::dump(llvm::raw_ostream&, llvm::DIDumpOptions) /home/linux/llvm-7/llvm/lib/DebugInfo/DWARF/DWARFCompileUnit.cpp:31:0 #22 0x00000000004482ad llvm::DWARFContext::dump(llvm::raw_ostream&, llvm::DIDumpOptions, std::array<llvm::Optional<unsigned long>, 27ul>)::'lambda0'(bool, char const*, llvm::DWARFSection, llvm::iterator_range<std::unique_ptr<llvm::DWARFCompileUnit, std::default_delete<llvm::DWARFCompileUnit> >*>)::operator()(bool, char const*, llvm::DWARFSection, llvm::iterator_range<std::unique_ptr<llvm::DWARFCompileUnit, std::defa...

Dear lmer-ers, My thanks for all of you who are sharing your trials and tribulations publicly. I was hoping to elicit some feedback on my thoughts on denominator degrees of freedom for F ratios in mixed models. These thoughts and practices result from my reading of previous postings by Doug Bates and others. - I start by assuming that the appropriate denominator degrees lies between n

...p;, >> llvm::DIDumpOptions) >> /home/linux/llvm-7/llvm/lib/DebugInfo/DWARF/DWARFCompileUnit.cpp:31:0 >> #22 0x00000000004482ad llvm::DWARFContext::dump(llvm::raw_ostream&, >> llvm::DIDumpOptions, std::array<llvm::Optional<unsigned long>, >> 27ul>)::'lambda0'(bool, char const*, llvm::DWARFSection, >> llvm::iterator_range<std::unique_ptr<llvm::DWARFCompileUnit, >> std::default_delete<llvm::DWARFCompileUnit> >*>)::operator()(bool, char >> const*, llvm::DWARFSection, >> llvm::iterator_range<std::unique_ptr&...

...) >>>> /home/linux/llvm-7/llvm/lib/DebugInfo/DWARF/DWARFCompileUnit.cpp:31:0 >>>> #22 0x00000000004482ad llvm::DWARFContext::dump(llvm::raw_ostream&, >>>> llvm::DIDumpOptions, std::array<llvm::Optional<unsigned long>, >>>> 27ul>)::'lambda0'(bool, char const*, llvm::DWARFSection, >>>> llvm::iterator_range<std::unique_ptr<llvm::DWARFCompileUnit, >>>> std::default_delete<llvm::DWARFCompileUnit> >*>)::operator()(bool, char >>>> const*, llvm::DWARFSection, >>>> llvm::it...

...::TCStrAggregate<NIds::NStr::TCTCStrTraits<NIds::NStr::TCStrTraits<char, NIds::NStr::CDefaultStrParams>, NIds::NStr::TCStrImp_Fixed<NIds::NStr::TCStrTraits<char, NIds::NStr::CDefaultStrParams>, 256ll> > >::f_AddFromIteratorUTF8<NIds::NStr::CStrIteratorUTF8>::'lambda0'(unsigned long)::operator()(unsigned long) const"}, {"_ZZN4NIds4NStr14TCStrAggregateINS0_13TCTCStrTraitsINS0_11TCStrTraitsIcNS0_17CDefaultStrParamsEEENS0_14TCStrImp_FixedIS5_Lx256EEEEEE21f_AddFromIteratorUTF8INS0_16CStrIteratorUTF8EEEvRxRKT_ENKSA_ISB_EUt0_clEm", "void N...