[R] About MC simulation of AR1 model in R

2016-11-27 Thread yz
I want to run MC simulation of AR1(auto-regression) matrix in R, which would be as residuals in linear mixed model. AR1 matrix with the following character: ��r,��c is the auto-correlation parameter in the row and column direction. And I wrote one function in R ( under following). But I run

Re: [R] In praise of "options(warnPartialMatchDollar = TRUE)"

2016-11-27 Thread David Winsemius
> On Nov 27, 2016, at 2:03 PM, Bert Gunter wrote: > > ... >> >> After either method, you can `cbind` to that dfrm object to your heart's >> content, because the cbind.data.frame method is dispatched. >> >> -- >> David. >> > > ... But of course, this is unnecessary

Re: [R] In praise of "options(warnPartialMatchDollar = TRUE)"

2016-11-27 Thread Bert Gunter
... > > After either method, you can `cbind` to that dfrm object to your heart's > content, because the cbind.data.frame method is dispatched. > > -- > David. > ... But of course, this is unnecessary anyway, as: "The cbind data frame method is just a wrapper for data.frame(..., check.names =

Re: [R] Fitting data using Generalized Pareto Distribution

2016-11-27 Thread Rui Barradas
Hello, Why quantile(train, 0.9) ? If you use quantile(train) it seems to fit the data much better. You haven't posted a data example so I've made up one. library(eva) # needed for rgpd() library(extRemes) set.seed(1) train <- rgpd(1e3, scale = 0.9, shape = -0.4) thresh90 <-

Re: [R] In praise of "options(warnPartialMatchDollar = TRUE)"

2016-11-27 Thread David Winsemius
> On Nov 27, 2016, at 11:21 AM, Chris Evans wrote: > > Was about to reply just to sender but thought there were good byproducts of > this so, to all ... > > Many thanks Dr. Winsemius, and apologies for the HTML: tired sloppiness. My > bad. > > Aha. For the first time

Re: [R] In praise of "options(warnPartialMatchDollar = TRUE)"

2016-11-27 Thread Chris Evans
Was about to reply just to sender but thought there were good byproducts of this so, to all ... Many thanks Dr. Winsemius, and apologies for the HTML: tired sloppiness. My bad. Aha. For the first time I can really see a logic for dataFrame[['variable']] -- thanks for that. I will have to

Re: [R] In praise of "options(warnPartialMatchDollar = TRUE)"

2016-11-27 Thread David Winsemius
> On Nov 27, 2016, at 7:12 AM, Chris Evans wrote: > > I am just posting this to the list because someone else may one day waste an > hour or so because s/he has unknowingly hit a partial match failure using > "$". It's my folly that I did but I am surprised that >

Re: [R] Fitting data using Generalized Pareto Distribution

2016-11-27 Thread TicoR
Let the train be the data set consisting of numbers that I need to fit. Code is as follows: library(extRemes) thresh90 <- quantile(train, 0.90) model<-fevd(train,threshold =thresh90,type="GP") Model returns the following : Negative Log-Likelihood Value: 317.7561 Estimated parameters:

Re: [R] Fitting data using Generalized Pareto Distribution

2016-11-27 Thread Rui Barradas
Hello, Please provide us with a reproducible example. A data exampla would be nice and some working code, the code you are using to fit the data. Rui Barradas Em 27-11-2016 15:04, TicoR escreveu: I am trying to fit some data using Generalized Pareto Distribution in R using extRemes

Re: [R] Merra2 files

2016-11-27 Thread Ben Tupper
Hi, One of the best places to start a search for information is RSeek.org http://rseek.org/?q=Merra Ben > On Nov 26, 2016, at 8:04 PM, Alemu Tadesse wrote: > > Dear R users, > > I am wondering if someone has a script to download and read Merra 2 files. > > I really

[R] In praise of "options(warnPartialMatchDollar = TRUE)"

2016-11-27 Thread Chris Evans
I am just posting this to the list because someone else may one day waste an hour or so because s/he has unknowingly hit a partial match failure using "$". It's my folly that I did but I am surprised that options(warnPartialMatchDollar = TRUE) isn't the default setting. Here's a bit of

[R] Fitting data using Generalized Pareto Distribution

2016-11-27 Thread TicoR
I am trying to fit some data using Generalized Pareto Distribution in R using extRemes package(https://cran.r-project.org/web/packages/extRemes) I am able to get the parameters for the distribution. How would I get the simulated values for the model using the parameters? [[alternative

[R] R Course in Dublin (January 30th-February 1st, 2017) Intoductory -> Modern

2016-11-27 Thread Antony Unwin
An R course from introductory to modern will be given by Louis Aslett (Oxford University, author of the packages PhaseType and ReliabilityTheory) and Antony Unwin (author of the book “Graphical Data Analysis with R” CRC Press 2015 http://www.gradaanwr.net ). The

Re: [R] Partial Fraction Decomposition

2016-11-27 Thread Rolf Turner
On 27/11/16 22:09, WRAY NICHOLAS wrote: I am not aware that R really does symbolic manipulation in the way which you want, but I would have thought that this process would make a very pretty little problem to solve for you yourself, in that you could write a programme to analyse the original

Re: [R] Partial Fraction Decomposition

2016-11-27 Thread WRAY NICHOLAS
I am not aware that R really does symbolic manipulation in the way which you want, but I would have thought that this process would make a very pretty little problem to solve for you yourself, in that you could write a programme to analyse the original expression as a string of characters, whose