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 and

Re: [R-es] consulta graficas para GLM

2015-08-06 Thread Jose Luis Cañadas Reche
Hola. Prueba con el en�simo paquete de R , effects, algo asi. library(effects) efectos - allEffects(glm6) plot(efectos) Saludos El 06/08/15 a las 09:09, Luis Fernando Garc�a escribi�: Estimados amigos y expertos del R, Les escribo para hacerles una pregunta que parece un poco sencilla

Re: [R-es] consulta graficas para GLM

2015-08-06 Thread Jose Luis Cañadas Reche
Supongo que si quieres hacer el gr�fico 1, puedes hacer esto. efectos - Effect(c(Presa, Exposici�n1), glm6 ) plot(efectos) El 06/08/15 a las 09:09, Luis Fernando Garc�a escribi�: Estimados amigos y expertos del R, Les escribo para hacerles una pregunta que parece un poco sencilla pero me

[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 transform the letters in

[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

[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

[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

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:

[R] adding a node in an existing xml file

2015-08-06 Thread hs
I have an xml file like this: ?xml version=1.0 encoding=utf-8 ? kml xmlns=http://www.opengis.net/kml/2.2; Document id=root_doc Schema name=sql_statement id=sql_statement SimpleField name=Name type=string/SimpleField SimpleField name=Description type=string/SimpleField /Schema

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

[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 =

[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

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,

Re: [R-es] Duda interpolación (package ' gstat ')

2015-08-06 Thread rubenfcasal
Hola de nuevo, Un par de detalles técnicos... Antes de nada comentar que el paquete gstat no es computacionalmente muy eficiente calculando las predicciones kriging (calcula el estimador mcg de la tendencia utilizando la expresión explícita, etc...), entre otras cosas requiere la

[R-es] consulta graficas para GLM

2015-08-06 Thread Luis Fernando García
Estimados amigos y expertos del R, Les escribo para hacerles una pregunta que parece un poco sencilla pero me ha costado mucho. Estoy tratando de graficar los datos correspondientes a tasas de consumo de algunos organismos cuando están expuestos o no expuestos a un químico sobre tres tipos de

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 glennmschu...@me.com wrote: Hello All, I have

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 federico.calb...@helsinki.fi 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

Re: [R-es] Duda interpolación (package ' gstat ')

2015-08-06 Thread rubenfcasal
Hola Marcos, Parece que el problema es con el ajuste del variograma (sale plano?), sin más información no se exactamente que puede estar pasando... Si me envías el código completo y los datos lo miro con más detalle (e incluso te doy una alternativa no paramétrica con el paquete

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 bgunter.4...@gmail.com 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

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

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 bgunter.4...@gmail.com 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

[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.

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

2015-08-06 Thread peter dalgaard
On 05 Aug 2015, at 22:47 , Duncan Murdoch murdoch.dun...@gmail.com wrote: On 05/08/2015 4:36 PM, peter dalgaard wrote: On 05 Aug 2015, at 20:32 , Duncan Murdoch murdoch.dun...@gmail.com wrote: On 05/08/2015 2:15 PM, Jeff Newmiller wrote: So 3.1.3 to 3.2.0 was a major release? Yes.

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

2015-08-06 Thread peter dalgaard
On 06 Aug 2015, at 18:08 , peter dalgaard pda...@gmail.com 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,

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 djosseparf...@gmail.com wrote: Good morning, I would like to know how often per year is a new full version release of

Re: [R-es] Duda interpolación (package ' gstat ')

2015-08-06 Thread Marcos Bermejo
Sale plano sí. Ya se que sin tener los datos y el código es un poco difícil, pero es que mis datos ocupan mucho, es imposible. Seguiré mirando por internet. Muchas gracias Rubén. Un saludo, To: r-help-es@r-project.org From: rubenfca...@gmail.com Date: Thu, 6 Aug 2015 14:21:47 +0200