Re: [R] Sweave files generating miktex errors

2013-01-30 Thread Troy S
the .sty files somewhere to help the situation. Thanks guys. I appreciate the help! Troy On Wed, Jan 30, 2013 at 5:52 AM, Michael Friendly frien...@yorku.ca wrote: On 1/29/2013 6:04 PM, Troy S wrote: Ok, yes I realize it. So let me try to fix it: I removed the sty files, and set TEXINPUTS

Re: [R] Sweave files generating miktex errors

2013-01-30 Thread Troy S
Duncan, Alas, I had left out the --pdf flag when I tried the solution. I am glad I have an answer now, and a better understanding of the problem. Troy On Wed, Jan 30, 2013 at 9:50 AM, Duncan Murdoch murdoch.dun...@gmail.comwrote: On 13-01-30 11:00 AM, Troy S wrote: Michael, I added C:\R

[R] Sweave files generating miktex errors

2013-01-29 Thread Troy S
Dear useRs-- I have been using Sweave with miktex for years, but on a new install on Windows XP, miktex seems to be hung up on single quotes. See example below. Digging through stackexchange, I found using \usepackage[noae]{Sweave} in the tex file solved the problem. My questions are: --Why

Re: [R] Sweave files generating miktex errors

2013-01-29 Thread Troy S
Duncan, Copying all 3 sty files to my working directory solved the problem. I'm still curious why this was necessary in my new install. Duncan, many thanks! Troy On Tue, Jan 29, 2013 at 10:58 AM, Duncan Murdoch murdoch.dun...@gmail.comwrote: On 29/01/2013 1:37 PM, Troy S wrote: Dear useRs

Re: [R] Sweave files generating miktex errors

2013-01-29 Thread Troy S
murdoch.dun...@gmail.comwrote: On 13-01-29 4:55 PM, Troy S wrote: Duncan, Copying all 3 sty files to my working directory solved the problem. I'm still curious why this was necessary in my new install. Duncan, many thanks! You do realize that you have now done a Very Bad Thing, I hope

[R] Chi-squared test when observed near expected

2012-12-03 Thread Troy S
Dear UseRs, I'm running a chi-squared test where the expected matrix is the same as the observed, after rounding. R reports a X-squared of zero with a p value of one. I can justify this because any other result will deviate at least as much from the expected because what we observe is the

[R] Using mixed models to analyze Longitudinal intervention

2011-08-18 Thread Troy S
Dear R List, I am trying to use mixed models to analyze an intervention and want to make sure I am doing it correctly. The intervention is for lowing cholesterol and there are two groups: one with an intervention and one without. The subjects were evaluated a differing amount of time, so there

[R] Implementing step-wise linear regression

2011-01-23 Thread Troy S
Dear R fans, I am trying to do step-wise linear regression using the F-test to decide which variables to admit. Ewout Steyerberg suggests using the F-test for this purpose. I first build a model using no variables using lm(y ~ 1) and then using one variable that is a strong predictor using lm(y

Re: [R] Confidence Intervals for logistic regression

2010-08-07 Thread Troy S
originale- Da: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]per conto di Troy S Inviato: Friday, August 06, 2010 6:31 PM A: Michael Bedward Cc: r-help@r-project.org Oggetto: Re: [R] Confidence Intervals for logistic regression Michael, Thanks for the reply. I believe

[R] Confidence Intervals for logistic regression

2010-08-06 Thread Troy S
Dear UseRs, I have fitted a logistic regression using glm and want a 95% confidence interval on a response probability. Can I use predict(model, newdata, se.fit=T) Will fit +/- 1.96se give me a 95% of the logit? And then exp(fit +/- 1.96se) / (exp(fit +/- 1.96se) +1) to get the probabilities?

Re: [R] Confidence Intervals for logistic regression

2010-08-06 Thread Troy S
Troy S troysocks-tw...@yahoo.com Dear UseRs, I have fitted a logistic regression using glm and want a 95% confidence interval on a response probability. Can I use predict(model, newdata, se.fit=T) Will fit +/- 1.96se give me a 95% of the logit? And then exp(fit +/- 1.96se

Re: [R] Splom custom superpanels

2008-08-11 Thread Troy S
R Gurus-- I am trying to adopt this code. I do not care about colors, but would like the data in the upper plots jittered. Here is what I have. So far, the lower squares are blank. Troy print(splom(foo, upper.panel = function(x, y, ...) {panel.splom(jitter(x), jitter(y), ...)},

[R] Lattice : Single reference line/label at 0.05

2008-06-30 Thread Troy S
r-friends-- I am trying to make a lattice plot with a single label at 0.05. A reference line at y=0.05 as well... This is what I have so far. The label reads 0.0 and there is a grid instead of a single reference line. Can someone help please? Troy x-18:35 (pv-1/(x-17)) xydf-data.frame(x,

[R] Help with stat.table in Epi package,

2008-06-13 Thread Troy S
R Fans-- I am having problems with the following code. It worked under R 2.6.0 but not in 2.7.0. library(Epi) df - read.table( c:/Documents and Settings/Troy S/My Documents/debug_chisq_080613b.txt) summary(df) cvd agecat Min. :0. (0,40] :1 1st Qu.:0.