[R] likelihoods in SAS GENMOD vs R glm

2007-03-19 Thread Brian S Cade
-0.0784 42.81 .0001 n 1 0.0001 0.0163 -0.0318 0.0320 0.00 0.9941 Scale 0 1. 0. 1. 1. NOTE: The scale parameter was held fixed. Brian Brian S. Cade U. S. Geological

Re: [R] likelihoods in SAS GENMOD vs R glm

2007-03-19 Thread Brian S Cade
To answer my own question since no one else has responded: a little investigating shows that proc GENMOD eliminates the sum of the log (x !) in its calculation of the log likelihood for Poisson regression. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre

Re: [R] enhanced question / standardized coefficients

2007-02-07 Thread Brian S Cade
in trying to equate effects of variables measured on disparate scales. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL PROTECTED] tel: 970 226-9326 John Fox [EMAIL PROTECTED] Sent by: [EMAIL

Re: [R] heteroscedasticity problem

2007-02-07 Thread Brian S Cade
to use to estimate these relationships. If you have trouble getting started or finding recent pubs, contact me off list and I will be glad to point you in productive (pun intended) directions. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C

Re: [R] Looking for transformation to overcome heterogeneity of variances

2006-08-03 Thread Brian S Cade
. A better analysis but still not as enlightening as possible. My 2 pence. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL PROTECTED] tel: 970 226-9326 Paul Smith [EMAIL PROTECTED] Sent by: [EMAIL

Re: [R] MRPP in R

2006-05-05 Thread Brian S Cade
data frames, so you can cross back and forth between statistical packages as needed (though not as conveniently as if MRPP were an R package). Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL

Re: [R] non parametric MANOVA

2006-03-29 Thread Brian S Cade
friendly implementation in the USGS Blossom software (www.fort.usgs.gov/products/software/blossom/blossom.asp) or obtain the raw Fortran programs from Paul Mielke's web site at Colorado State University Statistics Department. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center

Re: [R] how to calculate the quantile of the weighted data

2006-03-13 Thread Brian S Cade
and nonparametric smoothing models also are available. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL PROTECTED] tel: 970 226-9326 Frank E Harrell Jr [EMAIL PROTECTED] Sent by: [EMAIL PROTECTED

Re: [R] How to get around heteroscedasticity with non-linear least squares in R?

2006-02-21 Thread Brian S Cade
relation. Why not go after estimating the conditional quantile functions directly with nonlinear quantile regression, function nlrq() in the quantreg package? Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email

[R] inverse prediction intervals for nonlinear least squares

2006-01-06 Thread Brian S Cade
in the nonlinear model. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL PROTECTED] tel: 970 226-9326 [[alternative HTML version deleted]] __ R-help

Re: [R] Broken links on CRAN

2005-12-05 Thread Brian S Cade
the file to my hard disk and then click on it Adobe has no problem opening it. I'm using Adobe Reader 7.0 and Adobe Acrobat Standard 7.0. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL PROTECTED

[R] modifying code in contributed libraries - changes from versions 1.* to 2.*

2005-11-22 Thread Brian S Cade
documentation, etc. didn't suggest any simple approach like I use to do. Any help would be welcome. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL PROTECTED] tel: 970 226-9326 [[alternative

Re: [R] subsetting with by() or other function??

2005-10-13 Thread Brian S Cade
taylr-74 B 1997 3.63570636 0.48153827 5690 taylr-74 B 1998 1.94630194 3.63570636 Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL PROTECTED] tel: 970 226-9326 Florence

Re: [R] subsetting with by() or other function??

2005-10-13 Thread Brian S Cade
does the transformed tmp only have 1089 rows instead of 8314 like the full data frame? Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL PROTECTED] tel: 970 226-9326 Dimitris Rizopoulos [EMAIL

[R] subsetting data frame using by() or tapply() or other

2005-10-13 Thread Brian S Cade
3.88410354 5687 taylr-74 B 1994 3.29861501 3.32472557 5688 taylr-74 B 1996 0.48153827 3.29861501 5689 taylr-74 B 1997 3.63570636 0.48153827 5690 taylr-74 B 1998 1.94630194 3.63570636 Brian Brian S. Cade U. S. Geological Survey Fort

Re: [R] subsetting data frame using by() or tapply() or other

2005-10-13 Thread Brian S Cade
My thanks to Marc Schwartz who provided the solution - put unlist() around the tapply() statement. Looks like it works. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL PROTECTED] tel: 970 226-9326

[R] subsetting with by() or other function??

2005-10-12 Thread Brian S Cade
(dbc), i.e, shift them down by 1 row replacing the first row with NA, within groups defined by locid nested within population. I thought I could do something using by(csss3,list(locid, population), function) but don't seem to be having any success. Any suggestions?? Brian Brian S. Cade U. S

[R] stopping a function

2005-02-23 Thread Brian S Cade
are not executing properly and it would be nice to stop them without killing the R session. I've been using taskmgr to end the R session but I then lose all objects created during the session. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO

Re: [R] Fitting nonlinear (quantile) models to linear data.

2004-02-04 Thread Brian S Cade
? Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: [EMAIL PROTECTED] tel: 970 226-9326