similar to: brookmeyer & crowley CI for median survival time

Hello all, I have a very simple set of data and I would like to analyze them with probit analysis. The data are: X Event Trial 100 8 8 75 8 8 50 6 8 25 4 8 10 2 8 0 0 8 I want to estimate the value of X that will give a 95% hit rate (Event/Trial) and the corresponding 95% CI. Anyone can offer some help? Thanks!! -

Dear R helpers, I have a 3D scatter plot that I have generated from scatterplot3d (which looks great- thanks!) and I can see that the points in my graph fall in a plane. Following the example 5 from 3D scatter plot (below) I have fitted a regression plane. Now what I would like to do is a rotation so that my new co-ordinate system is about the fitted plane (by finding the normal to the plane

Got rsync daemon working, but I am trying to set some variables in the rsyncd.conf file: here is my rsyncd.conf file: use chroot = false strict modes = false hosts allow = * log file = cygdrive/f/rsyncd.log uid = 0 gid = 0 # Module definitions [BackupArea] path = /cygdrive/f/backup/test/%computerName% <---------- This is how I want it to work... read only = false transfer logging = yes

Hi, Storer and Crowley (JASA 1985) presented an approach for approximating the changes in maximum partial-likelihood parameter estimates for the Cox model when a single observation is deleted. Is there an R implementation of this approach? Any help is greatly appreciated. Thanks. Best, Ravi. ---------------------------------------------------------------------------- ------- Ravi

Dear R users, I'm interested in measuring the stability of a heirarchical clustering, of the overall clustering and finding sub clusters (from cutting the heirarchical clustering at different levels) which demonstrate stability. I saw some postings on the R help from a while back about bootstrapping for clustering (using sample and generating a consesus tree with a web based tool CONSENSE)

Dear list, I want to create a POSIX time vector as follows: day <- as.character("110809") time.t <- 1:3600 t.min <- time.t %/% 60 t.sec <- time.t-t.min*60 DATE <- as.POSIXct(strptime(paste(day,t.min,t.sec),"%y%m%d %M%S")) Tail(DATE) The problem is that the last element (3600) returns a NA and I don't understand why. 600, 1200, 2400

Hi, Thanks in advance for the help. I have a few vectors and I need to get all the combinations of the vector values. Normally I would use a few for-loops to do that. e.g x1_c(1,3,5) x2_c(2,5,6) for (i in x1) { for (j in x2) { ... }} My problem is that I don't know how many vectors there are ahead of time, so I don't know how many for-loops I need. Any

I am using Rsync in a Windows environment and its working perfectly. I am running into an issue where on one of my Windows servers there is a folder called "share" rsync sees this and copies it exactly the way its supposed to. Here is my rsync command: rsync -auvrlt --progress --stats --delete "/cygdrive/v/SHARE/" "jaxback02::u/%computername%/SHARE/ The issue is

Hello everyone, Does anyone know if I can do Breslow Day Test for the homogeneity of odds ratio in R? Thanks! - Jacqueline

Hi all, I am trying to run the k-means cluster analysis using the function kmeans in the package cluster. The data are: x = c(-0.26, -0.23, -0.05, -0.20, 0.30, -0.84, -0.10, -0.12, 0.10, -0.31, -0.19, 0.18, -0.26, -0.23, -0.37, -0.23) I've got two different solutions when I ran this function over a few times: kmeans(x, centers=2) The first solution gives the following: $cluster [1]

Dear all, I have a question on the dose-response estimation with clustered/ correlated binary data. I would like to estimate the hit rate for a certain test at various concentration levels. The test is used on 5 subjects, and each subject is tested 20 times. If we assume that the 100 samples are independent, the hit rate estimate is unbiased, but the variance is under-estimated. The other

Dear R list, I have a package downloaded and looked up a function in there. Now I find that it uses C code (.C call) to do part of its job. However, when I wanted to look that part up by using getAnywhere() I was told that no object of that name could be found. Then I tried typing C.(" funct.name") which only caused R to crash. Is there any way to look up the function? I'm not

When running the cupsaddsmb command, I get the following: Since I am running as root, what do I not have access to and why? Any help would be appreciated. Glen Crowley <<...OLE_Obj...>>

Dear list, I have problem that I cannot solve and would like to ask your opinion. I tried to ask a few days ago already but got no answer and all my attempts to solve it by myself since then failed. Sorry for repeated posting! Here the problem broken down a bit. My problem basically is, that I want to use the elements of a character vector as names for objects and by recalling only the

> -----Original Message----- > From: Jacqueline Lee [mailto:jlee@isdomaininc.com] > Sent: Friday, September 23, 2005 11:46 AM > To: asterisk-users@lists.digium.com > Subject: channel offhook state > > > We are using a digium card (TDM400) with asterisk for our access to the > PSTN. Initially when the server starts, all the zap channels on the card > are in the

I need help with replacing NaN with zero (the value '0') in my dataset. The reason is that I can't get it to graph because of the NaN in the dataset. I have tried: data[is.nan(data)] <- 0 that others have suggested in the help archives but this does nothing so I am not sure what I am doing wrong. Thanks in advance for your help. [[alternative HTML version deleted]]

Hi all I have a multiple datasets that look like this MM Freq 1 30 2 35 3 54 4 33 5 27 6 13 7 25 8 29 9 40 10 32 11 36 12 23 I am plotting this using barchart (there is probably something better but it gives me the results I want) and I would like the x-axis to have the names of the months instead of the month numbers. So I have

Changes with nginx 0.8.33 01 Feb 2010 *) Security: now nginx/Windows ignores trailing spaces in URI. Thanks to Dan Crowley, Core Security Technologies. *) Security: now nginx/Windows ignores short files names. Thanks to Dan Crowley, Core Security Technologies. *) Change: now keepalive connections after POST requests are not

I have a data set that I have imported (not sure if that makes a difference) and I would like to calculate the sum of only specific columns. I have tried >colSums(dataset, by=list(dataset$col5), dims=1) and I get an error of unused arguments I have also tried >aggregate(dataset, by=list(dataset$col5), sum) and I get the error that sum is not meaningful for factors. I want to only calculate

Hi all I have a bit of a problem. I want to make a barplot of some data. My data is of a score that is separated by year and by a limit (above 3 and below 3 to calculate the score). Year Limit HSS 1999 ALT 0.675 1999 VFR 0.521 2000 ALT 0.264 2000 VFR 0.295 I would like to have a barplot with year on the x axis and HSS on the y axis and the two limits as two different colors to show the