:
> 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
FALSE TRUE TRUE TRUE TRUEFALSE
On 22 September 2011 13:09, Henrik Bengtsson 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 wrote:
>> Yes, Mik
rds(x, y, xlabel, ylabel, log) : NAs introduced by coercion
> 2: In min(x) : no non-missing arguments to min; returning Inf
> 3: In max(x) : no non-missing arguments to max; returning -Inf
>
> Cheers, Mike.
>
>
> On Thu, Sep 22, 2011 at 8:31 AM, Alex Olssen wrote:
>> D
e 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 wrote:
>> Dear R-help,
>>
>> I have a problem plotting maps using th
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 <- sys
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 sepa
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.
bin\\R")
as well as many variations 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
PLEAS
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 fun
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
uspicion that maybe our numerical
> gradients are misbehaving in cases with high local collinearity?
>
> -pd
>
> On May 9, 2011, at 13:40 , Alex Olssen wrote:
>
>> Hi Mike,
>>
>> Mike said
>> "is this it, page 1559?"
>>
>> That is the
22 | .0496865 .0348052 1.43 0.153 -.0185305 .1179034
>> g23 | .0979991 .0174397 5.62 0.000 .0638179 .1321803
>> g33 | -.0257405 .0113869 -2.26 0.024 -.0480584 -.0034226
>> --
>>
>> In R I cannot ge
et positive
log-likelihoods in these system estimation.
Kind 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 res
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 differe
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 t
(e)%*%e)/(2*sigma2))
return(-logl)
}
p <- optim(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]
sigma
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
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
as.vector(x, mode) :
cannot coerce type 'builtin' to vector 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
19 matches
Mail list logo
