similar to: problem with length()

I'm new to Hmisc and trying to get the following to work, but if I un- comment the y-scale list (in order to get a log-scale for the hazard ratio), the error bars become strangely large. The dataframe is simply ODS output from TPHREG in SAS. Can someone point me towards what I'm sure is a naive error? Many thanks! tphreg1_parameterestimates1987$LCL <-

If I take the following simple data: YEAR MONTH DAY WEIGHT.KG 2003 10 6 1.2 2003 10 12 1.2 2003 10 16 1.3 and format the date data and plot it: dates <- strptime(paste(DAY,MONTH,YEAR),"%d%m%Y") plot(c(min(dates),max(dates)),c(0,max(WEIGHT.KG)), xlab="Date",ylab="Weight (kg)",type="n") lines(dates,WEIGHT.KG) points(dates,WEIGHT.KG) I find that the

I am trying to specify a legend title to be other than the variable name, but I find that the legend splits because scale_shape() takes effect but scale_colour() does not. Can someone spot my error? Here's some toy code that produces the problem on my system (R2.8.1, windows, today's CRAN package updates): library(ggplot2) a <- sample(letters[1:5],100,replace=T) b <-

Hi all, I have a question regarding the estimation of a two-sided truncation regression model using maximum likelihood estimation. I want to estimate a linear model where the dependent variable? is greater 0 and upper truncated at 1. Does anyone know a R-package that contains a routine that could do truncated regression for this situation? With many thanks for any hint Albrecht Dr. Albrecht

Problem: Domain logons temporary fail for a few minutes on all smb clients on the network. Reason: When using WINS support in samba-1.9.17p5 and the client request a resolv of a host-adress the sys_gethostbyname() system call in nmbd will block. When resolving a non-existant hostname or communicating with a slow nameserver the entire nmbd functionality will be unavailable causing this

Author: Darren Moffat <Darren.Moffat at Sun.COM> Repository: /hg/zfs-crypto/gate Latest revision: 957d30a3607ed9f3cbe490da5894d1e1b2104033 Total changesets: 28 Log message: resync onnv_105 partial for 6713916 Files: usr/src/Makefile.lint usr/src/Targetdirs usr/src/cmd/Makefile usr/src/cmd/Makefile.cmd usr/src/cmd/acctadm/Makefile usr/src/cmd/acctadm/acctadm.xcl

Release Announcements --------------------- This is the latest stable release of the Samba 4.18 release series. Changes since 4.18.3 -------------------- o? Douglas Bagnall <douglas.bagnall at catalyst.net.nz> ?? * BUG 15404: Backport --pidl-developer fixes. o? Samuel Cabrero <scabrero at samba.org> ?? * BUG 14030: Named crashes on DLZ zone update. o? Bj?rn Jacke <bj at

Release Announcements --------------------- This is the latest stable release of the Samba 4.18 release series. Changes since 4.18.3 -------------------- o? Douglas Bagnall <douglas.bagnall at catalyst.net.nz> ?? * BUG 15404: Backport --pidl-developer fixes. o? Samuel Cabrero <scabrero at samba.org> ?? * BUG 14030: Named crashes on DLZ zone update. o? Bj?rn Jacke <bj at

In lattice #toy data library(ggplot2) library(lattice) x <- rnorm(100) y <- rnorm(100) k <- sample(c("Weak","Strong"),100,replace=T) j <- sample(c("Tall","Short"),100,replace=T) w <- data.frame(x,y,j,k) xyplot(y~x|j+k,scales=list(y=list(relation="free"))) will give you a scale in each subplot with a range equal to the range of y

I'm able to do this as follows, but am wondering if anyone knows a simpler way which still avoids explicit loops? > (mystring <- letters[1:5]) [1] "a" "b" "c" "d" "e" > unlist(sapply(mystring[-length(mystring)], + function(x) paste(x,mystring[(grep(x,mystring)+1):length(mystring)],sep=""))) a1 a2 a3

How can I look up the aspect ratio of a plot, so I can use that to correctly adjust the angle of text which is supposed to be parallel to a line in the plot? The following example code works for a 1:1 aspect ratio, but puts the text at the wrong angle if the plot region is short and wide or tall and narrow. I can't find a par() component containing the plot aspect ratio. It will be for

Hi, We have the following environment for live-migration with non-shared stroage between two nodes, Host OS: RHEL 6.3 Kernel: 2.6.32-279.el6.x86_64 Qemu-kvm: 1.2.0 libvirt: 0.10.1 and use "virsh" to do the job as virsh -c 'qemu:///system' migrate --live --persistent --copy-storage-all <guest-name> qemu+ssh://<target-node>/system The

The order() function allows you to specify multiple vectors, which are used successively to break ties. If I want to use many vectors to break ties (say, 25 or more), that are columns of a matrix or elements of a list, does anyone know a shortcut to do this without passing 25 arguments to order()? -- Levi Waldron post-doctoral fellow Jurisica Lab, Ontario Cancer Institute Division of Signaling

For example, subtracting 1:2 from the rows of a two-column matrix: > t(apply(matrix(1:6,ncol=2),MARGIN=1,function(y) y - 1:2)) [,1] [,2] [1,] 0 2 [2,] 1 3 [3,] 2 4 > sweep(matrix(1:6,ncol=2),MARGIN=2,1:2,FUN="-") [,1] [,2] [1,] 0 2 [2,] 1 3 [3,] 2 4 Is there a logic to this difference, or is it just a quirk of the history of these

I'm trying to plot pspline'd explanatory variables from coxph() models as per Therneau and Grambsch (Modeling Survival Data: Extending the Cox Model). They have s-plus functions for this at: http://www.mayo.edu/hsr/people/therneau/book/sfunction/gamterms.s http://www.mayo.edu/hsr/people/therneau/book/sfunction/plotterm.s I'd like to make these plots in R, but they make use of non-R

Hello, I'm attempting to port OpenSSH to IBM's S/390 mainframe. Things have gone well but I have come a little unstuck with the internal PRNG. Although the commands in ssh_prng_cmds are being executed the select() seems to be returning 0 and therefore the cose is assuming that the forked process has timed out. This could be a difference in the way that select is implemented on OS/390.

Hi STATISTICA has a function called "Two-way joining" (see http://www.statsoft.com/TEXTBOOK/stcluan.html#twotwo) and the reference material states that this is based on the method as published by Hartigan (found this paper: http://www.jstor.org/pss/2284710 through wikipedia). What is the relationship (if any) between the "heatmap" function in R and this technique? Is there an

I've got some songs that I need OGG-ified as high quality as possible for my site, and I'm not sure how to go about doing it. I've checked out the tools available, but I'm wondering if there's anything else I should know about the process (above and beyond how to run the apps). Any advice in this area would be greatly appreciated. These suckers have got to be PERFECT...

I have used S-Plus for Windows for some time and have a large collection of personal programs. I have recently installed R-1.1.1 under Linux. I can read the .ssc files from S-Plus into emacs and save to R, but this is apt to be tedious. Any shortcuts to recommend? Will linux R read the _Data files from Windows S-Plus? Would Windows R serve as a useful intermediate transfer? -------------- next

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 1599 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20010102/cf0e43dd/attachment.bin