similar to: Extracting and using fitted values and residuals with missing data

Hi Hadley and everyone, here's a patch for ggplot2 that fixes the behavior of opts(legend.position={left,top,bottom}). If you try the following code in an unmodified ggplot2 options(warn = -1) suppressPackageStartupMessages(library("ggplot2")) data <- data.frame( x = c(1, 2, 3, 4, 5, 6), y = c(2, 3, 4, 3, 4, 5), colour = c(TRUE, TRUE, TRUE, FALSE, FALSE, FALSE))

Hi all: following the previous discussion, it looks like plot(x) with y=NULL still does not work correctly. If one tries for example plot(1:5) it works, but already for plot(runif(100)) it does not. I posted the proposed correction for plot.POSIXct and plot.POSIXlt before. Please voice your opinions whether the following fix for plot.default could be reasonable? I include the full function and

I don't usually do much with graphs in R, and this is my first time adding a line of best fit. Hopefully this is an easy problem to solve. I'm looking at a variable called soloKills along the range 5:28. Here are all my commands, in script form: range=5:28 graph=soloKills title="Solo kill/death data" xlabel="Number of deaths/1 game" ylabel="Mean number of kills/1

Hi all, As long as I've used R, add=TRUE hasn't worked in contexts like this: f <- function(x) x^2 X <- seq(0, 1, by=1/4) plot(f, col="blue") plot(X, f(X), col="red", type="l", add=TRUE) I attached a fix for version 2.6.2. Cheers, Andrew -------------- next part -------------- diff --git a/src/library/graphics/R/plot.R

This is a feature request for mosaicplot, including the suggestion for the code changes (mosaicplot.R) and docs (mosaicplot.Rd) by Wolfram Fischer and myself, diff'ed to 1.6.0 beta (2002-09-14): Uwe Ligges .../src/library/base/R/mosaicplot.R: ==================================================== 11a12,13 > ### Changes by W. Fischer and U. Ligges: > ## - Deparsing x in for main

Hi! How to find out which plot function is used when i call plot(hist.default(1:10,plot=F)) and all works fine ? The reason why I would like to know it is that after loading some self written R functions > plot(hist.default(1:10,plot=F)) Error in xy.coords(x, y, xlabel, ylabel, log) : x and y lengths differ > traceback() 5: stop("x and y lengths differ") 4: xy.coords(x, y,

Hi, Here are a few errors I found as well as a few comments. 1) In the man page of par: lty: The line type. Line types can either be specified as an integer (0=blank, 1=solid, 2=dashed, 3=dot- ted, 4=dotdash, 5=longdash, 6=twodash) or as one of the character strings `"blank"', `"solid"',

I have a need to build a time series and there are a couple of aspects about the time series object that are confusing me. First it seems that ts.union is not doing what I would expect. For example: x0 <- rep(0,10) x1 <- rep(1,10) xt0 <- ts(x0, frequency=10) xt1 <- ts(x1, frequency=10) st2 <- ts.union(xt0, xt1) > xt2 Time Series: Start = c(1, 1) End = c(1, 10) Frequency = 10

I am currently trying to use R to construct a regression model to explain output based on temperature. I have combined my output and temp data into a notepad file. there is no problem with loading the data into R. > data.df output temp 1 850 17 2 849 17 3 905 17 4 925 17 5 1043 19 6 1104 20 7 1097 18 8 979 19 9 926 18 10 1133 18 ~~ ~~ 240 1124

Hi, I am Using R.2.9.2 and I am trying to make nicer some graphics, which are constructed with plotmer.fnc. In particular, I have two questions: 1. Is is a way to decide in which order we want the labels on the x-axis. For now, the reference value label is the first one, and the other labels appear in the alphabetical order. I would like them to be sorted in the ascending way. Here is my

Hi, I am new to R. I am trying to figure out how to graph expression data from the TCGA database. If I understand correctly the expression data I have downloaded is from a microarray using the AgilentG4502A. I've had trouble reading into R in the level I, level II, and the gene expression analysis data using >dat<-read.table("C:\\file.txt", header=T, row.names=1) for

termplot() does not carry subsetting over to carriers that are in the environment but not in the data frame. This generates a "subscript out of bounds" error. > data(ToothGrowth) > logdose <- log(ToothGrowth$dose) > tooth.lm <- lm(len ~ logdose, data=ToothGrowth) > termplot(tooth.lm) ## Works fine > toothVC2.lm <- lm(len ~ poly(dose,2),

Good day all, I'm trying to plot a continuous line plot, which is followed by two forecast points eg. one forecast point is 12 months out, and another 24 months out from the last date of the line plot. In my attempts so far, the second plot (the forecast points) is scaled against a new axis scale, thus the two plots are not directly comparable (I need the forecast points to be scaled

I'm having trouble getting fun="cloglog" to work with plot on a survfit object. Here are the data I used for the commands that follow. days status 2 0 2 0 5 1 9 0 14 1 16 0 16 0 17 0 29 1 30 0 37 1 37 0 39 1 44 0 44 0 58 0 60 1 67 1 68 1 82 1 82 1 86 0 86 0 89 1 93 0 97 1 100 0 100 0 100 0 > library(survival) Loading required package: splines > eg1.km <-

HI, sorry but i don't understand how to make a function with as.function() formula<-"2+3*x" formu<-as.symbol(formula) > formu 2+3*x formul<-as.function(alist(x=,formu)) curve(formul,1,5,col="blue") Error in xy.coords(x, y, xlabel, ylabel, log) : x and y lengths differ > typeof(formul) [1] "closure" and not plot the curve function, Why?

Hi i want do a line graph. My y axis contains numeric values. My x axis contains non numeric statements. This is what i want the graph to look like. When i try to plot this graph on R it comes up with the following error message: "Error in plot.window(...) : need finite 'xlim' values In addition: Warning messages: 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by

Dear all, I am trying to plot a matrix I have? as an image str(matrixToPlot) ?num [1:21, 1:66] 0 0 0 0 0 0 0 0 0 0 . ?that contains only 0s and 1s, where the xlabel will be Labeled as str(xLabel) ?num [1:66] 1e+09 1e+09 1e+09 1e+09 1e+09 ... and the yLabels will be labeled as str(yLabel) ?num [1:21] -88 -87 -86 -85 -84 -83 -82 -81 -80 -79 ... I have found on the internet that I can do

A table without names displays like a vector: > unname(table(2:3)) [1] 1 1 1 and preserves the table class (as with unname in general): > dput(unname(table(2:3))) structure(c(1L, 1L), .Dim = 2L, class = "table") Does that make sense? R is not consistent in its treatment of such unname'd tables: In plot, they are considered erroneous input: >

Hello, I am hoping for some advice regarding warning/error messages I received when running a Cox regression # message 1 - obtained while creating a plot of residuals > plot (NV.zph, main = "groupNUSM - UNFIT", var= 'groupNUSM') Warning messages: 1: In approx(xx, xtime, seq(min(xx), max(xx), length.out = 17)[2 * : collapsing to unique 'x' values 2: In

I am trying to plot survival curves using the following code as an example: >rs1799964.coxph<-(coxph(Surv(sassurvmonths,status)~age+stage+rs1799964_TNFA,method="efron")) >plot(rs1799964.coxph,lyt=c(1,3),xlab="Survival in Months",ylab="Proportion Surviving") I am gettingthe following error message: >Error in xy.coords(x, y, xlabel, ylabel, log) :