similar to: Matrix of dummy variables from a factor

Does anyone know what could be causing the following error message on startup of R: Error in .addMethodFrom(def,argName[1],class[1],fromClass) : object "*tmp*" not found. I'm using R 1.7.0 and Windows 2000. This happened shortly after I installed 1.7.0. Once it occurs, R freezes. >From then on R will always freeze with this error message. I reinstalled 1.7.0 and it worked

xyplot doesn't seem to want to label my x-axis with dates but instead puts the day-number for each date. begdate is the number of days since January 1, 1960 and was initially created by library(date) ... polls$begdate<-mdy.date(begmm,begdd,begyy) I create a new dataframe (pollstack) which includes begdate. In the process begdate seems to lose its date attribute so I redo it as: >

A newbie question: I'm trying to decide whether to run a maximum likelihood estimation in R or Stata and am wondering if the R mle routine is reasonably robust. I'm fairly certain that, with this data, in Stata I would get a lot of complaints about non-concave functions and unproductive steps attempted, but would eventually have a successful ML estimate. I believe that, with the

Dear list In HLM, one can specify a covariate at one of the "levels". For example, if the data structure are repeated observations nested within students nested within schools, school size might be a covariate that is used at level 3, but not at the other levels. In HLM this is rather easy to do. However, how can one specify a covariate in R for only one of the levels? I have a

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

Is there any way to echo comments from an R source file into an SWeave->LaTeX document? Example: # Npop is population total # Npoph0..Npoph2 are stratum totals # Npoph is vector of stratum totals Npop<-sum(to2000) Npoph0<-sum(to2000[bg==0]) Npoph1<-sum(to2000[bg==1]) Npoph2<-sum(to2000[bg==2]) Npoph<-c(Npoph0,Npoph1,Npoph2) In the final LaTeX document, I'd like the

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

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

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

Hi, df <- read.table("data.txt", header=T); library(nls); fm <- nls(y ~ a*(x+d)^(-b), df, start=list(a=max(df->y,na.rm=T)/2,b=1,d=0)); I was using the following routine which was giving Singular Gradient, Error in numericDeriv(form[[3]], names(ind), env) : Missing value or an Infinity produced when evaluating the model errors. I also tried 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 have a very simple question that I have not been able to solve by reading the manual. When I produce a stripchart with two straight columns of dots representing individual observations, one representing one group of subjects and the other representing another, the columns wind up at the far left and far right sides of the plot, and I'd like them to be closer together, about as far

I am embarrassed to have to ask this but can anyone tell me of a Tukey multiple comparisons procedure available for R? I have looked through the search page, through the FAQ, and in the index of V&R (1999), and I still can't find such a thing. I see there is a ptukey function and a qtukey function but that is as far as I got. Do I need to roll my own? -- Douglas Bates

I have a tiny R script performing two tasks, the first one of which may contain error. mammon(12)% cat z.R version options(show.error.messages = FALSE) try(b <- log("foo")) ## task 1 1 + 2 ## task 2 Running in batch mode from a solaries machine, the second task never got started; see below: mammon(13)% R BATCH --vanilla -q z.R mammon(14)% cat z.Rout >

Hi all, I'm looking for a package that will estimate Markov models and provide transition probabilities. I'm not speaking of MCMM estimation packages. Thanks in advance, David -- ======================================================================= David Kaplan, Ph.D. Professor Department of Educational Psychology University of Wisconsin - Madison Educational Sciences, Room 1061

Hi, I am using nls library df <- read.table("data.txt", header=T); library(nls); fm <- nls(y ~ a*(x+d)^(-b), df, start=list(a=max(df->y,na.rm=T)/2,b=1,d=0)); coef(fm); q(); When i am using the above routine i am getting the following error Error in nlsModel(formula, mf, start) : singular gradient matrix at initial parameter estimates Can some one help me in this.

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

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