similar to: saving estimates from a for loop for later use

Perhaps someone will have a solution to my more general problem, but here is the specific one: I used the round() function to round some estimates to 3 decimal places. I then sent put the rounded estimates in a matrix and used latex() to make a LaTeX table from them. However, in my table, there are estimtes which only have 2 decimal places. I assume that the third decimal place in these numbers

I have a dataset with 4 years of students, and normally I want to estimate things using each individual year, so I have a for loop as follows for (i in 1:4){} However, the only way I know how to calculate estimates using all four years of data is to put the estimations outside of the loop. Is there anyway to make a for loop that uses all four years at once, then uses each individual year?

To anyone who uses the Match() function in the Matching library... How do you go about deciding how many matches you will use? With my data, my standard errors generally get smaller if I use more matches. Speaking of standard errors, when correcting for heteroscedasticity, how many matches do you use (this is the Var.cal option). It seems to me that it might make sense to use the same number

I am using R for Unix and want to make some LaTeX tables. I have already played around in R for Windows and have succeeded in making tables that I want using the following code: latex(Estimates, file='out.tex', rowlabel='',digits=3) However, when I use this code in Unix, I can never find the file "out.tex". I assumed that R would send the file to whatever directory I

I would like to use a for loop to run estimations on 12 different subsets of a dataset. I have the basics of my script figured out, but I am having problems getting the loop to name some files as I would like. Here is a sketch of my code: sub.dataset <- c(101, 201) #Assume I only have two subsets which I call 101 and 201 for (i in 1:length(sub.dataset)) { #Estimation commands

Dear R users, I am trying to become a convert to R (from Stata), and I'm having some growing pains. Is there anyone would be willing to answer a few basic questions I have? Of course I am using the relevant help guides... Thanks, Brian Quinif ps. If there is anyone at UGA on this list, I'd love to hear from someone local.

When I try to install hte Econometrics view I get the following error: CRAN task view Econometrics not available in: install.views("Econometrics") I have already install the ctv package and loaded it before trying to install the above... Any ideas as to what's going on?

Consider this: > y <- matrix(1:8, ncol=2) > is.matrix(y[-c(1,2),]) [1] TRUE > is.matrix(y[-c(1,2,3),]) [1] FALSE > is.matrix(y[-c(1,2,3,4),]) [1] TRUE It seems like an inconsistent behaviour: - with 2 or more rows we have a matrix - with 1 row we do not have a matrix and - with 0 rows we have a matrix again I just stumbled on this behaviour, because I had a problem with my

I have a question about how to reference variables in a dataframe. Normally, after I have read in some Stata data using the following command all <- read.dta('all.dta') Whenever I want to use the variable "sat.vr1" in the "all" data frame, I do so using all$sat.vr1 However, I'd like to be able to use the sat.vr1 variable without the "all$" (as well

Does anyone out there use dcolumn=TRUE in the latex() function in the Hmisc library? I would like to line up the data in a latex table I'm making using latex(), but I'm having some issues with this feature. Since there is no description of it in the help, I thought that it might be incomplete or something like that. On a related note, if anyone knows how to get the end result I want

Does anyone else find that using the Var.calc option (for heteroscedasticity consistent std. errors) in Match() (from the Matching library) slows down computation of the matching estimator by a lot? I don't really understand why when I use this option it slows down so much, but for me it does significantly. I want to use the heteroscedasticity consistent std. errors in my project, but as long

Forgive my ignorance, but how I can take a look at the code for the latex function in the Hmisc library? I tried just typing "latex" but all I got was this: > latex function (object, title = first.word(deparse(substitute(object))), ...) { if (!length(oldClass(object))) oldClass(object) <- data.class(object) UseMethod("latex") } What should I do?

Dear R users, I am using the Match function to generate nearest neighbor matching estimators. On that front I am ok, but where I am having problems is getting my output into nice LaTeX tables. Here is some basic code imitating the estimation I'm doing. No problem there. library(Matching) #make up some data X1 <- matrix(rnorm(1000*5), ncol=5) Y1 <- as.vector(rnorm(1000)) Tr1 <-

Hello all, Does anyone know how to constrain/force specific coefficients when running lm()? I need to run recresid() {strucchange package} on the residuals of forecast.lm, but forecast.lm's coefficients must be determined by parameter.estimation.lm I could estimate forecast.lm without lm() and use some other kind of optimisation, but recresid() requires an object with class lm.

I am ashamed to be asking this question, but I couldn't find the solution anywhere. Searching for "if" and "R" is not very productive... I cannot get a simple if statement to work. I have data on college students. I want to make a string variable that has the names of the years. That is, when the year variable i is equal to 1, I want to have a variable called years

Dear list colleagues, I'm trying to come up with a test question for undergraduates to illustrate comparison of means from a complex survey design. The data for the example looks roughly like this: mytest<-data.frame(harper=rnorm(500, mean=60, sd=1), party=sample(c("BQ", "NDP", "Conservative", "Liberal", "None", NA), size=500,

Okay in this plea for help I''m going to repeat some of what i posted before but with a larger amount of background info in the hope that i can get a decent grip on ferret before it wriggles away.. Firstly, what does installing the acts_as_ferret plugin actually do? I install it and add it to my model and then the index is automatically generated and a few methods are added to it and

I have found out that the way to have a break apart a long caption in the way that one desires is to have a short caption for use in the list of tables that is not broken with \\ I can manually adjust the LaTeX output from the Hmisc function latex() to get the tables how I want them. However, I would like to automate the creation of these tables. Inside the latex() function, I tried this:

I have to build an array through both activerecord and also through a ferret index. I was hoping to find a way of sorting the array as i combine the two so def advanced_search(search_text, store, format, sortby) # find items in ferret index items = Item.find_by_contents(search_text) # and now find all the items from a certain store or category items_from_sql = Item.find_by_sql("SELECT *

Hello. I am getting results from sem that are not correct (that's assuming that the results from my AMOS 4.0 software are correct). sem does not vary some of the parameters substantially from their starting values, and the final estimates of those parameters as well as the model chisquare value are incorrect. I've attached some code that replicates the problem. The parameters in