search for: jasa

...library(survival) ## Example from help page of survdiff ## Expected survival for heart transplant patients based on ## US mortality tables expect <- survexp(futime ~ ratetable(age=(accept.dt - birth.dt), sex=1,year=accept.dt,race="white"), jasa, cohort=FALSE, ratetable=survexp.usr) Error in floor(temp) : Non-numeric argument to mathematical function sessionInfo('survival') R version 2.12.1 Patched (2010-12-18 r53869) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=German_Switzerland.1252 LC_CTYPE=Ger...

I have linux box (does nat and firewall for small network) connected to dsl. I want to set priorities for protocols (that nothing could disturb web browsing). This is my rules (eth0 connected to internet): /sbin/tc qdisc del dev eth0 root /sbin/tc qdisc add dev eth0 root handle 1 htb default 30 r2q 100 /sbin/tc class add dev eth0 parent 1: classid 1:2 htb rate 900Kbit burst 15k /sbin/tc class

A user asked me about this and I can't figure it out. tmt% R R Under development (unstable) (2018-04-09 r74565) -- "Unsuffered Consequences" Copyright (C) 2018 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) > library(survival) > data(cgd0) Warning message: In data(cgd0) : data set ?cgd0? not found ---- The data set is present and can be

...Product: swfdec Version: unspecified Platform: x86 (IA32) OS/Version: Linux (All) Status: NEW Severity: normal Priority: medium Component: library AssignedTo: swfdec at lists.freedesktop.org ReportedBy: jasa.david at gmail.com QAContact: swfdec at lists.freedesktop.org There's changing text in file I'll attach. Problem is, that when "new" line appears, the old one doesn't go away and characters "stack" on each other. -- Configure bugmail: http://bugs.freede...

Dear R community, I have two questions regarding fitting GLMM using maximum likelihood method. The first one arises from trying repeat an analysis in the Breslow and Clayton 1993 JASA paper. Model 3 of the epileptic dataset has two random effects, one subject specific, and one observation specific. Thus if we count random effects, there are more parameters than observations. When I try to run the following code, I get an error saying: "Error in mer_finalize(ans) : q = 295 &...

Hi, Storer and Crowley (JASA 1985) presented an approach for approximating the changes in maximum partial-likelihood parameter estimates for the Cox model when a single observation is deleted. Is there an R implementation of this approach? Any help is greatly appreciated. Thanks. Best, Ravi. ---------------------------...

...Loading required package: splines > ## Example from help page of survdiff > ## Expected survival for heart transplant patients based on > ## US mortality tables > expect <- survexp(futime ~ ratetable(age=(accept.dt - birth.dt), + sex=1,year=accept.dt,race="white"), jasa, cohort=FALSE, + ratetable=survexp.usr) Error in floor(temp) : Non-numeric argument to mathematical function > sessionInfo('survival') R version 2.12.1 Patched (2010-12-18 r53869) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=German_Switzerland.1252 LC_CTYPE=Germ...

We have put a new package sparsenet on CRAN. Sparsenet fits regularization paths for sparse model selection via coordinate descent, using a penalized least-squares framework and a non-convex penalty. The package is based on our JASA paper Rahul Mazumder, Jerome Friedman and Trevor Hastie: SparseNet : Coordinate Descent with Non-Convex Penalties. (JASA 2011) http://www.stanford.edu/~hastie/Papers/Sparsenet/jasa_MFH_final.pdf We use Zhang's MC+ penalty to impose sparsity in model selection. This penalty parametrizes a fam...

We have put a new package sparsenet on CRAN. Sparsenet fits regularization paths for sparse model selection via coordinate descent, using a penalized least-squares framework and a non-convex penalty. The package is based on our JASA paper Rahul Mazumder, Jerome Friedman and Trevor Hastie: SparseNet : Coordinate Descent with Non-Convex Penalties. (JASA 2011) http://www.stanford.edu/~hastie/Papers/Sparsenet/jasa_MFH_final.pdf We use Zhang's MC+ penalty to impose sparsity in model selection. This penalty parametrizes a fam...

...e x-min and x-max of a data set?) Rather than speculate what causes it, I am including a short R snippet and a short data set that demonstrates it. is this a bug or a feature? * it would be nice if lowess was a little better documented. I have easy access to Becker-Chambers-Wilks, but not to JASA and AS. I wish "?lowess" would tell me a little more about the method. thanks in advance. regards, /iaw -------- test.R test <- read.table(file="test.txt", sep="\t", header=T); plot( test$x, test$y, log="xy"); lines( lowess( test$x, test$y, f=0.5 ) );...

