Re: [R] Does "coeftest" correctly use weights from "svydesign" in "svyglm" object?

hi, that setup is not correct. see examples in https://github.com/ajdamico/asdfree/tree/master/European%20Social%20Survey On Feb 8, 2017 11:54 PM, "André Grow" wrote: > Dear all, > > > > I am using data from the European Social Survey (ESS) and I would like to >

Re: [R] GLM and POST HOC test INTERPRETATION

Your questions are basically statistical and therefore OT here, although some kind soul may respond. I would strongly suggest that you consult with a local statistical expert, as you seem to be out of your depth statistically. Cheers, Bert Bert Gunter "The trouble with having an open mind is

[R] GLM and POST HOC test INTERPRETATION

Dear colleagues, I am analyzing a data set of 68 values (integers). In some treatments (exactly 6) the values are "zero". Because I record 0 in my measurement (or really a small value below zero) My experiment is designed in such a way that I record values for 6 treatments at 2 times.

[R] simplifying output of mapply

Hi I'm trying to use mapply to generate graphs but I cant figure out how to get the program to only save the graph to the array and not the other output. plt<-function (x,y) { p<-dfs[dfs$region==x,c("region",y)] colnames(p)<-c("v1","v2")

[R] Does "coeftest" correctly use weights from "svydesign" in "svyglm" object?

Dear all, I am using data from the European Social Survey (ESS) and I would like to calculate country-level cluster-robust standard errors for a regression model in R that includes country fixed effects and employs the design weights that come with the ESS. To correctly use the weights, I

Re: [R-es] Visualización de un gráfico guardado como variable

Hola Patricio, Prueba con bxp01 <- boxplot(var~itv,data) bxp(bxp01) Un saludo, Klaus. On 08-Feb-17 18:48, patricio fuenmayor wrote: Hola. En un análisis estoy guardando un gráfico en una variable, ya que necesito utilizar las estadísticas que este genera: bxp01 <- boxplot(var~itv,data),

Re: [R] Gaussian Filter

On Tue, 7 Feb 2017 19:51:21 +0200 catalin roibu wrote: The simplest way to get answers to such questions is run an internet search on your terms: "R Gaussian filter." Google reports 6,600,000 hits. The hits you are interested in are near the top. It wastes less time

Re: [R-es] Visualización de un gráfico guardado como variable

Hola, Tienes que llamar de nuevo a la función "boxplot()". Solamente se genera de nuevo el gráfico con "print" si utilizas funciones gráficas de los paquetes "ggplot" y de "lattice". Saludos, Carlos Ortega www.qualityexcellence.es El 8 de febrero de 2017, 18:48, patricio fuenmayor <

Re: [R] rms::latex.anova broken?

In recent versions of rms on CRAN there was a non-downward compatible change. To get latex output for summary, anova, and print on fit objects you leave off file="" (because we're usually using knitr anyway) and use options(prType='latex') anova(f) # LaTeX output You can use

Re: [R] How do I best create a R procedure from a R file?

Rough sketch: lookup.vec <- elencositi$nome.sito names(ookup.vec) <- elencositi$indirezzo.sito dati$FONTE <- lookup.vec(dati$FONTE) This, however, assumes that elencositi has all the values that con occur. > On 8 Feb 2017, at 10:27, Luca Meyer wrote: > > Hi, > > I am

Re: [R-es] Visualización de un gráfico guardado c omo variable

Hola Patricio, Creo que puede ser usando print(bxp01) Saludos, Guillermo > Hola. > En un análisis estoy guardando un gráfico en una variable, ya que necesito > utilizar las estadísticas que este genera: bxp01 <- boxplot(var~itv,data), > específicamente los outliers. > Luego necesito mostrar la

[R-es] Visualización de un gráfico guardado como variable

Hola. En un análisis estoy guardando un gráfico en una variable, ya que necesito utilizar las estadísticas que este genera: bxp01 <- boxplot(var~itv,data), específicamente los outliers. Luego necesito mostrar la gráfica. Pensaría que solo es llamando a bxp01 pero se me imprime el objeto. Existe

Re: [R] Gaussian Filter

Dear all! ​I know that! I need only the Gaussian filter, as for the rest is simple. Sent with Mailtrack On 8 February 2017 at 09:58, catalin roibu wrote: > Dear

Re: [R] Differential equations

It's been a while, but I think I have gotten through this sort of situation by splitting the integration into intervals, i.e., you run from t=0 to t=20 witn initial condition c(100,0), yielding a value c(y1,y2), then you run from 20 to 40 with initial condition c(y1+100, y2), etc. -pd On 08

Re: [R] Gaussian Filter

Always cc the list. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Feb 7, 2017 at 11:14 PM, catalin roibu wrote: > Dear

[R] calculating the cumulative distribution function for a Von Mises Mixture model?

Hello Is there a package for calculating the cumulative distribution function for a Von Mises Mixture model? I am using vonMF to calculate kappa, mu and the weighting parameter for Von Mises Mixtures; I now need to calculate the cumulative distribution function (CDF) for these mixtures. I

[R] Differential equations

Hi, I'm working on a system of 2 differential equations. My initial condition (t=0) is c(100,0) and i'm using lsoda function (from package deSolve) to solve it. My system reprensents the evoution of drug concentration in two compartments throug time. Problem is I would like to model a repeated

Re: [R] How do I best create a R procedure from a R file?

Hi one possibility is to use merge but it generates new data. Match can result in required data, but maybe there are also other options. First some toy data edat<-data.frame(letters, LETTERS) set.seed(111) dati<-data.frame(FFF=sample(letters[8:15], 35, replace=T))

Re: [R] pls package - validation

Bert Gunter writes: > However, if I understand correctly, using pls or anything else to try > to fit (some combination of) 501 variables to 16 data points -- and > then crossvalidate with 6 data points -- is utter nonsense. You just > have a fancy random number generator!

[R] fine tune changepoint

Dear all, I am using the package changepoint to identify a cut-off between positive and negative values for an assay. I followed the vignette of the package and the results are replicated, so I tried with a dataset of my choice. In the assay, the positive have a value above 0.8 and the negative

[R] How do I best create a R procedure from a R file?

Hi, I am working on the following file: > str(elencositi) 'data.frame':641 obs. of 2 variables: $ indirizzo.sito: chr "10ahora.com.ar" "abceconomia.co" "accmag.com" " actu.orange.fr" ... $ nome.sito : chr "10ahora" "ABC economia" "Acc Magazine" "Orange Actu" ... > head(elencositi)