similar to: Multivariate time-series

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

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

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, 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, I have been searching and can't seem to find a simple command that will allow me to sample from a multivariate normal distribution with known covariance matrix. I am likely missing something. Thanks in advance. David -- =========================================================================== David Kaplan, Ph.D. Professor Department of Educational Psychology University of

Hi all, I'm running White's test to correct for non-constant error variance and I am using the Anova function in the package CAR. My command structure is > Anova(scireg3, white.adjust="hc3") where scireg3 is an object from lm. I get the message "Error in SS[i] <- SS.term(names[i]) : nothing to replace with" What does this mean and how do I fix it. Thanks

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

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

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,

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

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

Hi all, I'm looking for a program that will take the predicted probabilities from a logistic regression using glm{stats}, dichotomize them according to a threshold that I can control, and then use them to form sensitivity, specificity, false pos and false neg rates. Thanks in advance. David -- ======================================================================= David Kaplan, Ph.D.

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, I'm getting this error message using Joe Shafer's "NORM" package. Error in storage.mode(x) <- "double" : (list) object cannot be coerced to type 'double' I'm not sure what this means. I get this message when running s <- prelim.norm(filename) prelim.norm is used to develop some summary statistics, etc. before implementing the EM

Hi, Are there any incompatibilities with R 2.3.1 and SPSS 14.0 with regard to the read.spss command? 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 email: dkaplan

Hi all, I'm moving to a Mac platform and have installed R and all seems fine. My question concerns moving my R objects over. On my pc, which file contains my R objects and once I copy that, where do I copy it within the Mac? I hope that was clear. Thanks in advance, David -- =============================================================== David Kaplan, Ph.D. Professor