[R] help plotting glmm values

2015-08-06 Thread Luis Fernando García
Dear Fellows, I´m sorry if my question is very basic, but I'm still new in the field of R and the GLMM. I have made a following glmm, and I need to plot the predicted mean values from the model and the CI, with this barplot +/- the CI, using the prey as explanatory factor. Although I have been

Re: [R] function ndq()

2015-08-06 Thread Uwe Ligges
Use an internet search engine. The one I used told me that the book "Statistics and Data with R: An Applied Approach Through Examples" by Yosef Cohen, Jeremiah Y. Cohen contains a function definition of nqd in section 1.11.4. Best, Uwe Ligges On 06.08.2015 23:31, kle...@sxmail.de wrote: H

Re: [R] R error

2015-08-06 Thread andymac
Hi, I had a similar problem running the 32-bit version of 3.2.0 on Windows 8.1. For me the problem was user permissions. I re-ran the R program as an administrator rather than a user and it now installs packages just fine. Cheers, Andy -- View this message in context: http://r.789695.n4.nab

[R] Potencial bug in R.adjust ("holm" method)

2015-08-06 Thread Iker Vaquero Alba
   Hello all,    I am doing some Bonferroni correction analyses with R.adjust function. I have a spreadsheet with 24 columns, each with 5 values. When I use the "holm" method, it gives me adjusted figures for all the original values except from the ones in the 4th row of each column. I mean, the

[R] APARCH models

2015-08-06 Thread SophieO
Dear all, I need to estimate the coefficients of an APARCH that fits some dataset. I have found to functions : * garchFit(formula=~aparch(1,1),...) from the fGarch package * GSgarch.FitAIC(mMAX = 0, nMAX = 0, pMAX = 1, qMAX = 1, ...) from the GEVStableGarch package The problem is that they gave

[R] adding a node in an existing xml file

2015-08-06 Thread hs
I have an xml file like this: http://www.opengis.net/kml/2.2";> sql_statement 2013/11/16 00:02:16 2013/11/16 00:02:16 70.3 65562 3.611564270451177,39.118047789677419 I'm able to read the x

Re: [R] Lattice: raw data and prediction of a non linear fitted function

2015-08-06 Thread François Collin
Hi, Thank you for your answer. However, I think you misunderstood my problem. The data and fit I provided are just here to illustrate 2 points: - I have few raw data - and a fitted model I want to represent. It is a simulated dataset, no more, which match to my requirement: a non linear response,

[R] Problem with survfit when model has interaction with strata

