similar to: Displaying Iteration Count

Hello, I created an array to hold the results of a series of simulations I'm running: d.eta <- array(0,dim=c(3,3,200)) <simulation goes here and populates the array but it's not important> Then I tried to save the results using this: save(d.eta,file="D:/Simulation Results/sim 9-23-11 deta") When I later tried to reload them using this: d.eta <-

Hello, I have a longitudinal dataset where each individual has a different number of entries. Thus, it is of the following structure: x <- runif(12) id.var <- factor(c(rep("D1",4),rep("D2",2),rep("D3",3),rep("D4",3))) dat <- as.data.frame(x) dat$id.var <- id.var dat > dat x id.var 1 0.9611269 D1 2 0.6738606 D1 3

Hello, Is there a way to get a single TRUE or FALSE statement from comparing two vectors? For example, c(1,2,3) == c(1,2,3) produces TRUE TRUE TRUE where I would like it to produce only TRUE for use in an if statement. Likewise, when two vectors are not exactly identical (in all elements) I would like a single FALSE result, as opposed to c(1,2,3) == c(1,2,5) TRUE TRUE FALSE Any ideas?

Hello, I hate to ask a question which is directly addressed in the documentation, but can someone please give me an example of how to change the maximum number of iterations used by lmer. I'm having a hard time understanding this: control a list of control parameters. See below for details. control a named list of control parameters for the estimation algorithm, specifying only

Hello, I'm clearly confused about the merge function. In the following r <- merge(x,y,all.x=T,all.y=F) my y vector has only unique values (no duplicates). So I don't understand how this can ever generate an r which is of greater length than x. I thought the default behavior was only matching rows are included, but that using all.x=T included rows with unmatched x's as well. If

Hello, I'm attempting to install the splm package from R-Forge. https://r-forge.r-project.org/R/?group_id=352 The page says, "In order to successfully install the packages provided on R-Forge, you have to switch to the most recent version of R..." It later says "To install this package directly within R type: install.packages("splm",

Hello, I am using {plm} to estimate panel models. I want to estimate a model that includes fixed effects for time and individual, but has a random individual effect for the coefficient on the independent variable. That is, I would like to estimate the model: Y_it = a_i + a_t + B_i * X_it + e_it Where i denotes individuals, t denotes time, X is my independent variable, and B (beta) is the

Hello, A colleague is trying to do a fairly complicated power analysis for a project. The project would be evaluating random assignment to one of three conditions within each of 8 sites. The dependent variable would be binary (we do not care at this point whether it would be analyzed with logit or probit). We can simulate the data in a bunch of iterations and do analyses in each iteration, but

Hello, I'd like to plot an empirical cumulative distribution function, except instead of the fraction of values < x, I'd like the fraction of values > x. I think this can be done using the ecdf function in {Hmisc}. I installed the package and loaded it. However, when following the example given in the documentation, I get an error: x <- rnorm(100) ecdf(x,what='1-F')

Hi, I would like to estimate a VAR of the form: Ay_t = By_t-1 + Cy_t-2 + ... + Dx_t + e_t Where A is a non-diagonal matrix of coefficients, B and C are matricies of coefficients and D is a matrix of coefficients for the exogenous variables. I don't think the package {vars} can do this because I want to include contemporaneous cross-variable impacts. So I want y1_t to affect y2_t and I

I understand that Rcmdr needs the tcltk package in order to run. However, when I attempt to using the library (Rcmdr) command, I get the message below. I have attempted to start X in R as well as using starting it from the Utilities folder on my Mac to no avail. I can't figure out if there is an issue with X11 or with R itself (I just updated to 1Is there a terminal command I can use

Hello, I'm using dmnorm from the package {mnormt} and getting strange results. First, according to the documentation, dmnorm should return a vector of densities, and I'm only getting one value returned (which is what I would expect). I've been interpreting this as the joint density of all values in the x vector (which is what I want). Should a vector of densities be returned, and if

I am trying to use winProgressBar, however I find that although setWinProgressBar updates the value on the screen, getWinProgressBar does not return this value. E.g > pb <- winProgressBar() > setWinProgressBar(pb, 0.2) > getWinProgressBar(pb) [1] 0 I tried the same with tkProgressBar, and it is working as I would expect from the documentation. > pb <- tkProgressBar() >

Hello, I'm working on an structural VAR using the var command to estimate and the svar command on the resultant object (package: vars). I want to constrain coefficients to equal one another, but that value to be estimated. So for the A matrix, I want A[2,1]=A[1,2] to be my constraints. Can this be done with this package? If so, how? If not, is there another package that it might be done with

This problem can be seen by the following commands: > pb <- winProgressBar(max=1000, label='0') > b <- 1 > setWinProgressBar(pb, b, label=b) This set of commands (on windows of course, XP in this case) causes R to crash. This is not strictly a bug since the documentation states that the label argument should be a character string and using as.character(b) does work

I am trying to do a tutorial improvement on the wiki but I am having some problems with the auto formatting. Using httdp.conf produces http://httpd.conf when filtered. Using *.so produces .so and bold text there after. Is there any shorthand for the user for making text appear as it should without the linking and re-writes? Or does the admin have to change the filters in Instiki? Tesla

Dear all, I have a column within a dataframe of values which range between 1 and 2. I want to display graphically the distribution of these values (i.e. are they clustered towards either exteme? Or spread evenly?). What is a good way of doing this in R? I've tried a few things, including using the 'hist' command, but receive the following error message: > hist(urban.long[3])

Hi, Rsync is still hanging under cygwin for me. I'm using rsync 2.6.0 and i've tried cygwin-1.5.6-1 as well as a snapshot of 1.5.7. I've noticed that some people are having success. I shouldn't have to do anything interesting in order to make this work, right? *********************************************************************** This message is intended only for the use of

I want to use ActiveRecord outside of Rails to do processing on large tables. I want to be able to process the rows as they are retrieved rather than having to store the results in an intermediate array. Any suggestions? Looks like I could hack the adapter to provide an accessor to the underlying DBI handle. Anything cleaner than that?

Firewall users, My apologies as I''m not on this list, so please respond directly as well as to the list. I did try to search the archives and didn''t find any hits, although the search did not like searching for terms with underscores in them (both clear_firewall and ip_forward). I was trying to understand why, when running shorewall stop, even though it echoes IP