similar to: simple coding question

Hello. I am looking for a function that will allow me to paste rows together without regard for the numbers of columns in the datasets to be joined. The only columns where it matters if they are aligned correctly are at the beginning - the rest of the columns represent differing numbers of ICD9 (disease) codes reported by each person(record) at a health visit. They are in no particular order.

Hello everyone, I have a .csv file with the following format: uniqueID SubjectID Distance_miles Tag 1 1001 5.5 3 2 1001 7 1 3 1001 6.5 1 4 1001 5 1 5 1002

Kirsten, this may not be the most elegant approach, but I think I would create dummy variables on the fly for each disease code you have, then use the "with" function to get means. Give this a try: # setup some fake data as an example icd9 <- c(rep("dis1",4), rep("dis2",5), rep("dis3",3)) n <- length(icd9) exposure <- rnorm(n) working <-

I would like to code records in a dataset with a 1 if any of the columns 9-67 contain a particular code, and zero if they don't. I've been working with "subset" and it seems that something like subset(data, data[9:67]--"12345") would work, but I have been unsuccessful so far. It seems like a simple problem - any help is appreciated!

Hello: I am working on getting some statistics related to clinical trials and stuff. I have to work with ICD9 codes. Is anyone aware of any R method that deals with ICD9 codes verification and manipulation. Thanks Vishwanath

I am interested in locating a script to implement a sampling scheme that would basically make it more likely that a particular observation is chosen based on a weight associated with the observation. I am trying to select a sample of ~30 census blocks from each ZIP code area based on the proportion of women in a ZCTA living in a particular block. I want to make it more likely that a block will

Dear R-people: François Michonneau's method to obtain the special character plus/minus works on Windows 2000 professional. Many Thanks to François for his work! Phil Smith Centers for Disease Control and Prevention Atlanta, GA _____ From: François MICHONNEAU [mailto:francois.michonneau@gmail.com] Sent: Monday, March 20, 2006 1:58 PM To: r-help@stat.math.ethz.ch; Smith, Phil

Kia ora I'm a using lme (from nlme package) with data similar to the Orthodont dataset and am getting positive log-likelihoods (>100). This seems usual and I wondered if someone could offer a possible explanation. I can supply a sample dataset if requested, but I feel almost certain that this question has been asked and answered recently. However, I can find no trace of it in the mail

I have a dissimilarity dataset with the form: 1 1 dissimilarity value 1 2 ... 1 3 1 4 2 2 2 3 2 4 ... I would like to do nonmetric multidimensional scaling with this data, but I am having trouble using this format. I would like to either find a function that accepts this format or find a way to easily convert this format to a matrix for use with existing functions. Thanks!

Hi all, I'm using the medical risk package to determine the risk stratification Index based on the ICD9 codes. Although, I have been successful in using it, I'm unable to interpret the output. here is the sample code # Calculate RSI for each patient ("id") in dataframe cases <- data.frame(id=c(1,1,1,2,2,2,2,2), +

Excuse the simple question... I'm not sure what I'm doing wrong with predict, but let me use this example: Suppose I do: dat<-matrix(c(0,0,10,20),2,byrow=T) lm1<-lm(dat[,2]~dat[,1]) Suppose I want to generate the linearly-interpolated y-values between the point (0,0) and (0,20) at every unit interval. I thought I just do: predict(lm1, data.frame(seq(0,10,1))) to get

Kia ora list members: I'm having a little difficulty getting the correct standard errors from a glm.object (R 1.9.0 under Windows XP 5.1). predict() will gives standard errors of the predicted values, but I am wanting the standard errors of the mean. To clarify: Assume I have a 4x3x2 factorial with 2 complete replications (i.e. 48 observations, I've appended a dummy set of data at the

Hi again, I changed from 3.4.5 to 3.3.10 after an AD upgrade to 2008R2 and prinitng is working mor ore less but it is incredibly sloooow. My server has no remarkable load so it is not due to hardware bottlenecks. Are there any ways to 'tune' 3.3.10 or any tricks to make it faster ? cheers ~christoph -- /* Christoph Beyer | Office: Building 2b / 23 *\ * DESY

Dear all, I have a data.frame which has names as following. [1] "XG1" "YG1" "XEST" "YEST" [2] "XNOEMP1" "XNOEMP2" "YNOEMP1" "YNOEMP2" [3] "XBUS10" "XBUS10A" "XBUS10B" "XBUS10C" [4] "YBUS10" "YBUS10A" "YBUS10B" "YBUS10C" [5]

Hi Jerry et al, will the 2.2.6 release show the printers in the 'add-new-printer-wizzard' on the client like the WINDOWS server do ? I'm not sure if you know what I mean, so here is a little illustration: ----------------------------- | NTSERVER | | \\ntserver\printer1 | | \\ntserver\printer2 | | \\ntserver\printer3 | | SAMBASERVER |

Hi, my windows folks migrated to AD 2008 R2, resulting in the following error message when trying to join the domain: [HOST] /etc $ /opt/csw/bin/net ads join -U <USER> Enter <USER>'s password: [2009/07/01 11:51:28, 0] libads/sasl.c:ads_sasl_spnego_bind(819) kinit succeeded but ads_sasl_spnego_krb5_bind failed: Strong(er) authentication required Failed to join domain: failed

Hi there, I use 3.5.2 with netbios aliases in security mode ADS. When I connect to the primary netbios name from a windows host I get the samba machine right away. Connecting to the netbios alias (which is on a differen IP address) I have to type name and pwd to get the connection. DNS entries are OK. Is there any way to get rid of the pwd question, I had the same setup with 3.0.20 for a

There appears to be an issue with the decreasing=TRUE option on the order() function that indicates either a bug or perhaps a design flaw (potentially flawed because I would suggest the majority of users would expect different behaviour). # demonstration of problem: x <- c(2,1,3,4,5) order(x) order(x, decreasing=TRUE) order(x) correctly reports the order as: 2 1 3 4 5 I expected the result

Hello. Suppose you define a new S4-class, say > setClass("track", representation(x="numeric", y="numeric")) Don't worry if you have a deja vu, it's from the help page. Your new class is said to have a fixed structure: two slots, x, and y, and that should apply to all objects you construct as members of that class. > tr <- new( "track" )

Hi, I found some hints but no solution that worked for me actually, I try to compile 3.3.5 on Solaris 10 and it does not find libtalloc: Linking shared library bin/libtalloc.so.1 /usr/ccs/bin/ld: cannot open linker script file /scratch/samba-3.3.5/source/exports/libtalloc.so.1: No such file or directory gmake: *** [bin/libtalloc.so.1] Error 1 Obviously because it's not there :( :