similar to: bootstrapped cox regression in rms package (non html!)

Hi, I am trying to convert a colleague from using SPSS to R, but am having trouble generating a result that is similar enough to a bootstrapped cox regression analysis that was run in SPSS. I tried unsuccessfully with bootcens, but have had some success with the bootcov function in the rms package, which at least generates confidence intervals similar to what is observed in SPSS. However, the

Hello, I am trying to plot two sets of boxplots together. These are estimates of two experiments and?seven?factors. The results of the two experiments I want to plot as boxplots stacked to each other. Therefore I plot first the results of the first experiment; and next with the add option the second set of boxplots. The boxplots are plotted at 'at = 1:7 - 0.15 for the first experiment and

Hello, Warning, I'm guessing only those who have used the Hmisc package's summary.formula function with LaTeX will be able to offer much help here. I am using the Hmisc package's summary.formula function to produce tables for a LaTeX report. The "latex" function in the same package supports longtables in LaTeX. Ideally, I would like for page breaks in the LaTeX output

Hey there, I know a lot of you have already started a side-project or company using Ruby and a lot of you want to. Wanted to share about FounderDating (no, it''s not romantic) since I know finding the right person can be a big challenge. FD is an invite-only, online network for entrepreneurs to start side projects & to connect with cofounders. I see some of you are already on but

Hi, I wish to compute multivariate test statistics for a within-subjects repeated measures design with anova.mlm. This works great if I only have two factors, but I don't know how to compute interactions with more than two factors. I suspect, I have to create a new "grouping" factor and then test with this factor to get these interactions (as it is hinted in R News 2007/2), but

Hi all, I had some trouble in?regrouping factor levels for a variable. After some experiments, I have figured out how I can recode to modify the factor levels. I would now like some help to understand why some methods work and others don't. Here's my code : rm(list=ls()) ###some trials in recoding factor levels char<-letters[1:10] fac<-factor(char) levels(fac) print(fac) ##first

Dear R Users, I am trying to write a wrapper around summarize and xYplot from Hmisc and am having trouble understanding how to pass arguments from the function I am writing to the nested functions. There must be a way, but I have not been able to figure it out. An example is below. Any advice would be greatly appreciated. Thanks, Dan # some example data df=expand.grid(rep=1:4,

Dear List, I am using function recode, from package car, within sapply, as follows: L3 <- LETTERS[1:3] (d <- data.frame(cbind(x = 1, y = 1:10), fac1 = sample(L3, 10, replace=TRUE), fac2 = sample(L3, 10, replace=TRUE), fac3 = sample(L3, 10, replace=TRUE))) str(d) d[, c("fac1", "fac2")] <- sapply(d[, c("fac1", "fac2")], recode, "c('A',

Stefan Roepcke <stefan.roepcke at metagen.de> writes: > Hey, > > hopefully there is an easy way to solve my problem. > All that i think off is lengthy and clumsy. > > Given a data.frame d with columns VALUE, FAC1, FAC2, FAC3. > Let FAC1 be something like experiment number, > so that there are exactly the same number of rows for each level of FAC1 > in the

Hey, hopefully there is an easy way to solve my problem. All that i think off is lengthy and clumsy. Given a data.frame d with columns VALUE, FAC1, FAC2, FAC3. Let FAC1 be something like experiment number, so that there are exactly the same number of rows for each level of FAC1 in the data.frame. Now i would like to scale all values according to the center of its experiment. So i can apply s

Dear all, I'm using the mgcv library by Simon Wood to fit gam models with interactions and I have been reading (and running) the "factor 'by' variable example" given on the gam.models help page (see below, output from the two first models b, and b1). The example explains that both b and b1 fits are similar: "note that the preceding fit (here b) is the same as

Dear all, I'm trying to understand how the strucchange package is working and I have been looking at the examples given for the Fstats() function. The first example (Nile), shows one peak in the F-stats and one breakpoint is estimated, that can be plotted using the following code ## Nile data with one breakpoint: the annual flows drop in 1898 ## because the first Ashwan dam was built

R version 2.10.1 (2009-12-14) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with

On the Digium's site it says: The Wildcard TDM400P is a half-length PCI 2.2-compliant card while for other cards it says: The TE411P is for use only with a 3.3 volt PCI slot. Does the TDM400 not only fits, but also functions in a 3.3V only slot? >From what i detected so far, is that some MOBO manufactures have pci-slots that provide 3.3 Volt AND 5.0 Volt, thus can handle all kind of

Can anyone inform me if voip can be used on a IPv6 network? Does any hard phones/soft phones/Asterisk support it? Google told me that there was/is a bounty on it, but that expired august last year. Furthermore, there used to be a patch (Bernhard Schmidt), but that one is about a year old. I presume it can't be used on recent versions of "*" Hans -- pgp-id: 926EBB12

Hi all, Just curious, Quite a while a go, i was checking for supported SW-platform. AFAIR, it was RHES and SLES Now it's only RHES-4 and FC-3 or FC-4. Not a single syllable about CentOS or SLES-9 or SLES-10 It probably just runs fine, but any chance of getting support for their *-enterprise version? (just in case of, if one needs it) Hans -- pgp-id: 926EBB12 pgp-fingerprint: BE97 1CBF

Just curious, I noticed that with SetLanguage() you can change it into a lot of other languages. Yes one can record them easily enough with "record", but don't like to re-invent the wheel.. Browsed through a lot of google-pages but failed to find any other languages (except for FR and ES on the digium site) Any pointers for german, dutch, greek, italian, .... prompts? Hans --

First encounter with * Just downloaded & installed aah-2.7 Started up AMP, but i can not find any reference towards isdn. I presume there has to be some configuration done for my Eicon-Diva-pro. Does aah actually support isdn-bri? On the mail-archive i found some references, but these are rather old ( they speak about the coming release of aah-2.1) aah-handbook (version 1.6) doesn't

Hello, Using a fixed effects linear model (with lm), I can get exact p-values out of the AVOVA table, even if they are very small, eg. 1.0e-200. Using lme (linear mixed effects) from the nlme library, it appears that there is rounding of the p-values to zero, if the p-value is less than about 1.0e-16. Is there a way we can obtain the exact p-values from lme without rounding? used commands:

Hi all, I have just started playing with speex, and come up with the following code, which just encode a frame of 160 shorts, and the decode it. For some reason the decoded sample is completely different than the original one. is my code wrong? If so what? Or is it a reasonable which depends of values that weren't correctly set? Thanks, Andre #include <stdio.h> #include