similar to: Simple question, I think

Hi all, In using the identify command, I get the following message > plot(hatvalues(scireg3)) > abline(h=.0154,lty=2) # plots a reference line at (k + 1)/n > identify(1:1165, hatvalues(scireg3),row.names(sciach)) Error in xy.coords(x, y) : 'x' and 'y' lengths differ which doesn't allow me to see the observation number when I scroll over with the mouse. What

Let me try to be more specific. The x y coordinates are different because of NAs in the dataset. In this analysis, a set of hat values (a measure of influence in regression) is given for each observation. On the basis of the regression that was run to get these hat values, the sample size was 1164 (one removed due to NA). The length of the data set is 1165. If I remove the NA from the

Hi all, I am trying to draw a Kaplan-Meier curve and I found online that Kaplan - Meier estimates are computed with a function called km in the event package. Is there an update for that because when I choose to download packages in R,. there is no package called event, even though I have selected all the repositories. Thanks in advance, Eleni [[alternative HTML version deleted]]

Hi, I am attempting to graph a Kaplan Meier estimate for some claims using the survfit function. However, I was wondering if it is possible to plot a cdf of the kaplan meier rather than the survival function. Here is some of my code: library(survival) Surv(claimj,censorj==0) survfit(Surv(claimj,censorj==0)~1) surv.all<-survfit(Surv(claimj,censorj==0)~1) summary(surv.all) plot(surv.all)

[This email is either empty or too large to be displayed at this time]

Hi all, My apologies for a very simple question. I just downloaded R 2.6.0. I want to bring in all of the objects from 2.5.0 that I see when I type ls(). I have no idea how to do that. Thanks in advance. David -- ======================================================================= David Kaplan, Ph.D. Professor Department of Educational Psychology University of Wisconsin - Madison

Hi all, I'm looking for a contributed package that can provide a detailed account of missing data patterns and perhaps also provide imputation procedures, such as mean imputation or hot deck imputation and the like. Is there anything out there? Thanks in advance, David -- =========================================================================== David Kaplan, Ph.D. Professor

Greetings, How to I read in SPSS .sav files into R. Thank you. David ===================================== David Kaplan, PhD Professor of Education School of Education University of Delaware Newark DE 19716 Voice: 302-831-8696 Fax: 302-831-4110 email: mailto:dkaplan at udel.edu Homepage:

Greetings, I'm interested in generating data from various bivariate or mulitivariate distributions (e.g. gamma, t, etc), where I can specify the parameter values, including the correlations among the variables. I haven't been able to dig anything up on the faq, but I probably missed something. A nudge in the right direction would be appreciated. David --

Hi all, I'm looking for R packages that estimate multivariate time-series models or vector-autoregression (VAR) time-series models. Thanks David -- =========================================================================== David Kaplan, Ph.D. Professor Department of Educational Psychology University of Wisconsin - Madison Educational Sciences, Room, 1061 1025 W. Johnson Street Madison,

Hi R experts, I have a data recoding problem I cant get my head around - I am not that great at the subsetting syntax. I have a dataset of longitudinal toxicity data (for multistate modelling) for which I want to also want to do a simple Kaplan-Meier curve of the time to first toxic event. The data for 2 cases presently looks like this (one with an event, the other without), with id representing

I can't seem to figure out why I'm getting this error. The output is copied right off the screen. Notice how in some cases the back slash is missing. In other cases, it can't read a file that I know is there. Thanks in advance > library(foreign) > hrout <- read.spss("c:\\hrab200.sav") Error in read.spss("c:\\hrab200.sav") : unable to open file

Hi When I use the read.table function with header = T, I notice that it gives me the variable names along the top as I expect. But, when I then attempt an analysis, e.g. regression, it doesn't recognize the variable names. Am I missing a step. Thank you David [[alternative HTML version deleted]]

Hi, Is there a package that will estimate simple Markov models and hidden Markov models for discrete time processes in R? Thanks in advance, David -- =============================================================== David Kaplan, Ph.D. Professor Department of Educational Psychology University of Wisconsin - Madison Educational Sciences, Room, 1061 1025 W. Johnson Street Madison, WI 53706

Hi all, I have started using LaTeX for writing papers and I have heard that R works well with LaTeX. I'm specifically interested in how I can have LaTeX read in R generated graphics - for example graphs formed by matplot, or other such processes. Does anyone out there use LaTeX and can point me in the right direction? Thanks David --

Hi, I'm plotting 5 autocorrelation plots on one page. Using par(mfrow=c(3,2)) everything comes out fine. However, for each plot, it prints a title on top of each plot that says Series followed by the variable name used in the plot. I want to suppress those titles, but I also want a general figure title on the bottom of the page. I've looked at the Murrell book as well as the acf

Greetings all, I have two questions. I have a data set that is arranged in the example below. I wish to obtain a plot of the performance of each ID over Year on v1. It's not clear how I set this up? ID Year V1 1 1980 1 1 1981 2 1 1982 6 1 1983 4 2 1980 5 2 1981 5 2 1982 5 2 1983 6 Also,I would like to transpose the data to have the

Greetings all, I have two logistic plots coming from two calls to plogis. The code is .x <- seq(-7.6, 7.6, length=100) plot(.x, plogis(.x, location=0, scale=1), xlab="x", ylab="Density", main="Logistic Distribution: location = 0, scale = 1", type="l") abline(h=0, col="gray") .y <- seq(-7.6, 7.6, length=100) plot(.x, plogis(.x,

Greetings, I'm trying to simply reorder a data frame on the row numbers. So, for example, instead of getting 1,2,3,4,5,6,7,8,9,10,11, ... 100 ..., I get instead 1, 10, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 11, ... I've tried commands such as df <- df[order(rownames(df)),] and and have substituted the order command with sort and sort.list to no avail. Any advice

Greetings all, I have three files that I would like to merge. The first is a student level file that contains the student ID, their school ID, and their country ID. The second is the school file that contains the school ID and country ID. The third is the country file with a country ID. I would like to merge the three together using the common country ID. Of course, what this would mean