2015-08-06 Thread Sebastian Pölsterl
Dear all, I encountered a problem when running survfit using a coxph model that contains an interaction with a strata variable (see the attached example). I hope this is the right place to report problems. > source("~/survfit-example.R") Error in scale.default(x2, center = xcenter, scale = FA

[R] function ndq()

2015-08-06 Thread klerer
Hello everyone!I had to notice that it is better to be curious about packages and functions contained therein, as I (yet new to R) recognised that there is a number of 'convenient' R packages simplifying and possibilitating a lot. At the moment, I have 'retrofitted' (install packages not deliver

[R] Piecewise (segmented) linear regression with center section slope constraint

2015-08-06 Thread Drew Morrison
Hi, I'm working on a way to predict the electricity consumption of electrically heated buildings as a function of outdoor air temperature. I've identified a three-segment linear model as a candidate for a good fit, with the slope of the center section constrained to zero. I'm working with the segm

Re: [R] testing whether two character vectors contain (the same) items in the same order

2015-08-06 Thread Bert Gunter
Boris: You may be right, but it seems like esp to me based on the op's non-description of likelihood of coming from the same noisy process. My response would be: seek local statistical help, as your replies indicate a good deal of statistical confusion. Cheers, Bert On Thursday, August 6, 2015

Re: [R] Lattice: raw data and prediction of a non linear fitted function

2015-08-06 Thread Kevin Wright
I find this type of lattice graph to be much easier to create via the use of layers in the latticeExtra package. See the "carmer.density" dataset example in the "agridat" package, which shows the use of layers. It doesn't have "groups", but it is not too hard to add a groups argument to a layer a

Re: [R] testing whether two character vectors contain (the same) items in the same order

2015-08-06 Thread Boris Steipe
You are looking for what is known as the "Cayley distance" between vectors - an edit distance that allows only transpositions. RSeek mentions PerMallows (https://cran.r-project.org/web/packages/PerMallows/PerMallows.pdf) and Rankluster (https://cran.r-project.org/web/packages/Rankcluster/Rankcl

Re: [R] Knit R and Book Publishing

2015-08-06 Thread boB Rudis
https://github.com/hadley/adv-r is how it was done. On Thu, Aug 6, 2015 at 8:33 AM, Bert Gunter wrote: > I would have thought that the first place to look was R Studio support > site. You will find a lot of (Imo well done) docs there as well as links to > Hadley's and Yihui's books and online doc

Re: [R] Release schedule (was (no subject) )

2015-08-06 Thread peter dalgaard
> On 06 Aug 2015, at 18:08 , peter dalgaard wrote: > >> >> Except when it doesn't: "It is intended to have a final patch release >> of the previous version shortly before the next major release." > > Argh. Needs fixing... :-p > Done. (Of course, by the canonical definition of "major", Maci

Re: [R] Release schedule (was (no subject) )

2015-08-06 Thread peter dalgaard
> On 05 Aug 2015, at 22:47 , Duncan Murdoch wrote: > > On 05/08/2015 4:36 PM, peter dalgaard wrote: >> >>> On 05 Aug 2015, at 20:32 , Duncan Murdoch wrote: >>> >>> On 05/08/2015 2:15 PM, Jeff Newmiller wrote: So 3.1.3 to 3.2.0 was a major release? >>> >>> Yes. We do have the oddity (se

Re: [R] (no subject)

2015-08-06 Thread Adams, Jean
A quick internet search for R version history yielded this page https://cran.r-project.org/bin/windows/base/old/ Jean On Wed, Aug 5, 2015 at 4:55 AM, Djossè Parfait wrote: > Good morning, > > I would like to know how often per year is a new full version release of R. > > > > Thanks >

Re: [R] testing whether two character vectors contain (the same) items in the same order

2015-08-06 Thread Federico Calboli
> On 6 Aug 2015, at 15:40, Bert Gunter wrote: > > Define "goodness of match" . For exact matches, see ?"==" , all.equal, etc. Fair point. I would define it as a number that tells me how likely it is that the same (noisy) process produced both lists. BW F > > Bert > > On Thursday, Aug

Re: [R] testing whether two character vectors contain (the same) items in the same order

2015-08-06 Thread Bert Gunter
Define "goodness of match" . For exact matches, see ?"==" , all.equal, etc. Bert On Thursday, August 6, 2015, Federico Calboli wrote: > Hi All, > > let’s assume I have a vector of letters drawn only once from the alphabet: > > x = sample(letters, 15, replace = F) > x > [1] "z" "t" "g" "l" "u"

Re: [R] Knit R and Book Publishing

2015-08-06 Thread Bert Gunter
I would have thought that the first place to look was R Studio support site. You will find a lot of (Imo well done) docs there as well as links to Hadley's and Yihui's books and online docs. Bert On Wednesday, August 5, 2015, Glenn Schultz wrote: > Hello All, > > I have finished my first R pack

[R] testing whether two character vectors contain (the same) items in the same order

2015-08-06 Thread Federico Calboli
Hi All, let’s assume I have a vector of letters drawn only once from the alphabet: x = sample(letters, 15, replace = F) x [1] "z" "t" "g" "l" "u" "d" "w" "x" "a" "q" "k" "j" "f" "n" “v" y = x[c(1:7,9:8, 10:12, 14, 15, 13)] I would now like to test how good a match y is for x. Obviously I can

[R] Supervised Learning for Text Classification

2015-08-06 Thread Anshuk Pal Chaudhuri
Hi All, The current process which I am doing: 1. Reading Unstructured data, Understand what text means what. e.g. a phrase like "JOHN SMITH" this is a customer name. or a phrase like "CLAIM DURATION" is a reason type field. Currently, doing this manually (no issue over here). 2.