search for: estimt

...haps someone will have a solution to my more general problem, but here is the specific one: I used the round() function to round some estimates to 3 decimal places. I then sent put the rounded estimates in a matrix and used latex() to make a LaTeX table from them. However, in my table, there are estimtes which only have 2 decimal places. I assume that the third decimal place in these numbers was equal to zero, so the round() function simply cut off the superfluous zero. However, I want all of my estimates to have the same number of digits. Is there a way to force R to keep all three digits, eve...

Here's one possibility, assuming muhat and sigmahat are estimtes of mu and sigma from N iid draws of N(mu, sigma^2): tStat <- abs(x - muhat) / sigmahat pValue <- pt(tStat, df=N, lower=TRUE) I'm not quite sure what df tStat should have (exercise for math stat), but given fairly large N, that should make little difference. Andy > From: Ross Clem...

I need to find out the difference between the way R calculates weighted regression and standard regression. I want to plot a 95% confidence interval around an estimte i got from least squares regression. I cant find he documentation for this ive looked in ?stats ?lm ?predict.lm ?weights ?residuals.lm Can anyone shed light? thanks Chris. -- View this message in context: http://n4.nabble.com/weighted-least-squares-vs-linear-regression-tp1387957p1387957.h...

Dear all, Suppose I have a linear mixed-effects model (from the package nlme) with nested random effects (see below); how would I present the results from the random effects part in a publication? Specifically, I?d like to know: (1) What is the total variance of the random effects at each level? (2) How can I test the significance of the variance components? (3) Is there something like an

Hi, I am computing the Nelson-Aalen (NA) estimate of baseline cumulative hazard in two different ways using the "survival" package. I am expecting that they should be identical. However, they are not. Their difference is a monotonically increasing with time. This difference is probably not large to make any impact in the application, but is annoyingly non-trivial for me to just

...for classification ctest Classical Tests eda Exploratory Data Analysis gee Generalized Estimating Equation models integrate numerical integration leaps Regression subset selection locfit Local Regression, Likelihood and Density Estimtion. mva Classical Multivariate Analysis polymars (polychotomous) regression based on Multivariate Adaptive Regression Splines polynom A collection of functions to implement a class for univariate polynomial manipulations spatial fun...

...for classification ctest Classical Tests eda Exploratory Data Analysis gee Generalized Estimating Equation models integrate numerical integration leaps Regression subset selection locfit Local Regression, Likelihood and Density Estimtion. mva Classical Multivariate Analysis polymars (polychotomous) regression based on Multivariate Adaptive Regression Splines polynom A collection of functions to implement a class for univariate polynomial manipulations spatial fun...

Dear list members, I try to simulate an incomplete block design in which every participants receives 3 out of 4 possible treatment. The outcome in binary. Assigning a binary outcome to the BIB or PBIB dataset of the package SASmixed gives the appropriate output. With the code below, fixed treatment estimates are not given for each of the 4 possible treatments, instead a kind of summary

...difference. There are ties in failure times in the data set. Consequently, it matters which method is used to handle ties in "coxph". The default is "efron". If I use method="breslow", there is no difference between the 2 different ways of computing the Nelson-Aalen estimte. require(foreign) gb <- read.dta("GB.dta") # Green & Byar data; N = 483 # Method 1 (note the use of Breslow estimator) fit1 <- coxph( Surv(time, status=="Cancer" | status=="CVD" | status=="Other") ~ 1, method="breslow", data=gb) h1...