similar to: R beginner - Error in as.vector(x, mode)

Hello R users, I am trying to estimate the parameters of a bimodal normal distribution using moments matching, so I have to solve a non-linear system of equations. How can I solve the following simple example? x^2 - y^2 = 6 x ? y = 3 I heard about nlsystemfit, but I don?t know how to run it exactly. I have tried the following code, but it doesn?t really work: f1 <-y~ x[1]^2-x[2]^2-6 f2

Hello, I want to solve a nonlinear 3SLS problem with "nlsystemfit()". The equations are of the form y_it = f_i(x,t,theta) The functions f_i(.) have to be formulated as R-functions. When invoking "nlsystemfit()" I get the error Error in deriv.formula(eqns[[i]], names(parmnames)) : Function 'f1' is not in the derivatives table

I've got a package I would like to send out for testing and have noticed a "problem" with the \tabular portion of the help "compiler" the code for inserting vertical lines, in the tabular environment, causes the build to choke, but not the check command. C:\>rcmd build --binary -docs=all nlsystemfit \tabular{|l|c|c|c|}{ \hline Method \tab Instruments \tab

Hello everyone, I have estimated system of three linear equations with one non-linear restrictions with nlsystemfit. I was wondering how I can calculate the R-squared (or some alternative coefficient of determination) for the whole system. This is automatically given by linear systemfit but not by nlsystemfit. I can get the values for each of the equations separately, but apparently not for

Hi, I would like to estimate a quadratic almost ideal demand system in R which is estimated usually by nonlinear seemingly unrelated regression. But there is no such function in R yet but it is readily available in STATA (nlsur), see B. Poi (2008): Demand-system estimation: Update, Stata Journal 8(4). Now I am thinking, what is quicker learning to "program" STATA which seems not really

Dear R-help, The documentation for systemfit shows that logLik() can be used to obtain loglikelihood values from linear systems estimated by systemfit(). It seems to me that logLik() cannot be used for nlsystemfit(). Does anyone know of any other packages that might let me obtain the loglikelihood of a model estimated with nlsystemfit()? Kind regards, Alex Olssen

Hi everybody, I'm a newbye with lots of problems :). I'm trying to use nlsystemfit, but I recieve two error messages whose origin that I don't understand. 1) When I try to reproduce the example reported in the systemfit package manual, that is library( systemfit ) data( ppine ) hg.formula <- hg ~ exp( h0 + h1*log(tht) + h2*tht^2 + h3*elev + h4*cr) dg.formula <- dg ~ exp( d0

I've been working on a new package and I have a few questions regarding the behaviour of the nlm function. I've been (for better or worse) using the nlm function to fit a linear model without suppling the hessian or gradient attributes in the objective function. I'm curious as to why the nlm requires 31 iterations (for the linear model), and then it doesn't work when I try to add

Dear R-help, Apologies if this isn't exactly the right place for this question. I am trying to run R in emacs using ESS. I have done this successfully many times - it works right now on my home computer in windows 7 and in ubuntu. The problem is that when I start emacs the regular splash screen doesn't show. In the mini-buffer I get the message "Finding all versions of R on your

I'm using nls to fit periodic gene-expression data to sine waves. I need to set the upper and lower boundaries, because I do not want any negative phase and amplitude solutions. This means that I have to use the "port" algorithm. The problem is, that depending on what start value I choose for phase, the fit works for some cases, but not for others. In the example below, the fit works

Hey, I want to implement a structural model with the package systemfit with some linear and nonlinear constraints. How to implement linear restrictions is clear. Does anybody know how to set up nonlinear restrictions in the systemfit packages. For example: beta1 = beta2-(beta4/beta6) I look forward to your reply -- View this message in context:

Dear R-help, I am drawing a graph using xyplot. The axis labels only appear on the left and bottom axes - I used the option xyplot(scale = list(alternating = 1)) However while the labels no longer appear on the top and right border of the graph the ticks do! Any help would be greatly appreciated. Kind regards, Alex Olssen

Dear R-help, I have a problem plotting maps using the raster package when I use R on my workplace server, but not when I use my own desktop. I suspect the server version must be missing something and was wondering if anyone would have any ideas. The problem occurs after library(raster) f <- system.file("external/test.grd", package="raster") r <- raster(f) plot(f) On

Hi everyone, I am looking to do some manual maximum likelihood estimation in R. I have done a lot of work in Stata and so I have been using output comparisons to get a handle on what is happening. I estimated a simple linear model in R with lm() and also my own maximum likelihood program. I then compared the output with Stata. Two things jumped out at me. Firstly, in Stata my coefficient

Hi all, I'm trying to estimate model parameters 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,

Hi R-help, I am trying to implement a nonparametric bootstrap to find the standard errors of a simple statistics - the ratio of two scalars. I am having difficulty getting boot() to work correctly. I code a function to create the ratio of the relevant scalars. theta(data, i). When I call the function for my data and every observation appearing once, theta(test, c(1)), I get the correct

Hello R-sig-debian and (hopefully) Dirk: On Debian wheezy, I have the R packaging that CRAN (you) provide. I run into a little trouble while trying to fiddle with alternative BLAS. I know you and I went around on this last year and I think perhaps I've found something wrong in the framework, or I've just done something wrong. I installed the packages openblas-base and openblas-dev, and

Hi, I would like to estimate a linear expendire with Systemfit package. (method: "SUR") As someone could show me how to define the equations? Thanks. -- Think before you print ! ********************************************************************** Disclaimer: This e-mail may contain confidential informa...{{dropped:9}}

Dear R-helpers, I am trying to install the package "rgdal" using the command install.packages("rgdal") in R. I get the following error "Error: proj_api.h not found. If the PROJ.4 library is installed in a non-standard location, use --configure-args='--with-proj-include=/opt/local/include' for example, replacing /opt/local/* with appropriate values for your

hi, does someone tried bridge and or router functionality of linux kernel on multiprocessor system. Does multiprocessor system increase network performance or give ability to process higher traffic.. What about highest possible troughput someone have achieved ? etc.. tia _______________________________________________ LARTC mailing list / LARTC@mailman.ds9a.nl