similar to: dropterm in a function

I wanted to add some text to a plot and got (R-1.6.0, Linux): > text(x = c(1, 4), y = 5, labels = x) Error in text.default(x = c(1766, 1895), y = 5, labels = x) : Object "x" not found With the default value of 'labels': > text(x = c(1, 2), y = 5, labels = seq(along = x)) Error in seq(along = x) : Object "x" not found A scoping bug? :) But >

Hi there - I am using R1.5.1 on WinNT and the latest MASS (Venables and Ripley) library. Running the following code: >minimod<-glm(miniSF~gtbt*f.batch+log(mxjd),data=gtbt,family="binomial") >summary(minimod,cor=F) Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.91561 0.32655 2.804 0.005049 ** gtbtgt 0.47171

NAME Douglas Bates EMAIL bates@stat.wisc.edu VERSION 0.63.3 PLATFORM i386-unknown-linux SYSTEM Debian 2.1 CC/FC/MAKE egcs/g77/make NAME Martyn Plummer EMAIL plummer@iarc.fr VERSION 0.63.3 PLATFORM i386-unknown-linux SYSTEM Redhat 5.1 CC/FC/MAKE gcc/egcs-g77/make NAME Göran Broström EMAIL gb@stat.umu.se VERSION 0.63.3 PLATFORM

Consider > oj <- data.frame(x = c(TRUE, FALSE, NA)) > oj x 1 TRUE 2 FALSE 3 NA > summary(oj) x Mode :logical FALSE:1 TRUE :1 But > oj$x <- factor(oj$x) > summary(oj) x FALSE:1 TRUE :1 NA's :1 My point is that NA's should be reported for logicals like they are for other data types. Göran --- Göran

In the first line, use the dist function, found in library mva, to get the distance between each pair of rows. From this calculate an incidence matrix for which element i,j is true if row i in dat equals row j in dat (and false elsewhere). In the second line, for each row calculate the indices of the matching rows and take the minimum of those as the key. incid <-

Hi, I'm working on fitting a glm model to my data using Gamma error structure and reciprocal link. I've been using dropterm (MASS) in the model simplification process, but the F values from analysis of deviance tables reported by dropterm and anova functions are different - sometimes significantly so. However, the reported residual deviances, degrees of freedom, etc. are not different.

Hello, I'm having some difficulty understanding the useage of the "dropterm()" function in the MASS library. What exactly does it do? I'm very new to R, so any pointers would be very helpful. I've read many definitions of what dropterm() does, but none seem to stick in my mind or click with me. I've coded everything fine for an interaction that runs as follows: two sets

SunOS fluke 5.8 Generic_108528-07 sun4u sparc SUNW,Ultra-Enterprise R-1.4.0 Running from the command line R CMD check chlib results in the log file 00check.log: * using log directory `/home/woodstock/hoffmann/R/Sources/chlib.Rcheck' * checking for file `chlib/DESCRIPTION' ... OK * checking package directory ... OK * checking for sufficient/correct file permissions ... WARNING *

I want the indices i for which x[i] < 0 (say): > x <- c(1, -1, 3, 3, -2) > where.negative(x) [1] 2 5 Surely where.negative is something simple, but how? G?ran -------------------------------------------------------------- G?ran Brostr?m Department of Statistics tel: +46 90 786-5223 Ume? University fax: +46 90 786-6614 S-90187 Ume?, Sweden

Dear R fans, I have got a difficult sounding problem. For fitting a linear model using continuous response and then for re-fitting the model after excluding every single variable, the following functions can be used. library(MASS) model = lm(perf ~ syct + mmin + mmax + cach + chmin + chmax, data = cpus) dropterm(model, test = "F") But I am not sure whether any similar functions is

Hello, I am apparently confused about the use of an id parameter for an event history/survival model, and why the EHA documentation for aftreg does not specify one. All assistance and insights are appreciated. Attempting to specifiy an id variable with the documentation example generates an "overlapping intervals" error, so I sorted the original mort dataframe and set subsequent entry

To compare two different models, I've compared the result of using dropterm() on both. Single term deletions Model: growth ~ days + I(days^0.5) Df Sum of Sq RSS AIC <none> 2.8750 -0.2290 days 1 4.8594 7.7344 4.6984 I(days^0.5) 1 0.0234 2.8984 -2.1722 AND Single term deletions Model: growth ~ days + I(days^2) Df Sum

Hi, the link to the R code for repeated measurement analyses of J Lindsey is unfortunately not working. I am desperate for repeated measurements in R; could you please help me out. Sincerely, Dr. G. Stoet -- Dr. Gijsbert Stoet email: stoet@thalamus.wustl.edu Web: http://eye-hand.wustl.edu/lab/people/stoet.html Phone: (314)7474095 Fax: (314)7474370

Only three entries thus far. Please keep them coming in. -k NAME Martyn Plummer EMAIL plummer@iarc.fr VERSION 0.63.3 PLATFORM i386-unknown-linux SYSTEM Redhat 5.1 CC/FC/MAKE gcc/egcs-g77/make NAME Douglas Bates EMAIL bates@stat.wisc.edu VERSION 0.63.3 PLATFORM i386-unknown-linux SYSTEM Debian 2.1 CC/FC/MAKE egcs/g77/make NAME Thomas

Dear List, After including cluster() option the coxreg (from eha package) produces results slightly different than that of coxph (from survival) in the following time-dependent treatment effect calculation (example is used just to make the point). Will appreciate any explaination / comment. cheers, Ehsan ############################ require(survival) require(eha) data(heart) # create weights

On Thu, 18 Apr 2002, Martin Maechler wrote: > >>>>> "Jason" == Jason Turner <jasont@indigoindustrial.co.nz> writes: > > Jason> strptime() mysteriously adds a day to a date, unless the year > Jason> is specified. Tested on: > Jason> Linux (RedHat 6.0) - R version 1.4.1 and R-devel. > Jason> Windows - R version 1.4.1 >

On Thu, 19 Sep 2002 10:43:01 -0400 (EDT), you wrote in message <200209191443.KAA17404@falmouth.bwh.harvard.edu>: >Hi Duncan -- I am curious as to whether you are building/distributing >R 1.6.0 for windows? There's nothing urgent about it, I have one. >The question is how to identify a distribution URL if one exists. I will be building and distributing the final release.

Hi, I have sent this message to this list the July, 7th. It was about a problem in MASS package. Until now there is no change in the devel version. As the problem occurs in a package and not in the R-core, I don't know if the message should have been sent here. Anyway, I have added a copy to Pr Ripley. I hope it could have been fixed. Sincerely Marc Le 09/07/2017 ? 16:05, Marc Girondot via

Hello, I'm trying to use coxph() function to fit a very simple Cox proportional hazards regression model (only one covariate) but the parameter space is restricted to an open set (0, 1). Can I still obtain a valid estimate by using coxph function in this scenario? If yes, how? Any suggestion would be greatly appreciated. Thanks!!! [[alternative HTML version deleted]]

I have a problem with some functions handling dates, in packages 'date' and 'survival' (they seem to be identical). For instance, from the documentation, -------------------- mdy.date {survival} R Documentation Convert to Julian Dates Description Given a month, day, and year, returns the number of days since January 1, 1960. Usage mdy.date(month, day, year, nineteen = TRUE,