similar to: Aggregate with non-scalar function

Hi, I'm reading van Belle et al "Biostatistics" and trying to run a cox test using a dataset from: http://faculty.washington.edu/~heagerty/Books/Biostatistics/chapter16.html (Primary Biliary Cirrhosis data link at top of the page), I'm using the following code: --------------- start of code library(survival) liver <-

R 2.10 Windows 7 I am trying to plot three graphs on top of each other. I need to have the axises perfectly aligned. For some reason the ticks on the y axes are slightly off so they do not perfectly align. Can someone tell me how I can get the to overlay each other perfectly? I thought the yaxp parameter would solve my problem, but it does not. My data and code follows: >

Hello Everyone, I am relatively new to linux and Samba and I am having some difficulties setting Samba up as a Primary Domain Controller. Ihave read the Pdf document called "Using Samba as a PDC" by developrWorks and I have read most of the Howto documents relating to Samba. When I try to connect to my samba server with Win98se I get the message "no Domain server was able to

Hallo, I have got an existing data frame and want to add a new column. The existing data frame was created like this: > df <- rbind( c("Fred", "Mary", 4), c("Fred", "Mary", 7), + c("Fred", "Mary", 9), c("Barney", "Liz", 3), + c("Barney", "Liz", 5) ) > df

Updated my Athlon-xp 6-stable system last night, got an em watchdog timeout for the first time a few hours later, during a fairly high-traffic period. System is UP but does have device apic in the config. Any chance this is the recent race condition? Workaround? ifconfig em0 down, ifconfig em0 up seemed to cure it, at least for the moment. Thanks, Barney Wolff -- Barney Wolff I never

Hi. I''ve installed ruby, gems, rails into a user directory according to the installation instructions (for linux). I haven''t setup any path information nor setup any symlinks or any compilation options. When I try rails, I find that the system accesses controllers but that the corresponding views are not displayed. I wonder if there is something that I forgot to configure

Dear all, I'm looking for someone that help me to write an R function to simulate survival data under complex situations, namely time-varying hazard ratio, marginal distribution of survival times and covariates. The algorithm is described in the reference below and it should be not very difficult to implement it. However I tried but without success....;-( Below there the code that I used; it

I guess I understand it now (although it is surprising to me). The following is a valid model formula fred ~ wilma + offset(barney) that sets the model offset to barney. Given that this works, it would seem that one could remove the offset = argument from the glm call (and document the offset feature somewhere). Too bad that one can't set weights the same way. The anova bug when offsets

Hallo, I build a list by the following way: Lst = list(name="Fred", wife="Mary", no.children=3, cild.ages=c(4,7,9)) I know how I can extract the information one by one. But now I want to add a new entry which looks like name="Barney", wife="Liz", no.children=2, cild.ages=c(3,5) How can I add this information to Lst without overwriting the first entry?

Last week one of our clients reported trouble with a csv file I generated with write.table. He said that columns with quotes for character variables were rejected by their data importer, which was revised to match the way Microsoft Excel uses quotation marks in character variables. I explained to them that quoted character variables are virtuous and wise, of course, but they say Microsoft Excel

Hi all: I am a user of "JM" package. Here's the problem of "sample size". The warning is: Error in jointModel(fitLME, fitSURV_death, timeVar = "time", method = "piecewise-PH-GH") : sample sizes in the longitudinal and event processes differ. According to the suggestion of "missing data",I use the same data set(data_JM) without any

Hi - this is my first day with Ruby on Rails. I''m starting out by working through the Helloworld (Hello Ruby) example program printed up in "Agile Web Development with Rails." My steps are as in the text: 1. ruby script/generate controller Say 2. edit the proper file and add the hello method 3. fire up the browser using address http://localhost:3000 4. fire up the

I''ve installed the gem "meta_where" (and it reported back: "Succeeded") but when I use the example in a search method of a controller: Person.where(:some_field.matches => ''%something%'') there is an error that it doesn''t know the method ''matches'', so what still needs to be done to properly activate the gem and make its

Still very new to puppet, so go easy :-) Bit of a braindump I''m afraid, but I''d really appreciate any "Yep", "interesting", "No!!!" comments. I''ll start with an example - adding NTP support. This primarily consists of two configs: a) NTP clients - this is quite easy, it''s a templated file requiring the names of the NTP

On p. 299 of the current Rails book it says that "Active Record takes care of saving all the dependent child rows when you save a parent row." So, when I add a row to the table that has "belongs_to" in its model, shouldn''t the table with the "has_many" be updated with data from the foreign key? If not, how do I save data from a form into more than one

Hi R-users,   I wrote a code to evaluate double sum as follows:   ff2 <- function(bb,eta,z,k) { r <- length(z) for (i in 1:r) { sm1 <- sum((z[i]*bb/2)*(psigamma((0:k)+eta+1,deriv=0)/(factorial(0:k)*gamma((0:k)+eta+1))))  sm2 <- sum((besselI(z[i]*bb,eta)*log(z[i]*bb/2) - sm1)/besselI(z[i]*bb,eta))  sm2 } ff2(bb,eta,z,10)     but it gave me the following message:   >

Hi, I seem to be having trouble getting more than 10 hits from Index#search_each since upgrading to 0.10.2 (ie, this was working in 0.9.4). Maybe a bug, as the #search_each doesn''t seem to use the options parameter any more ? Thanks, Neville =========================================== require ''rubygems'' require ''ferret'' p Ferret::VERSION idx =

Hi Step 1: I create a data.frame called iolm. Step 2: I create a conditional subset i_wtr. Step 3: In this subset I add 0.3 to all values in the IOLM_AST column. Step 4: Now I am looking for the best way to Œmerge¹ the altered subset back into the original iolm data.frame ## STEP 1 >iolm ID IOLM_AST IOLM_AXIS 1 1 1.15 165.33 2 2 1.20 79.00 3 3 0.40

I will be upgrading my Samba server from 2.2.8a to 3.0.10. I currently have security set to 'share' and plan on migrating to 'ads' for improved authentication. I have one snag, though... I have remote users who reside in and are managed by a Windows domain that is not in my control. There is no trust relationship at all. If I use 'ads' security, can I add a 'valid

I am having a problem with my samba config on my web server related to file permissions. The problem seems quite subtle! The samba config file contains (Samba 2.0.7): ----cut here [global] workgroup = UKC server string = University Web Server security = domain encrypt passwords = Yes password server = KNOLE CORFE log file = /var/adm/log/samba.log