Hello I don't know whether strictly speaking this is an R-help province
question but anyway...
I am looking at the extRemes package and have some code given beneath.
When I run it on my own laptop it works fine, but when I run it on my uni
laptop (which is pretty new and I only was given a couple of months ago)
when it gets to the thresh.range line it goes into an unending loop (or
something) with the little red Stop button, and just sits there. I have R
4.0.4 on my uni laptop and R 4.0.5 on my own one. Could this be the reason?
Thanks, Nick Wray
install.packages("extRemes")
library(extRemes)
data(Fort)
names(Fort)
bmFort <- blockmaxxer(Fort, blocks = Fort$year, which="Prec")
names(bmFort)
plot(Fort$year, Fort$Prec, xlab = "Year",
ylab = "Precipitation (inches)",
cex = 1.25, cex.lab = 1.25,
col = "darkblue", bg = "lightblue", pch = 21)
points(bmFort$year, bmFort$Prec, col="darkred", cex=1.5)
# Fit a GEV distribution to annual maximum Precipitation
# in Fort Collins, Colorado, U.S.A.
fitGEV <- fevd(Prec, data = bmFort)
fitGEV
plot(fitGEV)
plot(fitGEV, "trace")
# Select a threshold for daily data from above.
threshrange.plot(Fort$Prec, c(0.1, 2), type = "PP")
# Fit the PP model to the daily data from above.
fitPP <- fevd(Prec, data = Fort, threshold = 0.395,
type = "PP", units = "inches")
fitPP
plot(fitPP)
atdf(Fort$Prec, 0.395)
extremalindex(Fort$Prec, 0.395, blocks=Fort$year)
dcFort <- decluster(Fort$Prec, 0.395, r = 9)
plot(dcFort)
[[alternative HTML version deleted]]
Hi Nick, Thanks for providing an example that doesn't seem to depend on external data! On Thu, 29 Jul 2021 09:07:12 +0100 Nick Wray <nickmwray at gmail.com> wrote:> I am looking at the extRemes packageIt's probable that there aren't many extRemes users on R-help, but if you contact the maintainer of extRemes (run maintainer('extRemes') to get their name and address), you might get help that way.> When I run it on my own laptop it works fine, but when I run it on my > uni laptop (which is pretty new and I only was given a couple of > months ago) when it gets to the thresh.range line it goes into an > unending loop (or something) with the little red Stop button, and > just sits there.If you interrupt it there, backtrace() will point you in the direction of seemingly endless loop. You can also set options(error = recover) and get a debugger prompt at the point where you interrupted the code or in any other frame. See ?browser and RShowDoc('R-lang') part 9 "Debugging" for more information.> I have R 4.0.4 on my uni laptop and R 4.0.5 on my own one. Could > this be the reason?I'm afraid it's easier for you to answer this question. It should be possible to install R on computers even without administrator access. Can you reproduce the issue with https://cloud.r-project.org/bin/windows/base/old/4.0.4/ on your own laptop or https://cloud.r-project.org/bin/windows/base/old/4.0.5/ on the university-issued one?> [[alternative HTML version deleted]]This mailing list strips away the HTML part of multipart/mixed messages, so to make sure that your messages look exactly as you type them, please compose them in plain text. -- Best regards, Ivan