Displaying 20 results from an estimated 6000 matches similar to: "global object in user defined function"
2006 Sep 07
3
graphics - joining repeated measures with a line
I would like to join repeated measures for patients across two visits using
a line. The program below uses symbols to represent each patient. Basically,
I would like to join each pair of symbols.
library(lattice)
patient <- c(1,2,3,4,5,6,7,8,9,1,2,3,4,5,6,7,8,9)
var <-
c(826,119,168,90,572,323,122,10,42,900,250,180,120,650,400,130,12,33)
visit <- c(1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2)
2007 May 30
2
control axis
I have an outlier that I would still like to display, but would prefer to
shorten the axis. For example, display 0% - 40%, and 90% - 100%. Is this
possible? I am using an xyplot.
Thanks
Murray
--
Murray Pung
Statistician, Datapharm Australia Pty Ltd
0404 273 283
[[alternative HTML version deleted]]
2006 Sep 11
1
graphics: y limit on xyplot
I would like to set the y axis limit of an xyplot using the object 'ylimit',
but receive this error:
[1] 990
Error in extend.limits(limitlist[[i]], axs = axs) :
improper length of lim
I get the same error if I use ylim.
library(lattice)
trellis.device(col = FALSE, theme = lattice.getOption("col.whitebg"))
name <- "Variable name"
symbols <-
2007 Jul 23
1
maths characters in labels & ylab padding
I have checked out the help files, but cannot find details on how to use
maths characters in ylab. Instead of m^2, I would like the 2 in superscript,
if possible. I would also like to place more padding on the label so that
the label is not obscured by the horizontal numbers.
y <- 1:10
x <- rnorm(10,50000,2000)
plot(x ~ y,
ylab = 'Y Label (m^2)',
las = 1,
type =
2006 Oct 26
3
Read.csv
Are there ways to load a csv file without row.names by read.csv? Thanks.
[[alternative HTML version deleted]]
2006 Oct 25
1
density plot text
Is there any way of adding text to a density plot? I have had a go using the
text() function but I think the error is because this function doesn't work
with densityplot().
Alternatively, I understand I can achieve pretty much the same result if I
plot a density kernel estimate using plot() (which allows text()), but I do
prefer densityplot().
Also, is it possible to specify the dimensions
2006 Nov 17
2
s.e. on interaction plots
Is it possible to add standard error bars to the means on interaction plots?
Thanks
Murray
--
Murray Pung
Statistician, Datapharm Australia Pty Ltd
0404 273 283
[[alternative HTML version deleted]]
2007 May 31
1
cox goodness of fit
Is there an implementation of the Cox-Snell residuals / Nelson-Aalen plot
for goodness of fit?
Or otherwise is there an appropriate Goodness of Fit diagnostic?
Thanks
Murray
--
Murray Pung
Statistician, Datapharm Australia Pty Ltd
0404 273 283
[[alternative HTML version deleted]]
2007 May 15
1
differentiate groups on barplot
To differentiate between groups on the barplot, I guessed that col =
colr[test$group] would have worked. How can I do this?
Many Thanks
Murray
test <-
structure(list(patient = 1:20, score = c(100, 95, 80, 75,
64, 43, 42, 40, 37, 35, 30, 29, 27, 26, 23, 22, 19,
18, 17, 16), group = c(1, 0, 1, 0, 1, 0, 1, 0, 1,
0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0)), .Names = c("patient",
2007 May 16
2
log rank test p value
How can I get the Log - Rank p value to be output?
The chi square value can be output, so I was thinking if I can also have the
degrees of freedom output I could generate the p value, but can't see how to
find df either.
> (survtest <- survdiff(Surv(time, cens) ~ group, data = surv,rho=0))
Call:
survdiff(formula = Surv(time, cens) ~ group, data = surv, rho = 0)
N Observed
2007 May 30
3
sizing and saving graphics in R
Dear R wizards,
I am seeking advice on graphics in R. Specifically, how to manipulate
the size and save a plot I have produced using the LDheatmap library.
I confess I am relatively new to graphics in R, but I would greatly
appreciate any suggestions you may have.
LDheatmap produces a coloured triangular matrix of pairwise
associations between 600 genetic markers in my dataset.
2006 Nov 27
1
x axis on sciplot
Is it possible to format the x axis, so that days are spaced
proportionately?
i.e. a larger space between 50 & 100, than 0 and 5?
Have not had any luck with axis.POSIXct().
Many thanks
Murray
try <-
structure(list(visit = structure(as.integer(c(1, 2, 3, 4, 5,
1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1,
2, 3, 4, 5)), .Label = c("Screening", "0",
2007 Jun 17
1
error bars on survival curve
I am using plot(survfit(Surv(time,status) ~...) and would like to add
error bars rather than the confidence intervals. Am I able to do this
at specified times? e.g. when time = 20 & 40.
leukemia.surv <- survfit(Surv(time, status) ~ x, data = aml)
plot(leukemia.surv, lty = 2:3,xlim = c(0,50))
#can i add error bars at times 20 & 40?
legend(100, .9, c("Maintenance", "No
2009 Jul 27
2
Superstring in text()
I'd like to paste a superstring with a number in an object.
Thanks for any help.
Murray
mycor <- cor(1:10,1:10)
plot(1:10,1:10)
text(8,2,paste(expression(R^2)," = ",mycor))
[[alternative HTML version deleted]]
2005 Nov 23
8
getting started, reading listing and saving data
Dear List
I am new to R and to the list and will try best as I can be clear and
concise. My apologies if anything I write contravenes the posting code
on this list. I would also like to say I have run through most of the
material on the R website before writing this email however, I am
stuck.
Here is what I want to do and what I have done
1. Read a comma seperated text file into R
I have used
2012 Nov 22
1
Efficiently creating/defining new variables transformations
I would like to add an extension to the current name of a variable to create
a new variable that is its sqrt transform. Each piece of the equation below
works independently, but the left side definition fails on run. I also tried
creating the variable name first, but ended up with an object that
toString() did not fix. Better ideas?
Example 1
2005 Apr 21
2
Deciles and R
Hi everyone,
I'm a new R user (if this is a really basic question, please do excuse
me...) and I'm having some questions regarding a deciles problem.
I have a variable which I need to categorize according to its deciles (X).
However, this categorization should be made into another variable (call it
NewVar).
Ex. for the quartiles case (just for the sake of exposition, since I need
2007 Feb 01
2
Losing factor levels when moving variables from one context to another
Hi, there
I'm currently trying to figure out how to keep my "factor" levels for a
variable when moving it from one data frame or matrix to another.
Example below:
vec1<-(rep("10",5))
vec2<-(rep("30",5))
vec3<-(rep("80",5))
vecs<-c(vec1, vec2, vec3)
resp<-rnorm(2,15)
dat<-as.data.frame(cbind(resp, vecs))
2004 Jun 05
1
Crash in OSX (PR#6940)
Full_Name: Murray Pung
Version: 1.9.0
OS: OSX Mac
Submission from: (NULL) (134.148.20.33)
Date/Time: 2004-06-05 12:32:30 +1000
OS Version: 10.3.4 (Build 7H63)
Report Version: 2
Command: R.bin
Path: /Library/Frameworks/R.framework/Resources/bin/R.bin
Version: 1.9.0 (R 1.9.0)
PID: 358
Thread: 0
Exception: EXC_BAD_ACCESS (0x0001)
Codes: KERN_PROTECTION_FAILURE (0x0002) at
2006 May 24
5
Joining variables
Hello,
If I have two variables that are factors or characters and I want to
create a new variable that is the combination of both what function can
I use to accomplish this?
Ex.
Var1 Var2
SA100055113 19851113
And I want
NewVar
SA10005511319851113
Thanks in advance.
Cameron Guenther, Ph.D.
Associate Research Scientist
FWC/FWRI, Marine Fisheries Research
100 8th Avenue S.E.
St.