Is it appropriate to use stepAIC (library MASS) with coxph objects (from library survival5) that use "cluster(id)"? It is my understanding that, when using "cluster(id)", we can test for sets of terms by using the methods in Wei et al., (1989; JASA, 84: 1065-1073), or as explained in pp. 53 and ff. of the survival.ps document. But if we use a likelihood ratio test instead (comparing sets of nested models), we can get different results. I started wondering, then, if using AIC, per se, could give misleading results. Am I getting confused here?...

hello everyone, i'm trying to simulate a 2-parameter integrated ornstein-uhlenbeck (IOU) process, but i'm not sure exactly where to start (which package, which function). the motivation is the paper by taylor et. al. (JASA 1994) "a stochastic model for the analysis of longitudinal aids data." the model they suggest consists of a combination of fixed and random effects, a stochastic process, and measurement error. the 2-parameter IOU process they suggest has cov[W(s), W(t)] = (sigma^2/(2alpha^3)[2alpha min...

...x) Now I wanted to make a plot of survival probability vs. the covariate, and the 95% confidence interval for the survival probability. It's just like a Kaplan-Meier Survival curve, except now the x axis represents the value of covariate, not the time. Someone gave me a reference to a paper in JASA by Gary (1992) for this type of plot, but I didn't have the access to the paper. So I am wondering if anyone knows how to do this in R or S-Plus? In addition, can anyone explain to me what are the following "type" options in predict.coxph() predicting? predict(fit,type='lp',...

...ts, that is, not only trying the simple abrupt permanent impact (yt = w It ) with the xreg option but also trying with a gradual permanent impact (yt= d * yt-1 + w * It ), following the filosophy of Box and Tiao (1975). Intervention analysis with applications to economic and environmental problems. JASA 70: 70-92. Does anybody know where could I find how to incorporate them using the arima comand (or other), or a statistical package which can incorporate it? Thanks, Berta. [[alternative HTML version deleted]]

...treatment effect by a surrogate in a survival model (Lin, Fleming, and De Gruttola 1997 in Statistics in Medicine). The paper mentioned that the covariance matrix matches that of the covariance matrix estimator for the marginal hazard modelling of multiple events data (Wei, Lin, and Weissfeld 1989 JASA), and is implemented in Lin's MULCOX2, SAS, and S-plus. Is this the way to fit such a model in R? Suppose I have variables: time, delta, treatment, and surrogate. Should I repeat the dataset (2x) and stack, creating the variables: time1 (time repeated 2x), delta1 (delta repeated 2x), treatmen...

Hi, I'm a newbie to audio compression. And I have several question about ogg. 1. How many bits a vorbis sample take ? does the sample size 16 bits ? like the mp3 ? CMIIW 2. How do a decoder recognize how many sample it takes in a second ? does it related with the code books ? 4. Does vorbis use VBR ? 3. What is function of the Windows ? and what is the output of the windows ? Thanks

...survival_Rdata $bladder [1] "bladder" "bladder1" "bladder2" $cancer [1] "cancer" $cgd [1] "cgd" "cgd0" $colon [1] "colon" $flchain [1] "flchain" $genfan [1] "genfan" $heart [1] "heart" "jasa" "jasa1" ------- snip ---- data(bladder2) # Warning message: # In data(bladder2) : data set ?bladder2? not found data(mgus1) # Warning message: # In data(mgus1) : data set ?mgus1? not found > data(mgus1) Warning message: In data(mgus1) : data set ?mgus1? not found > data...

...lo, My question likely got buried so I am reposting it in the hopes that someone has an answer. I have thought more about the question and modified my question. I hope tha my specific question is: I am attempting to create a bootstrap procedure for a finite sample using the theory of Rao and Wu, JASA (1988) that replicates within each strata (h) n_h - 1 times. To this end, I require a different sample size for each stratum. Unfortunately, it appears that the sample command which is used to obtain my resamples only allows a scalar for the sample size. i.e it does not currently work with a vecto...

Barrett & Ling, JASA, 1992, v.87(417), pp184-191 define general classes of influence measures for multivariate regression models, including analogs of Cook's D, Andrews & Pregibon COVRATIO, etc. As in univariate response models, these are based on leverage and residuals based on omitting one (or more) obser...

Dear All, I'm a beginner user in R and I would like to make a quadratic and plateau model in R. Can you help please with an example? Thanks so much -- Ahmed M. Attia Assistant Lecturer El-Khattara farm Station Agronomy Dept., Zgazig Univ., Egypt Visiting Scientist Haskell Agricultural laboratory Agronomy and Horticultural Dept., Univ. of Nebraska-Lincoln ahmedatia at zu.edu.eg