Re: [R] [FORGED] Re: merging-binning data

2015-11-04 Thread Rolf Turner
I have been vaguely following this thread and have become very confused given the complications that seem to have appeared. The original question was: On Tue, Nov 3, 2015 at 10:47 AM, Alaios via R-help wrote: Dear all,I am not exactly sure on what is the proper name

Re: [R] JAGS 4.x, rjags 4.x problems | Linux

2015-11-04 Thread Yvan Richard
Hi Evan The last version of rjags on CRAN is 4.4. Have you tried it? http://cran.stat.auckland.ac.nz/web/packages/rjags/index.html I did have the same problem before updating it, but it now works on my Ubuntu with the new JAGS version. Cheers, -Yvan

Re: [R] JAGS 4.x, rjags 4.x problems | Linux

2015-11-04 Thread Yvan Richard
According to the README file in the rjags source code (in the tar.gz file), you can apparently pass options to the configure script to specify where JAGS is installed. I'm confident you can solve your problem by following the instructions in that README file. On 5 November 2015 at 10:10, Evan

[R] mysterious error message not sure if its from within package?

2015-11-04 Thread Sharon Louise Every
Hi everyone, Sorry if this is not appropriate through this email list but I wanted to check if this was something wrong with R in my computer rather than within a package I’m using. I am using SIAR - SIBER and have for a few months now. I went back to check some figures with code that I had

Re: [R] JAGS 4.x, rjags 4.x problems | Linux

2015-11-04 Thread Evan Cooch
On 11/4/2015 2:08 PM, Evan Cooch wrote: Greetings -- This has also been posted on the jags forum, but since I suspect the problem is more 'R-related' than jags, will aos post here. Decided to 'upgrade' from jags 3.x.x to 4.x.x today, on my GNU/Linux boxes (which run latest RHEL). Here are

Re: [R] JAGS 4.x, rjags 4.x problems | Linux

2015-11-04 Thread Evan Cooch
On 11/4/2015 3:31 PM, Yvan Richard wrote: Hi Evan The last version of rjags on CRAN is 4.4. Have you tried it? http://cran.stat.auckland.ac.nz/web/packages/rjags/index.html I did have the same problem before updating it, but it now works on my Ubuntu with the new JAGS version. Cheers, -Yvan

Re: [R] JAGS 4.x, rjags 4.x problems | Linux

2015-11-04 Thread Evan Cooch
I wonder if this is the problem -- on my system(s), R is installed to /usr/lib64/R, whereas JAGS and associated files is in /usr/local/lib/JAGS. I have a hunch that the problem relates to where rjags expects to find things during install/compile, given how R was compiled, and where it is

Re: [R] JAGS 4.x, rjags 4.x problems | Linux

