[R] what is f_eval? It has broken the HH package

2017-01-12 Thread Richard M. Heiberger
I am preparing for the new semester and have downloaded and installed R-3.3.2 for Macintosh and Windows, and then the HH package and its dependencies for both. Everything fresh off CRAN. All works correctly on the Macintosh. On the Windows 10 I am getting an error > library(HH) Loading required p

Re: [R] Question about Cubist Model

2017-01-12 Thread Mxkuhn
> On Jan 12, 2017, at 5:37 PM, Lorenzo Isella wrote: > > Dear All, > I am fine tuning a Cubist model (see > https://cran.r-project.org/web/packages/Cubist/index.html). > I am a bit puzzled by its output. On a dataset which contains 275 > cases, I get non mutually exclusive rules. > E.g., in the

Re: [R] for loop in R

2017-01-12 Thread jennifer . sheng2002
That is right, Bob. Only one loop for now, since I do not know how to set up the 2nd loop. Any advice from the community? Thank you! Sent from my iPhone > On Jan 12, 2017, at 4:49 PM, Robert Sherry wrote: > > I only see one for loop in your code. I am wondering if you want a second for > l

Re: [R] for loop in R

2017-01-12 Thread Robert Sherry
I only see one for loop in your code. I am wondering if you want a second for loop based upon the length of newdata. I would also think that you do not need the second call to set.seed. Bob On 1/12/2017 4:44 PM, Jennifer Sheng wrote: Dear friends, I am working on a double loop using for. On

[R] Question about Cubist Model

2017-01-12 Thread Lorenzo Isella
Dear All, I am fine tuning a Cubist model (see https://cran.r-project.org/web/packages/Cubist/index.html). I am a bit puzzled by its output. On a dataset which contains 275 cases, I get non mutually exclusive rules. E.g., in the output below, rules 2 and 3 cover all the 275 cases of the data set a

[R] Question about Cubist Model

2017-01-12 Thread Lorenzo Isella
Dear All, I am fine tuning a Cubist model (see https://cran.r-project.org/web/packages/Cubist/index.html). I am a bit puzzled by its output. On a dataset which contains 275 cases, I get non mutually exclusive rules. E.g., in the output below, rules 2 and 3 cover all the 275 cases of the data set a

[R] for loop in R

2017-01-12 Thread Jennifer Sheng
Dear friends, I am working on a double loop using for. One level of loop is to predict N times for each subject, and the second level is to predict M times for the every subject, one subject after one subject. Please note every subject have different N or M rows of data. Any advice? Thank yo

Re: [R] Qvalue package: I am getting back 1, 000 q values when I only want 1 q value.

2017-01-12 Thread Jim Lemon
Hi Tom, >From a quick scan of the docs, I think you are looking for qobj$pi0. The vector qobj$qvalue seems to be the local false discovery rate for each of your randomizations. Note that the manual implies that the p values are those of multiple comparisons within a data set, not randomizations of

Re: [R] Match ISO 8601 week-of-year numbers to month-of-year numbers on Windows with German locale

2017-01-12 Thread Nordlund, Dan (DSHS/RDA)
See comments inline. > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Bob Rudis > Sent: Thursday, January 12, 2017 12:41 PM > To: David Winsemius > Cc: r-help@r-project.org > Subject: Re: [R] Match ISO 8601 week-of-year numbers to month-of-year > numbe

Re: [R] Match ISO 8601 week-of-year numbers to month-of-year numbers on Windows with German locale

2017-01-12 Thread Bob Rudis
Aye, but this: some_dates <- as.POSIXct(c("2015-12-24", "2015-12-31", "2016-01-01", "2016-01-08")) (year_week <- format(some_dates, "%Y-%U")) ## [1] "2015-51" "2015-52" "2016-00" "2016-01" (year_week_day <- sprintf("%s-1", year_week)) ## [1] "2015-51-1" "2015-52-1" "2016-00-1" "2016-01

Re: [R] Match ISO 8601 week-of-year numbers to month-of-year numbers on Windows with German locale

2017-01-12 Thread David Winsemius
> On Jan 12, 2017, at 8:14 AM, Janko Thyson wrote: > > Dear list, > > I'm experiencing problems with converting strings of the format > "-" (e.g. 2016-01, 2016-52) to proper POSIX dates which (I > think) I need in order to retrieve the month-of-the-year number. > > Simpler put: I'd like to

[R] [R-pkgs] new package: arsenal

2017-01-12 Thread Sinnwell, Jason P.
Dear useRs, We are pleased to announce a new package that has recently appeared on CRAN, called "arsenal" (version 0.1.2: https://CRAN.R-project.org/package=arsenal). The package is a toolkit for statistical summaries, and is streamlined to work wit

[R] Qvalue package: I am getting back 1, 000 q values when I only want 1 q value.

2017-01-12 Thread Thomas Ryan
Hi all, I'm wondering if someone could put me on the right path to using the "qvalue" package correctly. I have an original p value from an analysis, and I've done 1,000 randomisations of the data set. So I now have an original P value and 1,000 random p values. I want to work out the false discov

[R] Match ISO 8601 week-of-year numbers to month-of-year numbers on Windows with German locale

2017-01-12 Thread Janko Thyson
Dear list, I'm experiencing problems with converting strings of the format "-" (e.g. 2016-01, 2016-52) to proper POSIX dates which (I think) I need in order to retrieve the month-of-the-year number. Simpler put: I'd like to match week-of-the-year numbers to month-of-the-year numbers. Ideally,

Re: [R] installing rgl

2017-01-12 Thread ProfJCNash
I had this problem this week in Linux Mint (debian/ubuntu based) and needed to install some libgl* and libglu* packages. A search for "rgl glu.h" and "rgl gl.h" found the appropriate suggestions, though I probably installed a couple of unnecessary packages too. JN On 2017-01-11 03:14 PM, Dunc

Re: [R] Problem installing R 2.12.1 on Mac OS X Yosemite 10.10.5

2017-01-12 Thread Martin Maechler
> David Winsemius > on Wed, 11 Jan 2017 22:45:38 -0800 writes: >> On Jan 11, 2017, at 10:20 PM, David Winsemius wrote: >> >> >>> On Jan 11, 2017, at 4:15 PM, Florence Lui wrote: >>> >>> I'm trying to download version 2.12.1 of R in order for it to be com

[R] [R-pkgs] Announcing new package expint

2017-01-12 Thread Vincent Goulet
I am pleased to announce the immediate availability on CRAN of the new package expint: https://cran.r-project.org/package=expint. The exponential integral E_1(x) = int_x^\Inf exp(-t)/t dt, x real and the incomplete gamma function G(a, x) = \int_x^\Inf t^{a-1} exp(-t) dt, x > 0, a real are