search for: obst

...2007' and 'Apr 1, 2007') group by id,tshour order by id,tshour; I've pulled data from PostgreSQL into R, and have a dataframe containing a timestamp column, v, and o (both numeric). I created an grouping index for every 5 minutes along these lines: obsfivemin <- trunc(obsts,units="hours") +( floor( (obsts$min / 5 ) ) * 5 * 60 ) (where obsts is the sql timestamp converted into a DateTime object) Then I tried aggregate(df,by=obsfivemin,cor), but that seemed to pass just a single column at a time to cor, not the entire data frame. It work...

...Internal Medicine University of Virginia, Charlottesville, Virginia ------------------ ######################################################## # Compute the Hosmer-Lemeshow 'goodness-of-fit' test cd.full_model = glm(formula = Collaterals ~ CHF + Age + CABG + relevel (as.factor (num.obst.vessels),"one") + Current.smoker + DM + HTN + ace.inhibitor + MI, family = binomial(link = "logit")) hosmerlem = function(y, yhat, g=10) { cutyhat = cut(yhat, breaks = quantile(yhat, probs=seq(0, 1, 1/g)), include.lowest=TRUE) obs = xtabs(cbind(1 - y, y) ~ c...