[R] multiple model selection for vglm models

Hi list, I'm working with vglm (family tobit) models and I need to perform multiple model selection. Something similar to step or stepAIC, is there anything alike? I found that neither step nor stepAIC work for objects of class vglm. Has anybody worked with them and can give a hint on how to

Re: [R] rgl.snapshot only captures a small portion what's visible in the RGL device window on CentOS 7

On 11/01/2016 5:39 PM, Wittner, Ben, Ph.D. wrote: The problem discussed below was fixed on my computer with much help from Tom Callaway, of RedHat. My computer has an NVIDIA graphics card and CentOS 7 comes with an open-source driver for NVIDIA graphic cards called nouveau. When I replaced the

Re: [R] rgl.snapshot only captures a small portion what's visible in the RGL device window on CentOS 7

The problem discussed below was fixed on my computer with much help from Tom Callaway, of RedHat. My computer has an NVIDIA graphics card and CentOS 7 comes with an open-source driver for NVIDIA graphic cards called nouveau. When I replaced the nouveau driver with a driver from NVIDIA, the

[R] Propensity score matching with MatchIt

Hello R experts, I am trying to do Propensity score matching for a medical data with two types of surgery. But somehow I am getting Summary of balance for all data and the matched data exactly similar resulting the Percent Balance Improvement as zero. > surgery.data<-read.csv(file.choose(), head

Re: [R] embedding expression into title in R plot

On 1/11/2016 2:08 PM, William Dunlap wrote: I tend to use bquote, as in x_label <- bquote(bold(species) ~ (italic(N1))) plot(1:10,main=bquote("This is the expression for" ~ .(x_label) * "!")) Thanks -- I thought I'd tried something very close to this in my various attempts, but it

Re: [R] [Q] It it possible to create the data frame only non-zero data column?

Hello Ulrik, Ashis and Petr Thank you for replying. I'll use dat[dat==0] <- NA. Thanks 2016年1月11日(月) 16:17 PIKAL Petr : > Hi > > as.matrix is rather dangerous, it converts all values to lowest mode > which, if there is text column, is character. > > And if I

Re: [R] embedding expression into title in R plot

David -- On 1/11/2016 1:01 PM, David Winsemius wrote: On Jan 11, 2016, at 7:59 AM, Evan Cooch wrote: Suppose I've specified that the xlab for a plot is expression(bold(species~(italic(N1 In other words, I want the axis label to be bold, italic 'species (N1)'

Re: [R] different coloured axis title labels for different axes

On 1/11/2016 1:54 PM, William Dunlap wrote: The following shows how to get different colors for most features of a scatterplot: plot(1:11,log(1:11),ann=FALSE,axes=FALSE,col="pink",pch=16) box(col="gray") title(xlab="X Axis Label", col.lab="light blue") title(ylab="Y Axis Label",

Re: [R] OpenAir package help!!

Hi Diana, As far as I can see, polarPlot (openair) does not specify the colors for the labels. One possibility is to redefine the default foreground color and see what happens: par(fg="red") but this will almost certainly change other elements in the plot. Unfortunately the panel.levelplot

Re: [R] Order of factor levels

I left out the example: > set.seed(1) > df <- data.frame(x1 = rpois(1000,4), x2 = rpois(1000,8)) > helper_fun <- function(x) { + cut(x, breaks = unique(quantile(x, seq(0, 1, 1/10), na.rm = TRUE)), + include.lowest = TRUE) + } > df2 <- data.frame(lapply(df, helper_fun)) >

Re: [R] Help request from Ph.D. Students

Dear NECMETTİN ALPAY KOÇAK, Re: > Dear All, > I am Ph.D. student in Econometrics. My thesis is about "Linear Filtering on a > Time Series" which R has already a nice package, namely "Filter". This > package is really helpful for my study. But, I really help from you to create > two filter

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

Dear R users, I would like to announce a new package that has just the appeared on CRAN, called "venn" version 1.0: http://cran.r-project.org/web/packages/venn/ (binaries will appear in one or two days) Although there are quite a few packages that draw Venn diagrams, there are a number of

Re: [R] [FORGED] Package "boot" and Bca for %

This question is not for homework... This question is about a real question about population structure with different morphotypes...   Thanks. - Mail original - De: "Rolf Turner" À: "msg contact" Envoyé: Lundi 11 Janvier 2016

Re: [R] How to suppress console output when using choose.dir()?

Thank you. Found uninstalling PC bloatware 'Acer Portal' rectified the problem. -Original Message- From: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de] Sent: Sunday, 10 January 2016 6:51 PM To: Rob Grant; r-help@r-project.org Subject: Re: [R] How to suppress console output when

Re: [R] rgl.snapshot only captures a small portion what's visible in the RGL device window on CentOS 7

Dear Duncan, Thanks very much for your reply. I guess I'll start with R-SIG-Fedora. I'll keep you cc'd unless you tell me you would rather I not. -Ben -Original Message- From: Duncan Murdoch [mailto:murdoch.dun...@gmail.com] Sent: Friday, January 08, 2016 8:13 PM To: Wittner, Ben,

[R] OpenAir package help!!

Dear all, I just wondered if anyone could help me on how to change the color of labels of a polar coordinate in a polar plot using the openair package of R. Kindest regards, Diana __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see

Re: [R] Help request from Ph.D. Students

Try DLM package : dlm: Bayesian and Likelihood Analysis of Dynamic Linear Models M.O On Mon, 2016-01-11 at 13:07 +0100, Franklin Bretschneider wrote: > Dear NECMETTİN ALPAY KOÇAK, > > Re: > > > Dear All, > > I am Ph.D. student in Econometrics. My thesis is about "Linear Filtering on > > a

Re: [R] different coloured axis title labels for different axes

How about this: plot(0,0,xlab='',ylab='') mtext(side=1,line=3,text='x axis',col=4) mtext(side=2,line=3,text='y axis',col=2) -Dan On Mon, Jan 11, 2016 at 7:44 AM, Evan Cooch wrote: > Consider a simple plot of X vs Y. There are elements on the plot that > represent X, or

Re: [R] embedding expression into title in R plot

> On Jan 11, 2016, at 7:59 AM, Evan Cooch wrote: > > Suppose I've specified that the xlab for a plot is > > expression(bold(species~(italic(N1 > > In other words, I want the axis label to be bold, italic 'species (N1)' > > Now, I want the title for the plot to be

Re: [R] different coloured axis title labels for different axes

> On Jan 11, 2016, at 7:44 AM, Evan Cooch wrote: > > Consider a simple plot of X vs Y. There are elements on the plot that > represent X, or Y, that are presented in different colours (say, blue > for X, red for Y). Rather than use a legend, I would like to have the

[R] embedding expression into title in R plot

Suppose I've specified that the xlab for a plot is expression(bold(species~(italic(N1 In other words, I want the axis label to be bold, italic 'species (N1)' Now, I want the title for the plot to be have this label embedded in the title. Say, 'This is the plot for Species (N1)'. For a

[R] different coloured axis title labels for different axes

Consider a simple plot of X vs Y. There are elements on the plot that represent X, or Y, that are presented in different colours (say, blue for X, red for Y). Rather than use a legend, I would like to have the title label for the X-axis be in blue, and the title label for the Y-axis be in

Re: [R] different coloured axis title labels for different axes

The following shows how to get different colors for most features of a scatterplot: plot(1:11,log(1:11),ann=FALSE,axes=FALSE,col="pink",pch=16) box(col="gray") title(xlab="X Axis Label", col.lab="light blue") title(ylab="Y Axis Label", col.lab="light green") axis(side=1, at=c(2,3,5,7,11),

Re: [R] embedding expression into title in R plot

I tend to use bquote, as in x_label <- bquote(bold(species) ~ (italic(N1))) plot(1:10,main=bquote("This is the expression for" ~ .(x_label) * "!")) Bill Dunlap TIBCO Software wdunlap tibco.com On Mon, Jan 11, 2016 at 7:59 AM, Evan Cooch wrote: > Suppose I've

[R] Order of factor levels

Dear list, What is a better way relative to the one below to keep the order of factor levels created from cut()? Notice, I'm simply pasting letters to levels before converting to character so to keep the desired order of levels. This is not very elegant... I'm converting to character so I can

Re: [R] Order of factor levels

Don't use vapply() here - use lapply() instead and then leave cut's output alone. vapply() will combine its outputs to create a character matrix and data.frame will pull apart the character matrix into its columns. Skipping the matrix intermediary solves lots of issues. Bill Dunlap TIBCO

Re: [R] Order of factor levels

Here's a solution with dplyr my_cut <- function(x){ breaks <- quantile(x, seq(0, 1, by = 0.1)) y <- cut(x, breaks = breaks, include.lowest = TRUE) levels(y) <- paste(head(letters, length(breaks) - 1), levels(y), sep = ": ") return(y) } library(dplyr) mutate_each(df, funs = funs(my_cut))