1) you can try to make a symbolic link to libgfortran.so.3 and call it
libgfortran.so.1 Sometimes it works.
2) you can try to downgrade your version of gfortran
3) you can ask to the package maintainer to recompile FEAR with a more
modern gfortran.
Hope it helps
Running R version 2.13.1 under Ubuntu 10.10 (kernel version
2.6.35-30-generic) on a x86-64 laptop with gfortran 4.4, I'm trying to
install and use the FEAR package from
http://www.clemson.edu/economics/faculty/wilson/Software/FEAR/fear.html.
I downloaded the the FEAR 1.15 package for 64-bit Linux (
Sys.setenv(TZ="Etc/GMT+8") will use standard time only for time calculations.
---
Jeff Newmiller The . . Go Live...
DCN: Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (Solar/Batteries
On 25 September 2011 00:16, Felipe Nunes wrote:
> Hi Arne,
> my problem persists. I am still using censReg [version - 0.5-7] to run a
> random effects model in my data (>50,000 cases), but I always get the
> message.
> tob7 <- censReg(transfers.cap ~ pt.pt + psdb.pt + pt.opp + pt.coa + psdb.coa
>
On 24 September 2011 23:58, felipnunes wrote:
> Hi guys,
>
> did you figure out a way to solve your problem, Igor?
>
> I'm still trying to understand what's happening to my estimations. My data
> set is also big (>50,000) and I get the same error:
>
> Error in maxNRCompute(fn = logLikAttr, fnOrig
On Sat, Sep 24, 2011 at 4:55 PM, Felipe Nunes wrote:
> That was a great idea, Josh!
> I found a problem even when I tried a simple model. Do you have any idea
> why?
I forgot you were using a pdata.framecensReg hands off estimation
in part to plm, which wants between and within effects. Even
That was a great idea, Josh!
I found a problem even when I tried a simple model. Do you have any idea
why?
tob6 <- censReg(transfers.cap ~ 1, left=0, right=Inf, method="BHHH",
nGHQ=20, iterlim=1, data = pdata2)
Error in uniqval[as.character(effect), , drop = F] :
incorrect number of dimens
Hi guys,
did you figure out a way to solve your problem, Igor?
I'm still trying to understand what's happening to my estimations. My data
set is also big (>50,000) and I get the same error:
Error in maxNRCompute(fn = logLikAttr, fnOrig = fn, gradOrig = grad,
hessOrig = hess, :
NA in the ini
Hello Simone,
Given that your response variable is binary and, consequently, you should
use generalized models, just occurs to me a comparison between a Generalized
Linear Model (the model without the random effect) and a Generalized Linear
Mixed Model (the model with the random effect).
You coul
Hello all,
After beating my head against the wall for a few hours, I give up.
I have two problems.
1. If I use seq() to generate a POSIXct series that crosses a daylight savings
to stardard time change I get two values for the first day of change. The time
change is Nov 6,2011. The value for No
Or just ask let the second plot know to add it to the existing plot:
set.seed(2)
plot(ecdf(rnorm(50)))
plot(ecdf(rnorm(50,,2)), add = TRUE, pch = 2)
legend(-2.9, .95, legend = c("First ECDF", "Second ECDF"), pch = c(16, 2))
note the limits are set by the first plot, so watch out for that. You
ma
Can you fit a simpler model successfully? I would play around with
all your RHS variables to try to see if there is any particular one
that causes the estimation problems. Two ways I would do this would
be to drop them one at a time, or to start from an intercept only
model (i.e., transfers.cap ~
Hi Arne,
my problem persists. I am still using censReg [version - 0.5-7] to run a
random effects model in my data (>50,000 cases), but I always get the
message.
tob7 <- censReg(transfers.cap ~ pt.pt + psdb.pt + pt.opp + pt.coa + psdb.coa
+ pib.cap + transfers.cap.lag + pib.cap + ifdm + log(popula
Joseph Park att.net> writes:
>
> Perhaps this question is inappropriate or posted to the wrong list?
>
> Any guidance would be appreciated. Thanks.
>
[snip]
This seems like a perfectly reasonable question.
Sometimes people just don't bite.
Ben Bolker
_
Greetings,
I am trying to use Reference Class methods in Snowfall, using R 2.12.1
on Ubuntu Natty. Using then directly seems to work (stanza 2 below),
but using them indirectly does not (stanza 3 below). I get an "attempt
to apply non-function" error.
In addition to exporting the instance
What version of R are you using?
x = ecdf(rnorm(50)); y = ecdf(rnorm(50))
plot(x); lines(y, col=2)
works for me on R 2.13.1. Specifically, there is a lines.stepfun method
(albeit non-visible) which handles this.
Michael
On Sat, Sep 24, 2011 at 12:23 PM, Adel ESSAFI wrote:
> Hello,
>
> is ot
Perhaps this question is inappropriate or posted to the wrong list?
Any guidance would be appreciated. Thanks.
--- On Fri, 9/23/11, Joseph Park wrote:
> From: Joseph Park
> Subject: Cross Spectrum : Conversion of 2-D spectrum into a single complex
> array
> To: r-help@r-project.org
> Date: F
Hello,
is ot possible to draw two ecdf of vectors (say s1 and s2) on the same
figire with R.
plot function draws a new plot and there is no function like points or lines
to draw a second ecdf on the figure.
Regards
--
*PhD candidate in Computer Science
Address
3 avenue lamine, cité ezzahra, S
For my part, knowing all of my (mostly bad) habits, I would only donate if I
could do it online. And I would be happy to donate!
Cheers
David Cross
d.cr...@tcu.edu
www.davidcross.us
On Sep 24, 2011, at 10:43 AM, Joshua Wiley wrote:
> Dear R Users,
>
> I know issues like this have come up
Dear R Users,
I know issues like this have come up in the past. I am wondering how
many people would be interested in an online mechanism for donating to
the R Foundation or for becoming supporting members. Right now, the
link for it is somewhat buried on the R home page, and once you get
there,
It is not simply the answer, it is a list
str(integrate(dnorm, -1.96, 1.96))
List of 5
$ value : num 0.95
$ abs.error : num 1.05e-11
$ subdivisions: int 1
$ message : chr "OK"
$ call: language integrate(f = dnorm, lower = -1.96, upper = 1.96)
- attr(*, "class")= chr "i
Any answer to this?
I really need to compare a mixed model with binomial error against the same
model without the random effect. I would use anova() but I don't know how to
specify both models in order to make them comparable.
Thanks for any answer
Simone
--
View this message in context:
http:
Hi,
I need to do a calculation in R which involves adding the estimate of a
numerical integration with another number - however, if I add the number to
the output of the intergrate function (which when run on its own provides me
with a number with error), I get the message 'non-numeric argument to
2011/9/24 Muhammad Rahiz :
> Dear all,
>
> I'm having issues with the installation of the ncdf package. It returns a
> non-zero exit status. Can anyone suggest what I should do next? FYI, I do
> not have problems installing other packages.
>[...]
> /usr/bin/ld: /usr/local/lib/libnetcdf.a(attr.o): r
Please read instructions CAREFULLY.
This is the R-help mailing list, not R-help-request. People on this list are
not interested in an empty mail, much less one that is followed by an entire
digest message!
On Sep 24, 2011, at 13:56 , Mathew Brown wrote:
>
> Mathew Brown
> Institute of Biocl
Mathew Brown
Institute of Bioclimatology
University of Göttingen
Büsgenweg 2
37077 Göttingen, Germany
t: +49 551 39 9359
mathew.br...@forst.uni-goettingen.de
On 9/24/2011 6:00 PM, r-help-requ...@r-project.org wrote:
Send R-help mailing list submissions to
r-help@r-project.org
To subsc
Dear all,
I am not able to draw a 45 degrees line in a zoo plot. The time index is a
sequence of 15 min bins from 7 in the morning to 18:45 in the evening.
Thanks,
Marco
--
View this message in context:
http://r.789695.n4.nabble.com/45-degress-line-in-zoo-plot-tp3839182p3839182.html
Sent from t
Dear R-Users!
I am using the svm function (e1071 package) to classify two groups using a
set of 180 indicator variables. Now I am confused about the cross-validation
procedure.
(A) On one hand I use the setting cross=10 in the svm function to run 10
cross-validation iterations and to get an estim
I will be out of the office Sep 32-28, 2011, with limited email access.
For immediate help, please call the JPSM main number, 301-314-7911.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Szeptember 12-től 26-ig irodán kívül vagyok, és az emailjeimet nem érem el.
Sürgős esetben kérem forduljon Kárpáti Edithez (karpati.e...@gyemszi.hu).
Üdvözlettel,
Mihalicza Péter
I will be out of the office from 12 till 26 September with no access to my
emails.
In urgent cases please contact
Is there any way to use expression() in legend labels with ggplot2?
It seems that things like
> scale_shape_manual(value = c(
> x = expression(italic(x)),
> y = expression(italic(y))
> ))
don't work.
Thanks very much :)
--
Using GPG/PGP? Please get my current public key (ID: 0xAEF6A134,
Dear Dr. Snow,
Thank you for your reply.
1. Are you doing the 2 sample KS test? Comparing if 2 samples come from the
same distribution? -Yes, I am doing 2-sample KS test
2. With 3,000 points you will still likely have power to find meaningless
differences, what exactly are you trying
32 matches
Mail list logo
