similar to: Enquiry about 2nd-order interactions survival analysis

Hi experts, http://r.789695.n4.nabble.com/file/n4034318/Parametric_survival_analysis_2nd-order_efffect.JPG Parametric_survival_analysis_2nd-order_efffect.JPG As we know a normal survival regression is the equation (1) Well, I'ld like to modify it to be 2nd-order interaction model as shown in equation(2) Question: Assume a and z is two covariates. x = dummy variable (1 or 0) z = factors

Hello, I am new to R on Linux (and to LINUX), running a 1.9.1 on a Linux MDK10.0. When updating or installing different packages, I get messages like this (under), the package doesn't install. Can anyone help me? I know this is all to little information, but I don't know were to start. I think there is something wrong in how or were the C++ or fortran compilator is installed?? Thanks, Geir

Thanks to Roger Bivand, my problem was solved. binutils 2.14 was updatet to binutils 2.15, as indicated under. Geir Helge Systad Norwegian Institute for Nature Research Arctic Ecology, Polarmilj??senteret, N-9296 Troms??, Norway Adress aug-04 to feb-05: Albert Pettersonsvei 13, N-5750 Odda fax +47 85 03 82 14 phone +47 53 64 24 94 mobil +47 91 63 70 55 E-mail: geir.systad at nina.no

I thought I'd give this another shot before the official release. I tried building R-beta_2004-04-07 on the AIX system that I have access to, and it seemed to failed at lazy-loading survival. I'd very much appreciate any pointer on what to try or look for next. 1. I set OBJECT_MODE to 64 for building 64-bit binary. 2. I edited config.site with the following: CC="xlc_r"

Full_Name: D Kreil Version: 2.1.1 OS: HP-UX B.11.23 U ia64 0029870451 unlimited-user license Submission from: (NULL) (62.178.15.60) I sadly get a core-dump during installation, when it tries to lazyload package survival. Any idea what to do about this? Best regards, D. ... include -fPIC -g -O2 -c survindex2.c -o survindex2.o gcc -I/bi/common/src/R-2.1.1/include

ryusuke: It sounds like you need to have (D)COM server to be able to work on the background. As for the foreground, rcom is what you need and it appears to be working OK. I am out of my office right now but I'll back to work next week and I will be able to explain in more detail how the applications I tied together work with RExcel and R. You can download the latest staconnDCOM.latest.exe from

Hello, The likelihood includes two parameters to be estimated: lambda (=beta0+beta1*x) and alpha. The algorithm for the estimation is as following: 1) with alpha=0, estimate lambda (estimate beta0 and beta1 via GLM) 2) with lambda, estimate alpha via ML estimation 3) with updataed alpha, replicate 1) and 2) until alpha is converged to a value I coded 1) and 2) (it works), but faced some

Hi, I'm trying to get maximum likelihood estimates of \alpha, \beta_0 and \beta_1, this can be achieved by solving the following three equations: n / \alpha + \sum\limits_{i=1}^{n} ln(\psihat(i)) - \sum\limits_{i=1}^{n} ( ln(x_i + \psihat(i)) ) = 0 \alpha \sum\limits_{i=1}^{n} 1/(psihat(i)) - (\alpha+1) \sum\limits_{i=1}^{n} ( 1 / (x_i + \psihat(i)) ) = 0 \alpha \sum\limits_{i=1}^{n} (

Dear All I'm attempting to speed up my model building procedure, but need some help with the loops I've created...please bear with me through the explanation! My basic model call is something like: m0sulf.nlme<-nlme(conc~beta0*exp(-beta1*day)+beta2*exp(-beta3*day), data=m0sulf, fixed=(beta0+beta1+beta2+beta3~1),

Hi all, I am trying to create simulated data for exploring reclassfication measures in a logistic setting with two continuous predictors and I would like to set the average population probability of outcome rather than the logistic beta0. Is there a way to find a beta0 that will generate the desired average population probability of outcome given x,y and their odds ratios? #Here is an

I want to compute a new variable (newvar) based on the values of two other variables (t1freq, t2freq).The two variables (t1freq,t2freq) are contained in a dataframe - study1dat <- read.csv("c:\\study1rb.csv",header=T) . I gather this computation can be done using Apply and I have run the following example from the help menu ## Compute row and column sums for a matrix: x

Dear R-help, I am trying to estimate a Cox model with nested effects basing on the minimization of the overall AIC; I have two frailties terms, both gamma distributed. There is a error message (theta2 argument misses) and I don?t understand why. I would like to know what I have wrong. Thank you very much for your time. fitM7 <- coxph(Surv(lifespan,censured) ~ south + frailty(id,

I use "while" loop but it produces an errro. I have no idea about this. Error in "[<-"(`*tmp*`, i, value = numeric(0)) : nothing to replace with The problem description is The likelihood includes two parameters to be estimated: lambda (=beta0+beta1*x) and alpha. The algorithm for the estimation is as following: 1) with alpha=0, estimate lambda (estimate beta0

Hi all-- We are planning an intervention study for adolescent alcohol use, and I am planning to use simulations based on a hurdle model (using the hurdle() function in package pscl) for sample size estimation. The simulation code and power code are below -- note that at the moment the "power" code is just returning the coefficients, as something isn't working quite right. The

Hi there, I am currently working on my dissertation which is about testing the martingale hypothesis in the stock market using a methodology involving a range of variance ratio tests and multiple variance ratio tests. I contacted the author of a reference paper and I was told that the tests can be conducted using R programming language. Although I have gone through the theoretical background of

I have downloaded samba-2.2.2-sparc-solaris-2.8.pkg from your site and would like to install it to my unix server sun solaris 5.8. I am battling to use the pkgadd to install this package. Can you please help me in this regard as I have failed to find samba/docs/htmldocs/UNIX_INSTALL which is suppose to be instructions on how to install it on unix sun solaris. Thanking you in advance. Regards,

This is a minor question, but is there any difference between the two objects fred <- structure(1:10, dim = c(2,5)) and fred <- structure(1:10, .Dim = c(2,5)) Should I be using one rather than the other? Thanks, Jonathan. Jonathan Rougier Science Laboratories Department of Mathematical Sciences South Road University of Durham Durham DH1 3LE

Reid Spencer wrote: >Yes, that's right! > >In fact, shortly the process of doing that will get easier with the >llvmc (compiler driver) tool that I'm working on. You write your >compiler to generate either bytecode or LLVM assembly and a > > Does that mean the front end must decide between emiting bytecode for interpretting/JITing and LLVM assembly for native

On Thu, Aug 19, 2004 at 10:35:02AM +1200, Peter Ashford wrote: > Does that mean the front end must decide between emiting bytecode for > interpretting/JITing and LLVM assembly for native compilation? You > can't emit the one kind of output for either end target (interpretted > bytecode or native compilation)? Not at all, there's only one version of LLVM IR, but it can

Misha Brukman wrote: >On Thu, Aug 19, 2004 at 10:35:02AM +1200, Peter Ashford wrote: > > >>Does that mean the front end must decide between emiting bytecode for >>interpretting/JITing and LLVM assembly for native compilation? You >>can't emit the one kind of output for either end target (interpretted >>bytecode or native compilation)? >>