2015-11-04 Thread Rainer Hurling
Am 04.11.2015 um 20:08 schrieb Evan Cooch: > Greetings -- > > This has also been posted on the jags forum, but since I suspect the > problem is more 'R-related' than jags, will aos post here. > > Decided to 'upgrade' from jags 3.x.x to 4.x.x today, on my GNU/Linux > boxes (which run latest

[R] JAGS 4.x, rjags 4.x problems | Linux

2015-11-04 Thread Evan Cooch
Greetings -- This has also been posted on the jags forum, but since I suspect the problem is more 'R-related' than jags, will aos post here. Decided to 'upgrade' from jags 3.x.x to 4.x.x today, on my GNU/Linux boxes (which run latest RHEL). Here are the basic details: 1\ used R 3.2.2

Re: [R] merging-binning data

2015-11-04 Thread Boris Steipe
Whatever approach is "best" to define subsets depends completely on the semantics of the data. Your approach (a fixed number of equally spaced breaks) is the right one if the absolute ranges of the data is important. It should be obvious that either the top or the bottom group could contain

[R] Error in eval(expr, envir, enclos)

2015-11-04 Thread Alaa Sindi
Hi All, I am receiving this error Error in eval(expr, envir, enclos) : could not find function “LL” the following is in a for loop and all the variables are defined and have values. Prob[i,1]=log(sigma)^((1-M[i,1]*

Re: [R] Alternatives for explicit for() loops

2015-11-04 Thread Maram SAlem
Hi Jim, Thanks a lot for replying. In fact I'm trying to run a simulation study that enables me to calculate the Bayes risk of a sampling plan selected from progressively type-II censored Weibull model. One of the steps involves evaluating the expected test time, which is a rather complicated

Re: [R] merging-binning data

2015-11-04 Thread Boris Steipe
I would transform the original numbers into integers which you can use as group labels. The row numbers of the group labels are the indexes of your values. Example: assume your input vector is dBin nGroups <- 5 # number of groups groups <- (dBin - min(dBin)) / (max(dBin) - min(dBin)) # rescale

Re: [R] merging-binning data

2015-11-04 Thread Boris Steipe
I don't understand: - where does the "label" come from? (It's not an element of your data that I see.) - what do you want to do with this "label" i.e. how does it need to be associated with the data? B. On Nov 4, 2015, at 7:57 AM, Alaios wrote: > Thanks it works great

Re: [R] merging-binning data

2015-11-04 Thread Alaios via R-help
Thanks it works great and gives me group numbers as integers and thus I can with which group the elements as needed (which (groups== 2)) Question though is how to keep also the labels for each group. For example that my first group is the [13,206) RegardsAlex On Wednesday, November 4,

Re: [R] Error survreg: Density function returned an an invalid matrix

2015-11-04 Thread Therneau, Terry M., Ph.D.
Hi, I want to perform a survival analysis using survreg procedure from survival library in R for a pareto distribution for a time variable, so I set the new distribution using the following sintax: library(foreign) library(survival) library(VGAM) mypareto <-

Re: [R] merging-binning data

2015-11-04 Thread Alaios via R-help
you are right.by labels I mean the "categories", "breaks" that my data fall in.To be part of group 2 for example you have to be in the range of [110,223) I need to keep those for my plots. Did I describe it more precisely now?Alex On Wednesday, November 4, 2015 2:09 PM, Boris Steipe

[R] Splines with more restrictions

2015-11-04 Thread BXC (Bendix Carstensen)
I would like to generate a natural spline where the slope beyond the last knot (the right outer knot) is constrained to be 0. Is that available in some package around? _ Bendix Carstensen Senior Statistician Steno Diabetes Center Clinical Epidemiology Niels Steensens Vej 2-4

Re: [R] merging-binning data

2015-11-04 Thread Alaios via R-help
Thanks for the answer. Split does not give me the indexes though but only in which group they fall in. I also need the index of the group. Is the first, the second .. group?Alex On Tuesday, November 3, 2015 5:05 PM, Ista Zahn wrote: Probably

Re: [R] merging-binning data

2015-11-04 Thread Boris Steipe
The breaks are just the min() and max() in your groups. Something like sprintf("[%5.2f,%5.2f]", min(dBin[groups==2]), max(dBin[groups==2])) ... should achieve what you need. B. On Nov 4, 2015, at 8:45 AM, Alaios wrote: > you are right. > by labels I mean the

Re: [R] Error in eval(expr, envir, enclos)

2015-11-04 Thread William Dunlap
You need to show more (e.g., the entire function, the output of traceback() after the error) to get a definitive answer, but that error message can occur if you mistakenly use parentheses instead of square brackets when subscripting LL. E.g., > LL <- 11:20 > LL[3:4] [1] 13 14 > LL(3:4)

Re: [R] CRAN - Package "activity" - how do you convert time (hh:mm) into radians in R 3.1.2 on Windows 7 OS 64bit

2015-11-04 Thread jwd
On Wed, 4 Nov 2015 09:02:46 +0100 Martin Maechler wrote: > > Bert Gunter > > on Tue, 3 Nov 2015 14:54:46 -0800 writes: > > > ... should have said: with a 24 hour "clock". > > -- Bert > > Bert Gunter > > Hmm, thank

Re: [R] merging-binning data

2015-11-04 Thread Alaios via R-help
Thanks for your comments. Actually only the last group has a single element. The first group is always "full" of members and as that it works fine. Some constant spacing between the groups would be good as well and thus I will check quantiles. Thanks for the great support and time invested on

Re: [R] Error in eval(expr, envir, enclos)

2015-11-04 Thread Jeff Newmiller
You assert "all the variables are defined and have values" but something is wrong that you are not showing us. Please read [1] and make your example reproducible. Also post in plain text (a setting in your email software) so we don't have to decipher a mangled version of what you write on

Re: [R] CRAN - Package "activity" - how do you convert time (hh:mm) into radians in R 3.1.2 on Windows 7 OS 64bit

2015-11-04 Thread Jeff Newmiller
Think astronomy. https://en.m.wikipedia.org/wiki/Hour_angle --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go...

Re: [R] Error in eval(expr, envir, enclos)

2015-11-04 Thread William Dunlap
No, I meant for you to (a) show an example that anyone could run and reproduce your problem (b) show the output you get when calling traceback() after you ran into the error. The lines and error messages you show now are different than the ones you showed before, but I see you are call

[R] Group Lasso for Proportional Odds Model

2015-11-04 Thread Nussbaum Madlene
Dear R experts, Having an ordinal response (3 levels) and factors in the predictors set, I would like to fit Proportional Odds Group Lasso. I found several packages that fit the group lasso: grplasso, grpreg, gglasso and ordPens. However, they vary penalty and implementation in several ways,

[R] Count of included observations in sp.correlogram

2015-11-04 Thread Robert U
Dear Rusers, I’m tryingto figure out what I think is a pretty simple thing for anyone who knows about correlograms.I’ve a regular grid (say 5*5 points) with some quantity associated to eachpoint (count data). I’m trying to verify whether this quantity is regularly /randomly or “clusterdly”

Re: [R] Extract entries from matrix

2015-11-04 Thread Jorge I Velez
Dear Dr. Steipe, Thank you for the code. It does exactly what I needed. Best regards, Jorge.- On Wed, Oct 28, 2015 at 1:18 PM, Boris Steipe wrote: > Your code does not produce the matrix in your image. > The first three rows contain all-zeros and the last row is

Re: [R] problem in R-code

2015-11-04 Thread thanoon younis
Dear R-Users After correct some errors in the code below i have only this error "Error in if (BZ[i] < 0.25) { : missing value where TRUE/FALSE needed" how can i solve this problem please? N<-200;P<-9 #Sample size BZ=matrix(NA, nrow=N, ncol=1) W=matrix(NA, nrow=N,

Re: [R] CRAN - Package "activity" - how do you convert time (hh:mm) into radians in R 3.1.2 on Windows 7 OS 64bit

2015-11-04 Thread Martin Maechler
> Bert Gunter > on Tue, 3 Nov 2015 14:54:46 -0800 writes: > ... should have said: with a 24 hour "clock". > -- Bert > Bert Gunter Hmm, thank you Bert. I know nothing about measuring time in radians, and for me, 'radians' are angles, and hence we

[R] R Consortium projects

2015-11-04 Thread Hadley Wickham
Hi all, I'm very pleased to announce that the Infrastructure Steering Committee (ISC) of the R consortium is calling for proposals: https://www.r-consortium.org/about/isc/proposals In brief: * We want to fund projects that help the R community, broadly construed. * Currently, we are mostly

Re: [R] CRAN - Package "activity" - how do you convert time (hh:mm) into radians in R 3.1.2 on Windows 7 OS 64bit

2015-11-04 Thread Jim Lemon
Ah, Martin, how is it possible that you have never encountered the 24 hour cuckoo clock? http://www.lutececreations.com/fiche-schneider_black_forest_cuckoo_clock__24_hour_cuckoo_clock_mt_1139_9++de+la+marque-4297.html Jim On Wed, Nov 4, 2015 at 7:02 PM, Martin Maechler

Re: [R] problem in R-code

2015-11-04 Thread Jim Lemon
Hi thanoon, Well, you have generated a one column matrix of missing values (NA) and then tried to use those values in a logical test... Jim On Wed, Nov 4, 2015 at 8:02 PM, thanoon younis wrote: > Dear R-Users > > After correct some errors in the code below i have