similar to: Visualizing binary response data?

Hi All, I'm wondering if its possible to write degree in symbol. I would like y-label as "Temperature (degreeF)". where degree should be in symbols. Thanks in advance, #R Code library(lattice) data(barley) barchart(yield ~ variety | site, data = barley, groups = year, layout = c(1,6), ylab = "Temperature (degreeF)", scales = list(x =

Hi All, Can someone help me reshape following data: Var1 Var2 Val A X 1 A Y 2 A Z 3 B X 4 B Y 5 B Z 6 to some kind of matrix/tabular format (preferably as a matrix), may be like Var1 X Y Z A 1 2 3 B 4 5 6 Any help would be greatly appreciated, Kim [[alternative HTML version deleted]]

Hi All, I'm a SAS user but I'm very much interested in learning R. I use ODS system in SAS to make nice frequency tables. Is it possible to export the output of table() [in TABULAR FORM]? So, that I can use those directly for publications? Thank you. # R Code: library(datasets) Orange summary(Orange) # outputing this, not as it is, but in table for. is it possible? ~Kim [[alternative

Hi Kim and Others, Can anyone please help me on how can I edit default labels ( "Study 1", "Study 2",...) using forest() - same question as below. Which option I should try? It would be a great help. Thanks in advance, On Sun, Dec 13, 2009 at 7:14 PM, Kim Jung Hwa <kimhwamaillist@gmail.com>wrote: > Hi All, > > I'm using forest() from metafor package to

Hi All, I need a small help with following code: I'm trying to convert "dashed lines" to regular ones; and changing default "blue" border color to say "black"... but I'm doing it wrong and its not working. Can anyone help please. Thanks, Code: require(lattice) ex <- data.frame(x=1:10, y=rep(c("A","B"), 5)) bwplot(y~x, data=ex,

Hi All, I need a file which I can import to MS Word, I'm trying win.metafile(), but it does not seem to support multiple print commands at once (please see below). Is there an alternative to get plots file which can be used in powerpoint/word? # R code: # this does not work; but same thing works with pdf() library(lattice) win.metafile("test.wmf") p1 <- xyplot(decrease ~

Hi All, I'm fitting a Poisson regression. And I want to plot 95% Confidence Interval of Regression Estimates. After coming back to original scale (using following formula): exponential(estimate +/- 1.96*SE), at best I can get the output in the form of estimates, lower_limit, upper_limit values. As far I know forest() in metafor package needs input in the form of estimates and their

Hello! I would like to draw a graph like the following: http://www.optics.rochester.edu/workgroups/cml/opt307/spr04/pavel/plot_small.jpg Aim is to plot confidence intervals of treatments for X(=response1) and Y(=response2) axis simultaneously to visualize aggreement of confidence interval for two responses. Can anyone please provide me some direction to start with? Thanks! -- Kim.

Hello, I want to arrage the label according to my preference eg.. (va, vp, a, b, c) but don't know how to supress default ordering. Any suggestions? Please try the code below: n <- c("va", "vp", "a", "b", "c") p <- c(2, 2,1, 3,5) pm<- c(3,4,2,5,4) pn <- c(1,1,1,2,3) x<-data.frame(cbind(n,p,pm,pn)) library(ggplot2) qplot(x=n,

Hi All, I'm using glm() in R to perform Poisson regression, I'm wondering if its possible to get equivalent Type 1 / Type 3 Analysis (similar to one in PROC GENMOD). Thanks, Kim [[alternative HTML version deleted]]

Hello, I have a small statistics question, and as I'm quite new to statistics and R, I'm not sure if I'm doing things correctly. I am looking at two quantitative variables (x,y) that are correlated. When I divide the data set according to a categorical variable z, then x and y are more poorly correlated when z = A than when z = B (see attached figure). In fact x and y are two

Hi All, I'm trying to draw boxplots. I'm having a hard time to get "ticks labels" on multiple panels using 'alternating' option. # R Code: # May not be the best example, please just look into 'scales' option library(lattice) data(OrchardSprays) dta <- subset(OrchardSprays, OrchardSprays$rowpos %in% c(1,2,3)) # Original # This works fine, as you can see

Hi All, Can anyone please help me with getting a "single title" and "legend" for both the plots in the following R code. I'll eventually be using .wmf file. # R code: library(lattice) p1 <- xyplot(decrease ~ treatment, OrchardSprays, groups = rowpos, type = "a", main="Same title / legend", auto.key = list(space = "right",

Hi All, I'm trying a 3D-plot using wireframe() from Lattice package. Below is my code and sample data (read in "dta" object). I'm wondering if it is possible to make center grid (curved) line more prominent (or bold). The curved line I'm talking about is the center line drawn for a fixed value of log(B) [i.e. log(8)=0.90]. Any help would be highly appreciated. Thanks in

Hi All, I'm using 'tck' option to *reduce* the length of tick marks but it is not working, can anyone please tell me where I'm going wrong... require(graphics) require(grDevices) x <- seq(-10, 10, length= 30) y <- x f <- function(x,y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r } z <- outer(x, y, f) z[is.na(z)] <- 1 # 'bg' works but 'tck' is not showing

This is a question about removing the intercept in a binomial glm() model with categorical predictors. V&R (3rd Ed. Ch7) and Chambers & Hastie (1993) were very helpful but I wasn't sure I got all the answers. In a simplistic example suppose I want to explore how disability (3 levels, profound, severe, and mild) affects the dichotomized outcome. The glm1 model (see below) is

This email isn't asking for assistance, but I thought R-help readers would find it interesting. This week we offered a half-day introduction to R for researchers at Augsburg University. The response was astonishing. Although Augsburg has no medical faculty and no engineers, there was far too much demand, with interest from every faculty (barring theology, "for one small

R version 1.3.0 OS: SunOS 5.7, but I think the same problem occurs with Windows An incorrect formula seems to be used to calculate the pearson residuals for a generalized linear model with a binomial response. Here is a simple program which gives (a) the pearson residuals calculated directly, (b) the pearson residuals from glm, and (c) the deviance residuals from glm. The first and last

I have tried glm.nb in the MASS package, but many models (I have 250 models with different combinations of predictors for fish counts data) either fail to converge or even diverge. I'm attempting to use the negbin function in the AOD package, but am unsure what to use for the "random" term, which is supposed to provide a right hand formula for the overdispersion parameter.

Hi. I'm attempting to fit a logistic/binomial model so I can determine the influence of landscape on the probability that a box gets used by a bird. I've looked at a few sources (MASS text, Dalgaard, Fox and google) and the examples are almost always based on tabular predictor variables. My data, however are not. I'm not sure if that is the source of the problems or not because the