similar to: A Lattice Question

Hi, The data is x <- c(4,3,5); I need to translate it to y <- c(1,2,3,4,1,2,3,1,2,3,4,5); i.e. for each number in x, we need to generate 1:x and put it in y. The program need to evaluate this type of calculation for millions of times in simulation. Is there any elegant way to make this happen without iteration? Many thanks!! Best, Haoda [[alternative HTML version deleted]]

Dear all - Is there anyone know how to let VB or C# know where I install R automatically(i.e. auto detect R directory)? Many thanks! Best, Haoda

Hi, I am wondering if there is a simple way to plot MCF(mean cumulative function) for recurrent events in R? MCF (http://www.weibull.com/hotwire/issue57/relbasics57.htm) And do you have some other recommendation in visualize recurrent events? I did some research online. Some people use ggplot2 to plot MCF, but it requires some additional coding. http://user2010.org/slides/Shentu.pdf Your help

Hi - How can I add different notes in different panels? My data looks like ID Dose Visit Value 1 1 0 0.5 -6.5802e-02 2 1 0 1.0 2.4085e-01 3 1 0 1.5 -2.2907e-01 4 1 0 2.0 2.4074e-01 ... ... ... ... ... 270 45 30 3.0 -8.1316e-01 271 46 30 0.5 -2.0786e-01 272 46 30 1.0 -2.9336e-01 273 46 30 1.5 -5.5657e-01 274 46 30 2.0 -2.1659e-01

Hi - I am not R expert and I would appreciate your time if you can help me about my xyplot question. I would like to add text (p-value) in a 4 panels xyplot. I thought panel = function{} should work but I am not sure where I did it wrong. The error message from the following code is "Argument subscripts is missing with no default values" xyplot(GLG ~ PD | factor(TRT) , groups =

Hi there. I have now found that dist() does what I want in a specific case. But I have (out of many fields) in my Fragebogen.data file: Vpn Code Family Test 1 X1 m 45 2 X1 t 58 3 X2 m 44 4 X2 t 43 ... When I now do: Fbg <- read.table("Fragebogen.data", header=TRUE) dist(Fbg['Test']) I have the distances between everyone to everyone. Now I want to have the mothers (m in

Hi,I would like add a horizontal line at 126 (col=red) and a vertical line at 6.5 ( col= blue) in each panel .How should I use the panel.abline function in the following code I am using: ------------------ library(lattice) with(reg.dat.5,coplot(lbxglu~lbxgh|eth,panel=panel.smooth,xlab="ABC", ylab="FBG")) ---------------------------------------- Thanks a lot. Professor of

All - When I plot something like a<-rnorm(5) b<-rnorm(5) plot(a,b,col = "red") points(10,-10) The last point is missing because it is out of range of the first plot. I just try to switch from Matlab to R. In Matlab, it always can automatic adjust the xlim and ylim for such case. Is it possible auto adjust in R? Otherwise keep tracking xlim and ylim is really annoying.

Is there any way to control the format of the axis label numbers on a plot? More specifically, I have some plots that get axes with label numbers in exponential format, and I'd like to change that to non-exponential. Thanks!! -- M. Edward (Ed) Borasky, MS, MNLP, NST, FBG, PGS & PTA znmeb at borasky-research.net http://www.borasky-research.net

Anyone know how many simultaneous calls you can receive on a NuFone DID? -Mark -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20050310/e20d2eb8/attachment.htm

Is anyone here using Rspec and Ruby to test software written in other languages? I''m about to start an R language project, and may also be doing some work in another language (I''m not sure which one yet, but it''s one of the ones that''s part of the Axiom symbolic math package). I''d dearly love a test framework that can handle all three languages,

Hi, Is there anyone know if BUGS language allows the combination of variables as response such as Y[i] <- a*X1[i]+b*X2[i] Y[i] ~ dnorm(c,d) It seems doesn't work in my model. The problem is between two ######. The error message is > modelCheck("BayesBioMarker.BUGS") model is syntactically correct > modelData(paste("BUGS_data.txt",sep="")) data

My experience is very limited on Unix and Samba. Please forgive me if the answer is right in front of me. Problems accessing shares on my Samba 3.0.1pre1 running on Gentoo1.4 I'm stuck big time! I've been trying to figure out why I cannot access shares. The only share I can currently read & write to is: [OpenShare] comment = Temporary file space path = /tmp

Looking at http://rsync.samba.org/ftp/rsync/patches/fileflags.diff -- It looks like the changes from fileflags-fixes.diff patch were applied to the patch from http://rsync.samba.org/ftp/rsync/rsync-patches-3.0.0.tar.gz , but this entire chunk of the that original diff file was lost: > diff -up a/config.h.in b/config.h.in > --- a/config.h.in > +++ b/config.h.in > @@ -64,6 +64,9 @@

HI: I need ideas on how to make this code shorter (maybe with a second loop?). The code as it is works, but in this case I only have 14 samples, but it will become insane with more, so I need a way to make it more automatic. The problem is that the output from ts1, ts2, and so on is a vector with more than one value, so I do not know how to solve this. Thanks Prenewbie The code is the

Hi there, I'm rying to 'turn' an Schunk in an .Rnw file(Xemacs-21.4.13, ESS-5.2.8, R-2.1, miktex-2.4.1705). Has anyone got the isorot package to work with Sweave? JC example test.Rnw: \documentclass[a4paper]{article} \usepackage{Sweave} \usepackage{isorot} \rotdriver{dvips} \clockwise \title{Sweave Example 1} \author{apologies to Friedrich Leisch } \begin{document} \maketitle

Hello R Gurus: I would like to take a character string and split at the $ sign. I thought that strsplit would do it, but here are the results: > vv [1] "whine$ts1" > vv [1] "whine$ts1" > strsplit(vv,"$") [[1]] [1] "whine$ts1" Does anyone have any suggestions, please? Thanks, Edna Bell

Dear R-help, How could a cross-correlation plot be optimized such that the relationship between seasonal time-series can be studied? We are working with strong seasonal time-series and derived a cross-correlation plot to study the relationship between time-series. The seasonal variation however strongly influences the cross-correlation plot and the plot seems to be ?rather? symmetrical (max

How can I get a single attribute value of an object ? I jhave the tiemSeries object > ts1 Open 2003-10-09 02:00:00 1.27 2003-10-10 02:00:00 1.25 2003-10-13 02:00:00 1.27 2003-10-14 02:00:00 1.29 When I unclass ts1 I get: > unclass(ts1) list() attr(,"Data") Open 2003-10-09 02:00:00 1.27 2003-10-10 02:00:00 1.25 2003-10-13 02:00:00 1.27

HI, May be this helps: ?ts1<- ts(1:20) ?ts2<- ts(1:25) ts1[-(1:3)]<- ts1[-(1:3)]+ts2[1:17] ?as.numeric(ts1) # [1]? 1? 2? 3? 5? 7? 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 A.K. Hey everyone, I`m an absolut beginner in R and need some help for an exercise: I want to do ordinary calculations with 2 time series. The issue with this, that I want to use different elements of time