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

After I ran my program through Rubyscript2exe, the Random.exe file that it created is causing the following error: D:\Ruby projects>random.exe C:\DOCUME~1\JEDELS~1\eee\eee.Random.exe.2\app\Random.rb:59: undefined method ''ru n'' for Wxruby2::App (NoMethodError) from C:\DOCUME~1\JEDELS~1\eee\eee.Random.exe.2\bootstrap.rb:77:in `load'' from

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

Hi friends, I have following data and would like to plot this with barchart() availble with lattice package. RsID Freqs Genotype AAA 63.636 1/1 AAA 32.727 1/2 AAA 3.636 2/2 BBB 85.965 2/2 BBB 14.035 2/1 CCC 63.158 1/1 CCC 21.053 1/2 CCC 15.789 2/2 DDD 26.786 2/2 DDD 46.429 2/1 DDD 26.786 1/1 EEE 32.759 2/2 EEE 43.103 2/1

Hi folks Anybody running CentOS on Eee PC 1000H? In CentOS wiki[0] there is a draft about Eee PC 900 from Fabian Arrotin. Does the same howto work for Eee PC 1000H? So far he writes: > Waiting for an official way to host the necessary rpms ... Where to download WiFi and NIC driver? regards Sven [0] http://wiki.centos.org/HowTos/Laptops/Asus/Eeepc

Hi everyone. The function image() seems not to be correctly plotting some matrices that I give it. (I佲檓 using R 2.1.1) The following code creates a matrix (denoted x) with 1500 rows and 3 columns, with all entries 0 or 1, and then plots the image of this matrix. cc=runif(n=1500,min=0.1,max=1.2) ccc=ceiling(cc-1) dd=runif(n=1500,min=0.1,max=1.2) ddd=ceiling(dd-1)

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

Dear R Gurus, I am having a little difficulty with nlm. I've searched the archives and found nothing that tells me why this is occuring -- though there are some slightly similar issues. A simple example: lev2<-function(aaa,bbb,ccc,ddd,eee){ res<-aaa+bbb+ccc+ddd+eee res } nlm(lev2,p=c(32,4,5),ddd=45,eee=23) Error in f(x, ...) : argument "bbb" is missing, with no default