[R] Bayesian Analysis in GJR-GARCH (p, d) model with Student-t innovations

2018-01-18 Thread RESA MAE SANGCO
Good day Ma'am/Sir, I am Resa Mae R. Sangco a Master of Statistics student from the MSU- Iligan Institute of Technology located in Iligan City, Philippines. I am currently doing my thesis entitled “Bayesian Analysis in GJR-GARCH (p,d) model with Student-t innovations". In finding my posterior

Re: [R] error while loading ggplot2

2018-01-18 Thread shijin mathew via R-help
Thnaks Jeff. Below is the session info you requested. R version 3.4.3 (2017-11-30) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows >= 8 x64 (build 9200) Matrix products: default locale: [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252 [3]

Re: [R] reading lisp file in R

2018-01-18 Thread Martin Møller Skarbiniks Pedersen
Here are the beginning of a R program. I hope it can help you writing the rest of the program. Regards Martin M. S. Pedersen filename <- "university.data" lines <- readLines(filename) first <- T for (ALine in lines) { ALine <- sub("^ +","",ALine) ALine <- sub(")","",ALine, fixed = T)

[ESS] Sweave with ESS and AUCTeX

2018-01-18 Thread Sparapani, Rodney
Hi Gang: I’ve run into a problem that I can’t seem to figure out. I have a vignette that I am writing based on jss.sty (the JSS style) which comes with R. From the command line everything works as expected, i.e., “R CMD Sweave --pdf computing.Rnw” produces computing.pdf without issue.

Re: [R] Split charts with ggplot2, tidyquant

2018-01-18 Thread Joshua Ulrich
If you don't want to wait for a ggplot2 solution, here are two alternatives you can use right now: chartSeries(SPYxts) # or (with xts > 0.10 plot(SPYxts$SPY.Close) addSeries(SPYxts$SPY.Volume, type = "h") You might also try autoplot.zoo(), though I've never used it. On Thu, Jan 18, 2018 at

Re: [R] Perform mantel test on subset of distance matrix

2018-01-18 Thread Sarah Goslee
Hi Andrew, Yes, you cannot have NA values in your matrices. Instead, you could incorporate a model matrix. See Legendre, P. & Fortin, M.J. Vegetatio (1989) 80: 107. https://doi.org/10.1007/BF00048036 for ideas. Sarah On Sun, Dec 31, 2017 at 12:55 PM, Andrew Marx wrote:

Re: [R] Time-dependent coefficients in a Cox model with categorical variants

2018-01-18 Thread Jeff Newmiller
Offlist... for your information... It is unfair to suggest that the mailing list participants are at fault for using old software. Even if the mailing list participants use email programs that can handle HTML, any email that goes through the list gets the formatting stripped, which leaves it

Re: [R] Split charts with ggplot2, tidyquant

2018-01-18 Thread Eric Berger
Hi Charlie, I am comfortable to put the data in any way that works best. Here are two possibilities: an xts and a data frame. library(quantmod) quantmod::getSymbols("SPY") # creates xts variable SPY SPYxts <- SPY[,c("SPY.Close","SPY.Volume")] SPYdf <-

Re: [R] Time-dependent coefficients in a Cox model with categorical variants

2018-01-18 Thread Therneau, Terry M., Ph.D.
First, as others have said please obey the mailing list rules and turn of First, as others have said please obey the mailing list rules and turn off html, not everyone uses an html email client. Here is your code, formatted and with line numbers added. I also fixed one error: "y" should be

Re: [R] Split charts with ggplot2, tidyquant

2018-01-18 Thread Charlie Redmon
Could you provide some information on your data structure (e.g., are the two time series in separate columns in the data)? The solution is fairly straightforward once you have the data in the right structure. And I do not think tidyquant is necessary for what you want. Best, Charlie --

Re: [R] MCMC Estimation for Four Parametric Logistic (4PL) Item Response Model

2018-01-18 Thread ProfJCNash
If you have the expression of the model, package nlsr should be able to form the Jacobian analytically for nonlinear least squares. Likelihood approaches allow for more sophisticated loss functions, but the optimization is generally much less reliable because one is working with essentially

Re: [R] MCMC Estimation for Four Parametric Logistic (4PL) Item Response Model

2018-01-18 Thread Doran, Harold
I know of no existing functions for estimating the parameters of this model using MCMC or MML. Many years ago, I wrote code to estimate this model using marginal maximum likelihood. I wrote this based on the using nlminb and gauss-hermite quadrature points from statmod. I could not find that

Re: [R] reading lisp file in R

2018-01-18 Thread peter dalgaard
Yes, and the structure is obviously case-insensitive. More troublesome is probably that there can be multiple ACADEMIC-EMPHASIS entries, which can be tricky to tidify. Also one would need to figure out what is the meaning of lines like (DEFPROP BOSTON-COLLEGE0 T DUPLICATE) -pd > On 18 Jan

[R] MCMC Estimation for Four Parametric Logistic (4PL) Item Response Model

2018-01-18 Thread ALYSSA FATMAH MASTURA
Good day Sir/Ma'am! This is Alyssa Fatmah S. Mastura taking up Master of Science in Statistics at Mindanao State University-Iligan Institute Technology (MSU-IIT), Philippines. I am currently working on my master's thesis titled "Comparing the Three Estimation Methods for the Four Parametric

Re: [R] reading lisp file in R

2018-01-18 Thread Barry Rowlingson
The file also has a bunch of email headers stuck in the middle of it: . (QUALITY-OF-LIFE SCALE:1-5 4) (ACADEMIC-EMPHASIS HEALTH-SCIENCE) ) --- --- >From lebow...@cs.columbia.edu Mon Feb 22 20:53:02 1988 Received: from zodiac by meridian (5.52/4.7) Received: from

Re: [R] effects & lme4: error since original dataframenotfoundWASeffects: error when original data frame is missing

2018-01-18 Thread Gerrit Eichner
Thanks, John, for your hint! (Unfortunately, I was not aware of this vignette, but I am glad that I seem to habe been on the right track.) Indeed very helpful, in particular of course, the warning regarding the danger of overwriting already existing objects. That danger might be reduced by

Re: [R] request for code

2018-01-18 Thread Marc Schwartz
> On Jan 18, 2018, at 7:49 AM, Anjali Karol Nair wrote: > > Hi, > > I want to convert my MATLAB programs to R studio programs. > Kindly guide on the same. Hi, Using Google with a search phrase such as "MATLAB to R" will yield a number of possible resources for you such

Re: [R] request for code

2018-01-18 Thread jim holtman
a simple Google search turns up several possible choices. There is a package 'matconv' that might serve your purposes. Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it. On Thu, Jan 18, 2018 at 7:49 AM, Anjali

[R] request for code

2018-01-18 Thread Anjali Karol Nair
Hi, I want to convert my MATLAB programs to R studio programs. Kindly guide on the same. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help

[R] Web scraping different levels of a website

2018-01-18 Thread David Jankoski
Hey Ilio, On the main website (the first link that you provided) if you right-click on the title of any entry and select Inspect Element from the menu, you will notice in the Developer Tools view that opens up that the corresponding html looks like this (example for the same link that you

Re: [R] Letters group Games-Howell post hoc in R

2018-01-18 Thread Meyners, Michael
Apologies if I missed any earlier replies - did you check multcompLetters in package {multcompView}? It allows you to get connecting letters reports (if that's what you are after, I didn't check what exactly agricolae is providing here). May have to add some manual steps to combine this with any

[R] Web scraping different levels of a website

2018-01-18 Thread Ilio Fornasero
I am web scraping a page at http://catalog.ihsn.org/index.php/catalog#_r=1890=1=100==_by=nation_order==2017==s= >From this url, I have built up a dataframe through the following code: dflist <- map(.x = 1:417, .f = function(x) { Sys.sleep(5) url <-

Re: [R] reading lisp file in R

2018-01-18 Thread Peter Crowther
That's a nice example of why Lisp is both powerful and terrifying - you're looking at a Lisp *program*, not just Lisp *data*, as Lisp makes no distinction between the two. You just read 'em in. The two definitions at the bottom are function definitions. The top one defines the def-instance