similar to: startup bug (PR#7138)

Hi list, I was trying to use "predict.coxph" to calculate martingale residuals on a test data, however, as pointed out before http://tolstoy.newcastle.edu.au/R/e4/help/08/06/13508.html predict(mycox1, newdata, type="expected") is not implemented yet. Dieter suggested to use 'cph' and 'predict.Design', but from my reading so far, I'm not sure they can

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

For a Monte Carlo study I need to extract from an lme model the estimated standard deviation of a random effect and store it in a vector. If I do a print() or summary() on the model, the number I need is displayed in the Console [it's the 0.1590195 in the output below] >print(fit) >Linear mixed-effects model fit by maximum likelihood > Data: datag2 > Log-likelihood:

Hi If I add a specific machine to my hosts allow list in the global section of my smb.conf, then that machine can access any of the shares on that server. If I only add the machine to a hosts allow list in a specific share, when I try and connect to that share from that machine I get an error : 'Specified computer did not recieve your request' and in log.smb : 'Connection denied from

I'm reposting this problem (perhaps a bug) now I've got more information on it. This is another point of view of the situation and I hope someone could have run into the same trouble before (and solved it :-)) This is it: * with ntbackup 2000 I create a 22Gb .bkf file in the windows machine. * I can copy that file over a samba share and get correct info form the file in windows

Dear all, Has anyone worked on coding for calculating Lyapunov Exponent for a time series data? or any package is available for computing Lyapunov? Please advice and many thanks in advance. Catherine X Wang

Hi Gurus: In the coxph() objects in Survival package, there is an attribute called residuals. Usually, there are several kinds for censored survival data. I can't seem to find in the documentation as to which one this is calculating. Anyone knows? Karen _________________________________________________________________ [[alternative HTML version deleted]]

Is there an implementation of the Cox-Snell residuals / Nelson-Aalen plot for goodness of fit? Or otherwise is there an appropriate Goodness of Fit diagnostic? Thanks Murray -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

R-users -- I would be interested in tools in R to fit the following gam model: logit(p) = a + f(x1) + f(x2) + f(x1,x2), where f(x1,x2) defines a surface. I have looked into the mgcv library, but it seems only to fit models of the form: logit(p) = a + f(x1) + f(x2) Any ideas? Cheers, Dan =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Dan Powers Associate Professor,

Hi Folks, Hopefully this question has enough R and not too much stats to be appropriate for this list. Based on,* Hosmer et al. 2013. Logistic regression for matched case-control studies. Applied Logistic Regression *(eqtn. 7.8)*, *I am assessing GOF of conditional (or matched) logistic regression models with the *standardized Pearson residuals*. The authors define ?large? as delta chi-squared

Hey fellas: In the context of the gambler's ruin problem, the following R code obtains the mean duration of the game, in turns: # total.capital is a constant, an arbitrary positive integer # initial.capital is a constant, an arbitrary positive integer between, and not including # 0 and total.capital # p is the probability of winning 1$ on each turn # 1-p is the probability of loosing 1$ # N

R help list subscribers, I am a new user of R. I am attempting to use R to explore a set of equations specifying the dynamics of a three trophic level food chain. I have put together this code for the function that is to be evaluted by LSODA. My equations Rprime, Cprime, and Pprime are meant to describe the actual equation of the derivative. When I run LSODA, I do not get the output that

Is there any R function that computes the convolution of the double exponential distribution? If not, is there a good way to integrate ((q+x)^n)*exp(-2x) over x from 0 to Inf for any value of q and for any positive integer n? I need to perform the integration within a function with q and n as arguments. The function integrate() is giving me this message: "evaluation of function gave a

Hi all, I know pnorm() is able to calculate a sequence of probabilities by providing a vector of means, and a single number for location and standard deviation. For example, pnorm(0.5,c(0.5,2),1) will produce [1] 0.5000000 0.0668072. However, I wonder if its multivariate counterpart mvtnorm can do the same thing? Namely, if I provide a sequence of (vector-)means and a single vector for

Hello. I am making a program that allows a user to enter a formula to be used in coxph() and then does some other stuff to it, and also allows the user to set other options within coxph(). What is the proper way to pass a formula to a function along with a ...? I have a function: func <- function(coxformula, data = parent.frame(),...) { # Cut to relevant part of function coxmod <-

Hi all, I've uploaded STAR (Spike Train Analysis with R) on CRAN two days ago. The package is designed to analyze neuronal spike (action potential) trains. It uses S3 classes and methods and makes heavy use of other CRAN packages like gss, R2HTML, mgcv, survival. * Analysis of both spontaneous and stimulus evoked activity is implemented for single neuron spike trains as well as for many

Hi all, I've uploaded STAR (Spike Train Analysis with R) on CRAN two days ago. The package is designed to analyze neuronal spike (action potential) trains. It uses S3 classes and methods and makes heavy use of other CRAN packages like gss, R2HTML, mgcv, survival. * Analysis of both spontaneous and stimulus evoked activity is implemented for single neuron spike trains as well as for many

Hi, is there any method for goodness of fit testing of an (as general as possible) univariate distribution with parameters estimated, for normal, exponential, gamma distributions, say (e.g. the corrected p-values for the Kolmogorov-Smirnov or Chi-squared with corresponding ML estimation method)? It seems that neither ks.test nor chisq.test handle estimated parameters. I am aware of function

Hi I want to compare the survival curves in two groups. Because the hazards are not proportional (the curves cross) the log rank test or Cox proportional hazard test are not suitable. How should such curves be compared? Comands are welcome.... Thanks in advance

Hello, My question concerns model selection, stepAIC(), add1(), and coxph(). In Venables and Ripley (3rd Ed) pp389-390 there is an example of using stepAIC() for the automated selection of a coxph model for VA lung cancer data. A statistics question: Can partial likelihoods be interpreted in the same manner as likelihoods with respect to information based criterion and likelihood ratio tests?