search for: eqsystem

Displaying 8 results from an estimated 8 matches for "eqsystem".

2011 Jun 24
4
What does class "call" mean? How do I make class "formula" into a "call"?
I have a list called "tabs" that I would like to have the same structure as my list "eqSystem." The two look like they have the same format but they are different because when I look at their attributes, class(eqSystem[[1]]) is "call" but class(tabs[[1]]) is "formula". I want to have class(tabs[[1]]) as a call too. So what does "call" mean? And how do I ma...
2012 Oct 17
1
extracting and restricting coefficients
HiĀ  I want to fit two equations simultaneously EQ1<-Y1~X1+X2 EQ2<-Y2~X1+X2 eqsystem<-list(Y1HAT=EQ1,Y2HAT=EQ2) fitols<-systemfit(eqsystem, method="OLS", data=BB) How do I get coefficients for the first equation? R code How do I restrict coefficient of X2 in the first equation (say , restrict it to less thanĀ  zero). R code Your help is appreciated. Dereje [[alte...
2009 Apr 13
2
joint estimation of two poisson equations
...that unless it's entirely necessary. I'll post my solution to the list when I've worked it out. Regards, ~Owen # CODE FOR "sur" OPTION rm(list = ls()) library(Zelig) y1 = c(1,2,3,4) y2 = c(0,2,0,2) x = c(2,3,4,8) d = data.frame(cbind(y1, y2, x)) eq1 = y1 ~ x eq2 = y2 ~ x eqSystem = list (eq1, eq2) system_out = zelig(formula = eqSystem, model = "sur", data = d) summary(system_out) ----------------------------------------------------------------- # ERROR FROM REPLACING "sur" WITH "poisson" Error in switch(mode(x), `NULL` = structure(NULL, clas...
2011 Jun 28
2
How do I output all the R-squares of an SUR? summary(fitSUR$eq[[1:4]])$r.squared does not work
...library(systemfit) ## EQUATIONS: EQ1 <- Y1 ~ X1 + X2 + X4 + X7 + X10 ## equation 1 formula EQ2 <- Y2 ~ X2 + X3 + X5 + X8 + X10 ## equation 2 formula EQ3 <- Y3 ~ X5 + X6 + X7 + X9 ## equation 3 formula EQ4 <- Y4 ~ X1 + X3 + X4 + X6 + X9 ## equation 4 formula eqSystem <-list(form1 = EQ1, form2 = EQ2, form3 = EQ3, form4 = EQ4) fitSUR <- systemfit(eqSystem, method ="SUR", data=YX) ## How do I out put all the R-squares of the system without having to type a single line for each? summary(fitSUR$eq[[1]])$r.squared summary(fitSUR$eq[[2]])$r.squared...
2011 Nov 15
1
Estimating model parameters for system of equations
...ers in R for a pretty simple system of equations, but I'm having trouble. Here is the system of equations (all derivatives): eqAlgae <- (u_Amax * C_A) * (1 - (Q_Amin / Q_A)) eqQuota <- (p_max * R_V) / (K_p + R_V) - ((Q_A-Q_Amin)*u_Amax) eqResource <- -C_A * (p_max * R_V) / (K_p + R_V) eqSystem <- list(C_A = eqAlgae, Q_A = eqQuota, R_V = eqResource) I want to estimate u_Amax, Q_Amin, p_max and Q_Amin with the data I've collected using least squares. I've tried using systemfit but I'm not sure how to write out the equations (my attempt is above but that doesn't work sin...
2005 Dec 22
13
in-memory SQLite for testing
I am following the rails book, and arrived to the section about testing. The test database is configured as: test: adapter: sqlite3 database: ":memory:" and I get this error with the simple product test: % ruby test/unit/product_test.rb Loaded suite test/unit/product_test Started E/usr/local/lib/ruby/gems/1.8/gems/sqlite3-ruby-1.1.0/lib/sqlite3/ errors.rb:94:in
2006 Jan 12
10
Take Rails Studio Workshop?
In two weeks, Dave Thomas will be in Pasadena doing a three-day Rails workshop. I''m about to sign up for it. I figure three intense days with Dave and a few dozen other Rails developers will produce huge learning leaps for me. Then it occurred to me that it might be useful to ask the list if anyone has attended one of these workshops before and, if so, how valuable you found
2006 Jan 18
35
Pay rates for Rails developers
Does anyone have an idea of the rates being charged by Rails developers (or salaries for FTEs)? I''m curious to see whether rates will become comparable to those paid to more experienced Java/.NET developer types, or if rates will be lowered by the free/open source mentality, and the possible perception that Rails makes Web development "easy." Sometimes rates are more