similar to: Enquiry on Vrtest

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

Full_Name: Rob Strawderman Version: 1.9.1 OS: red hat linux 8 Submission from: (NULL) (128.253.217.211) The following demonstrates the basic problem: [martingale <99>]% rm .RData rm: remove regular file `.RData'? y [martingale <100>]% R R : Copyright 2004, The R Foundation for Statistical Computing Version 1.9.1 (2004-06-21), ISBN 3-900051-00-3 [deleted] >

Dear R-helpers, I am using a vrtest on time series data. My commands are as follows; read.table("B.txt",sep="\t",fill=TRUE, na.strings = "NA") require(vrtest) rm(list=ls(all=TRUE)) datamat <- read.table("B.txt",sep="\t",fill=TRUE, na.strings = "NA") column <- 1 nob <- nrow(datamat) y <-

Hello I am a new user of R software. I benefit from using vrtest-package. However, the codes provided by the aforementioned package, for example, calculate the test statistics for Lo and Mackinlay (1988) under the assumptions of homoscedasticity and heteroscedasticity without computing the value of the variance ratio itself.   I would be grateful if you could instruct me how to calculate the

Hello list. I am hoping for some help with a relatively simple problem. I have a data frame arranged as below. I want to be able to count the occurrence of each gene (eg let-7e) by Experiment. In other words how many times does a given gene crop up in the dataframe. I tried table but couldn't work out how to get the output I want. I have also considered rearranging this data into a list (by

Hi, How can I load the alpm(4) module for SMBus support on my ASUS Vintage AH-1 system? It appears the I/O range it uses is claimed by the acpi(4) driver. Can I override the mapping in some way i.e. tell ACPI not to claim the range? I don't see anything obvious about this in the acpi(4) man page. thanks, BMS

I apologize about doubling these up, I forgot the subject! I have a cisco VPN from router to router over a Data T-1. The ping times are consistently 32ms with random ping responses of 295ms -408ms about every 30 secs to a minute, I have jitter buffer enabled. The connection goes like this... Mitel SIP phone to Asterisk A, IAX trunked to Asterisk B and then T-1 to a PBX, the calls are internal

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]]

I guess my CentOS is screwed up because that function doesn't work. Only way to resize it was to go into the menu areas and set it there. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.centos.org/pipermail/centos/attachments/20060412/4f37d0af/attachment-0001.html>

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

I just downloaded samba-1.9.18p3.tar.gz and compiled it on a DEC 3000/300 running OSF/1 v3.0 (uname -mrsv gives: "OSF1 V3.0 347 alpha") using the built-in cc compiler. The compiler gives the following warning: Compiling nmbd_incomingdgrams.c /usr/lib/cmplrs/cc/cfe: Warning: nmbd_incomingdgrams.c, line 452: Unknown character \ ignored ((DEBUGLEVEL>=(5))?(Debug1

Thanks for the help guys. To resize my screen I had to go through the GUI part and change it but unlike Microsoft it makes me reboot in order for it to take affect which sucks. The Ctrl Alt +/- does not work. As for adding a printer someone finally gave me the link to the correct document which showed how to do it. I received a link where you can download the CentOS manual in PDF format:

Dear Rhelp I am first time user of R and I have the following questions: 1. How I can load my data (it is a matrix 10x1618) with the titles of the variables in the first row? 2. How I can run vrtest to calculate Variance Ratio Tests for Weak-form Market Efficiency? Thanks in advance. Jasim ____________________________________________________________________________________ the

Hello, I am looking for a code in R for the variance ratio test statistic (the Lo and Mackinlay version or any other versions). Does anybody have such a code they can share or know a library in which I can find this function? Basically I have a number of time series which I need to check for persistence. One other test I can use is the runs test in the tseries package. Any help will be greatly

WinXP-x32, R-21.13.0 Dear list, I have a problem that (I think) relates to the interaction between Windows and R. I am trying to scrape a table with data on the Hawai'ian Islands, This is my code: library(XML) u <- "http://en.wikipedia.org/wiki/Hawaii" tables <- readHTMLTable(u) Islands <- tables[[5]] The output is (first set of columns):

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 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,