Re: [R] Plot map by region using kasc and adehabitat

2011-09-22 Thread Alex Olssen
...@knmi.nl wrote:  On 09/21/2011 01:39 AM, Alex Olssen wrote: 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) I get

[R] raster plot is empty

2011-09-21 Thread 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 -

Re: [R] raster plot is empty

2011-09-21 Thread Alex Olssen
the versions OS, it's impossible to say. But my first guess (and it IS a guess) is that the server needs to be updated. Sarah On Wed, Sep 21, 2011 at 5:45 PM, Alex Olssen alex.ols...@gmail.com wrote: Dear R-help, I have a problem plotting maps using the raster package when I use R on my workplace

Re: [R] raster plot is empty

2011-09-21 Thread Alex Olssen
; returning -Inf Cheers, Mike. On Thu, Sep 22, 2011 at 8:31 AM, Alex Olssen alex.ols...@gmail.com wrote: 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

Re: [R] raster plot is empty

2011-09-21 Thread Alex Olssen
FALSE TRUE TRUE TRUE TRUEFALSE On 22 September 2011 13:09, Henrik Bengtsson h...@biostat.ucsf.edu wrote: Adding $0.01 to this: Maybe comparing outputs of capabilities() gives some clues. /Henrik On Wed, Sep 21, 2011 at 6:00 PM, Alex Olssen alex.ols...@gmail.com wrote: Yes

[R] Plot map by region using kasc and adehabitat

2011-09-20 Thread Alex Olssen
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) I get a raster map of profit per parcel for all regions. I want to plot

[R] xyplot remove ticks from top and right

2011-09-15 Thread Alex Olssen
appreciated. Kind regards, Alex Olssen __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] Emacs ESS finding all versions of R

2011-09-11 Thread Alex Olssen
on the R path, including \, /, and going to Rterm instead. Any help would be greatly appreciated. Thanks, Alex Olssen __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] Getting bootstrap statistic to work

2011-08-12 Thread Alex Olssen
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

[R] Installing package rgdal

2011-05-28 Thread Alex Olssen
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

Re: [R] maximum likelihood convergence reproducing Anderson Blundell 1982 Econometrica R vs Stata

2011-05-09 Thread Alex Olssen
estimation. Kind regards, Alex On 9 May 2011 19:04, peter dalgaard pda...@gmail.com wrote: On May 9, 2011, at 06:07 , Alex Olssen wrote: 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

Re: [R] maximum likelihood convergence reproducing Anderson Blundell 1982 Econometrica R vs Stata

2011-05-09 Thread Alex Olssen
regards, Alex On 9 May 2011 19:04, peter dalgaard  wrote: On May 9, 2011, at 06:07 , Alex Olssen wrote: 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

Re: [R] maximum likelihood convergence reproducing Anderson Blundell 1982 Econometrica R vs Stata

2011-05-09 Thread Alex Olssen
, Alex Olssen wrote: 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 and is your first coefficient -.19 or -.019? Presumably typos would be one

Re: [R] maximum likelihood convergence reproducing Anderson Blundell 1982 Econometrica R vs Stata

2011-05-08 Thread Alex Olssen
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 or more, which is exactly considering I am using

[R] maximum likelihood convergence reproducing Anderson Blundell 1982 Econometrica R vs Stata

2011-05-06 Thread Alex Olssen
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 can estimate the static model successfully in Stata but for

[R] maximum likelihood accuracy - comparison with Stata

2011-03-27 Thread Alex Olssen
(c(0,0,0,2), ols.lf, method=BFGS, hessian=T, y=y, x=x) OI - solve(p$hessian) se - sqrt(diag(OI)) se - se[1:3] beta - p$par[1:3] results - cbind(beta, se) results sigma2 - p$par[4] sigma2 Kind regards, Alex Olssen Motu Research Wellington New Zealand __ R

[R] recommended reading for manual maximum likelihood of a system of equations

2011-03-23 Thread Alex Olssen
also be a great learning experience. All suggestions will be highly appreciated! Thanks a lot, Alex Olssen Motu Economic and Public Policy Research Wellington New Zealand __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] nlysystemfit and loglikelihood values

2011-03-13 Thread Alex Olssen
the loglikelihood of a model estimated with nlsystemfit()? Kind regards, Alex Olssen __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

[R] R beginner - Error in as.vector(x, mode)

2011-03-10 Thread Alex Olssen
of type 'any' Kind regards, Alex Olssen __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained