similar to: inclusion criteria help

I have managed to create a data frame like this: > tsus_same_mean[1:10,] PID Grp Dist PercAln PercId 1 12638 Acidobacteria 0.000000000 1.0000000 1.0000000 2 87 Actinobacteria 0.000000000 0.9700000 0.9700000 3 92 Actinobacteria 0.008902000 1.0000000 0.9910000 4 94 Actinobacteria 0.000000000 1.0000000 1.0000000 5 189 Actinobacteria 0.005876733

Hi! I have a set of 49 pairwise comparisons that I have done. From this I would like to plot either histograms or the density plots of the values I get. Now, I can plot one histogram per comparison, but I have problems getting the output I want. When plotting like I normally would do: histogram(~percid | orgA_orgB, data = alldata) I get the histograms next to eachother in a boxlike shape.

Dear R-listers Is there anyone who knows why I get different eigenvectors when I run MatLab and R? I run both programs in Windows Me. Can I make R to produce the same vectors as MatLab? #R Matrix PA9900<-c(11/24 ,10/53 ,0/1 ,0/1 ,29/43 ,1/24 ,27/53 ,0/1 ,0/1 ,13/43 ,14/24 ,178/53 ,146/244 ,17/23 ,15/43 ,2/24 ,4/53 ,0/1 ,2/23 ,2/43 ,4/24 ,58/53 ,26/244 ,0/1 ,5/43) #R-syntax

I have been attempting to use cph models to get better calibration of my models for which I had originally used logistic regression. I tried running with 40 repetitions and got an error. I then tried 500 repetitions (thinking that the NaNs in the output below might be caused by that choice) and then let my computer crunch for several hours and got only the same error message and

Dear list members, I am using Mclust in order to deconvolute a distribution that I believe is a sum of two gaussians. First I can make a model: > my.data.model = Mclust(my.data, modelNames=c("E"), warn=T, G=1:3) But then, when I try to plot the result, I get the following error: > mclust1Dplot(my.data.model, parameters = my.data.model$parameters, what = "density")

Dear R users, I noticed a problem in the anova command when applied on a single coxph object if there are missing observations in the data: This example code was run on R-2.6.1: > library(survival) > data(colon) > colondeath = colon[colon$etype==2, ] > m = coxph(Surv(time, status) ~ rx + sex + age + perfor, data=colondeath) > m Call: coxph(formula = Surv(time, status) ~ rx +

Hi, When I run the survfit function, I want to get the restricted mean value and the standard error also. I found out using the "print" function to do so, as shown below, print(km.fit,print.rmean=TRUE) Call: survfit(formula = Surv(diff, status) ~ 1, type = "kaplan-meier") records n.max n.start events *rmean *se(rmean) median 200.000

Hello R People: When using the arima function with the seasonal option, are the seasonal options only good for monthly and quarterly data, please? Also, I believe that weekly and daily data are not appropriate for seasonal parm estimation via arima. Is that correct, please? Thanks, Sincerely, Laura Holt mailto: lauraholt_983 at hotmail.com download!

I have a one-to-many relation between my documents and statuses tables. I want to display only the documents whose states are not "Published". I haven''t figured out how to do this using the standard paginate method (maybe there''s a way using its :joins option). I''m also considering writing a custom paginator or using paginate_by_sql[1]. Does anyone have

Is there an R function that I can use to calculate the p-value from the Z statistics computed for the relationship between chance and observed proportions in predictions. More sprcifically I am refering to proportional chance criteria (Cpro). Details are in Huberty's book on Applied discriminant analysis, unfortunately our library has misplaced the book. I got some details from the

Dear R-Helpers Has anyone written a script for the dose proportionality assessement using the power model and the confidence interval criteria as described in - Gough at al. Assessment of dose proportionality: report from the statisticians in the pharmaceutical industry in Drug Information Journal,1995, 29:1039-1048 and/or - Smith et al. Confidence interval criteria for assessment of dose

Hi, I have a data frame that I want to split into two based on the values of a variable in it. The variable Y has numeric values ranging between 0 through 70. I want to plot the frequencies of another variable X in two different cases: - When Y = 0 and - When Y > 0 How does one go about doing this? In general, I want to do several analyses with this data frame that are a variation of the

Hi all! I'm running some simulations and I need to estimate some paramaters with optim( ), in some cases optim stops with the next message: "L-BFGS-B needs finite values of 'fn' " I would like to know how to include and "if" condition when this happen, could it be something like: myfun <- optim(....) # run my function

Hello, I'd like to clarify a few points on how the release process works. Each release must satisfy the following criteria for each supported target: * LLVM-GCC & LLVM must build in both release and debug mode. They must also build srcDir != objDir. LLVM-GCC must build with support for c, c++, fortran, and obj-c/obj-c++ (Mac only). * LLVM-GCC must bootstrap. * "make check"

http://bugzilla.mindrot.org/show_bug.cgi?id=624 djm at mindrot.org changed: What |Removed |Added ---------------------------------------------------------------------------- Status|NEW |RESOLVED Resolution| |WONTFIX ------- Additional Comments From djm at mindrot.org 2004-12-20 11:22

Hi guys, I've got a quick feature suggestion to solve a problem that I don't think is readily solvable having to do with an "after hours" message, playing only during (or rather outside) specified times. I know all about the helpful feature that already exists which allows you to tack on a cron-ish specification of time, date, month, etc. to the end of an "include

Bojan wrote: > Is there a list of criterias for a 1.0 release of Wine? I couldn't find > one. I was looking in the developers mailing list. But there's only a > collection of bugs that must be fixed. ... http://winehq.org/?issue=320 has a summary of Alexandre's talk at Wineconf recently. It said in part: "... a 1.0 release seems to be near... The current plan is to

Greetings and thank you CentOS team for very quick production and maint across all version number releases !!! In regards to 4.4, is the CentOS Server_CD for 4.4 already in production and testing? May we know the ETA for it on the mirrors please? If it will be quite awhile, may we know more about how it is produced and any release critera you have for it please? Thank you and kind regards! -

gcc is a great compiler, but as of gcc3.2.2, it can't compile Samba on Sparc with -O2. This has been marked as a high priority bug in the gcc bug tracking system, and should be fixed before the next release. To make such regressions less likely in the future, it would be helpful if somebody were to look at http://gcc.gnu.org/testing/ and submit a build-and-test guide for Samba to the