Search Archives

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

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

Dear R-help, I have a raster map which has a measure of profitability of land by parcel over several regions of geographic aggregation (think of counties). This data is stored in an ASCII file. If I type plot(profit...

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

Hi everyone, I am looking for some recommended reading. I want to read up on the estimation systems of linear equations using maximum likelihood? I have a strongly applied bias, I want to be able to do such estimation myself...

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

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

Hi everyone, I am new to R and keep getting the message Error in as.vector(x, mode) while trying to run nlsystemfit. Below is my exact code. The data is in Stata format because I only recently swapped to...

Dear Sarah, I am not sure how to make a reproducible example. The 3 lines of code are exactly what I type in either system. The code can be reproduced, but the result is unlikely to be reproduced - for example...

Peter said "Ahem! You might get us interested in your problem, but not to the level that we are going to install Stata and Tsp and actually dig out and study the scientific paper you are talking about. Please cite...

Hi there, I thought it would be useful to answer my own question for anybody else searching Google. An easy way to split a national map into smaller maps at lower levels of aggregation is to use the function mask...

Unfortunately the capabilities are identical too. For both systems I get jpeg png tiff tcltk X11 aqua http/ftp sockets TRUE TRUE TRUE TRUE FALSE FALSE TRUE TRUE libxml fifo cledit iconv NLS profmem cairo TRUE FALSE TRUE TRUE TRUE...

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

Thank you all for your input. Unfortunately my problem is not yet resolved. Before I respond to individual comments I make a clarification: In Stata, using the same likelihood function as above, I can reproduce EXACTLY (to 3 decimal places...

Wow that is really interesting, Sorry I was asleep when you emailed these. And yes, of course, I had been trying to implement model 18, not 18s, that was a typo, sorry. I will have a look at the code...

Yes, Mike you are correct. In fact I did use plot(r) The code is exactly the same in both cases - the axis labels have the same numbers, same with the legend - there is just no picture. sessionInfo() for my...

Hi Mike, Mike said "is this it, page 1559?" That is the front page yes, page 15*6*9 has the table, of which the model labelled 18s is the one I replicated. "did you actually cut/paste code anywhere...

Dear R-help, I am trying to reproduce some results presented in a paper by Anderson and Blundell in 1982 in Econometrica using R. The estimation I want to reproduce concerns maximum likelihood estimation of a singular equation system. I...

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