similar to: step function

Dear members, I would like a help for extracting the values from a step function (stepfun). >From help(stepfun) we have the following example: Y0<-c(1.,2.,4.,3.) y0<-c(1.,2.,3.,4.) sfun<-stepfun(1:3,y0,f=0) plot(sfun) Now, suppose instead I was given the object (*sfun*, say) from which I wanted to extract the values generated by the function *stepfun*. More precisely, I want to

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'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, I'm dealing with binary response data for the first time, and I'm confused about what kind of graphics I could explore in order to pick relevant predictors and their relation with response variable. I have 8-10 continuous predictors and 4-5 categorical predictors. Can anyone suggest what kind of graphics I can explore to see how predictors behave w.r.t. response variable... Any

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 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]]

Dear R developer, I am not quite sure, if I should post my concern as a wish to r-bugs at r-project.org. Thus, as recommended, I first send an email to you. My request is the following: I would appreciate, if it was possible to obtain a plot of a 'stepfun' with a strict interpretation of xlim. What I mean: sf <- stepfun(1:4, 1:5) plot(sf, xlim=c(0,10)) does not bound the function to

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

i've noticed a strange problem when plotting a stepfun. according to the documentation, the xlim parameter should bound the range of the function being plotted, and is returned as the extreme two values (i.e. first and last) in the vector t from the plot.stepfun call. instead, it plots beyond the desired range (although the limits are preserved for the viewing space). to reproduce: foo

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,

I am using Rprof() to help find ways to improve performance. I found a function whose total seconds and self seconds were large. I replaced it with something else. The something else had a small number of total seconds and self seconds. But the total time did not decrease. I don't understand how that could be, and would appreciate any suggestions. Thanks -Don Details, unfortunately

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]]

I believe the posting guide suggests that you not ask the question in the subject... plot(..., type="s") or ?stepfun Andy > From: john miller

On Tue, 30 Dec 2008, m.u.r. wrote in thread [R] plot.stepfun xlim: > foo <- stepfun(0.5, c(1, 0)); On Wed, Dec 31, 2008 at 4:18 AM, Prof Brian Ripley <ripley at stats.ox.ac.uk> replied: > Why are you adding two blank commands via the semicolons? The R parser (2.8.0 Windows) does not seem to have the concept of blank commands. After all: (using <NL> to represent Newlines

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,

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 list:) I have lots of values which I would like to get a histogram outline out of. An example of what I am talking about: testdata = runif(100) bbb = seq(0,1, by = 0.01) hist(testdata, breaks = bbb) I would like to get the outline of the resulting histogram. Now, I think that I can do this using the stepfun function. However, I am uncertain of how to get to the data the stepfun function

Hello, My specific system is > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 5.1 year 2002 month 06 day 17 language R Each of the hmisc and stepfun packages has a function

Hello, if I plot a horizontal line, e.g., plot(c(1,2),c(1,1),xlim=c(0,3),lwd=2,type="l") or plot(c(1,2),c(1,1),xlim=c(0,3),lwd=4,type="l") then the left end (1st example) or both ends (2nd example) of the lines are not rectangular but slanted on the graphical display (screen). That behavour first occurred when I was trying to plot a stepfun, e.g., y <-