search for: flom

...I try model4 <- polr(ypsxcat~committed + as.factor(sex) + as.factor(drugusey) + anycsw + as.factor(sex)*committed + as.factor(sex)*as.factor(drugusey)+as.factor(sex)*anycsw, data = duhray) I get a message that the polr function was not found. Any help appreciated thanks Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...1 9.6 1 48.9 1 12.9 0 63.0 0 51.8 etc. for several hundred rows. I would like a new data frame that has all the cases where In_2006 = 1, and those cases of In_2006 that are closest to those cases Thanks in advance Peter Peter L. Flom, PhD Statistical Consultant Website: www DOT peterflomconsulting DOT com Writing; http://www.associatedcontent.com/user/582880/peter_flom.html Twitter: @peterflom

...ith error code 1 in: shell(cmd, wait = TRUE) 3: cmd execution failed with error code 1 in: shell(cmd, wait = TRUE) 'rm' is not recognized as an internal or external command, operable program or batch file. I have no idea what these mean.....any help appreciated thanks Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...finzi.psych.upenn.edu/R/doc/html/packages.html rmutil is in the standard library..... If I ignore the message and try fitting a model with one of the functions in gnlm (e.g. fmr) I get an error that it coudn't find the function 'finterp' Any help appreciated Thanks Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...onstructed What I want to do is simulate this with normally distributed data for the X's. That is, generate (say) 1000 sets of X1 through X4, use those to caluculate 1000 var1 and var2, and then 1000 correlations between var1 and var2, and then plot those results. Any help appreciated Peter Flom Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...45, 40) new <- c(28, 32, 35, 36, 41, 38, 47, 48, 52) what I'd like is to plot reelect and new (maybe with smoothing) on a time scale, with the x axis labeled for date..... But I can't figure out how, despite looking at ?ts and various things there Thanks in advance Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...search on the effects of rounding on regression? e.g., when you ask people "How often have you _______?" you are more likely to get answers like 100, 200, etc. than 98, 203, etc. I'm interested in investigating this, but don't want to reinvent the wheel. thanks Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...ne. Has anyone programmed the Welch test for equality of variances? I tried RSiteSearch, but this gave references to t test and oneway.test, which are not quite what I need.....I need the Welch test itself, for use in a meta-analysis (to determine if variances are equal). TIA Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St http://cduhr.ndri.org www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...s on the tick marks (that is, the numbers) without shrinking anything else. I could not find this in either the Rhelp archives or in ?dotchart, which recmmends cex to avoid 'label overlap', but cex shrinks all the characters in the plot. Is there a way to do this? Thanks Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...ggested in Pinheiro & Bates produce errors; e.g. plot(m1.mod1, schoolnum~resid(.), abline = 0) gives an error "could not find function "NaAct". Searching the archives showed a similar question from 2007, but did not show any responses. Thanks for any help Peter ) Peter L. Flom, PhD Statistical Consultant Website: www DOT peterflomconsulting DOT com Writing; http://www.associatedcontent.com/user/582880/peter_flom.html Twitter: @peterflom

...llinearity 2) VIFs don't diagnose the number of collinearities and 3) No one has determined how high a VIF has to be for the collinearity to be damaging. He then develops and suggests using condition indexes instead, so I was wondering if anyone had programmed them. Thanks Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...dotchart(mol, labels = prof, main = 'Dot chart', xlab = 'Meaning of life score') I get a dot chart sorted by the values of prof; what I'd like is one sorted by values of mol. I looked at the help file for dotchart, but did not see anything. Thanks again in advance Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...his didn't work. What's the best way to do this? More specifically 1) How do I import any .sd2 file to R ? 2) Is it better to create and import small files as needed, or one huge file which would contain all the data? 3) Any other tips on how to do this? Thanks in advance Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax) -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

...did not work. e.g. with(gdsgraph, gdsbox(BReT3T5T)) produced an error that gds3level was not found but if I try with(gdsgraph, boxplot(BReT3T5T~gds3lev)) it works fine. Similar problems occurred when I tried to use data = What am I missing? Thanks Peter Peter L. Flom, PhD Brainscope, Inc. 212 263 7863 (MTW) 917 488 7176 (ThF) [[alternative HTML version deleted]]

...DRUGUSEY>=0.5 63 65.428570 1.4285710 14) SEX>=1.5 21 6.571429 1.1428570 * 15) SEX< 1.5 42 56.285710 1.5714290 * So, is this the 'exception that proves the rule'? Have I done something wrong? Or what? Any ideas or thoughts? Thanks in advance Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

Proposal That a new mailing list be established that pertains exclusively to R documentation. The purpose of the list would be to discuss weak sections of the documentation and establish fixes for those weak spots. Pro If it works, there would be better documentation. It would be an excellent opportunity for newish and/or less technical people to contribute to R. In some respects such people

...1 on a Windows machine without installing MASS and without getting any errors, but I did. Third, thanks to all who replied, and especially to Andy Jaworski who showed me how to solve the problem, namely to get the package VR from CRAN and then type library(MASS). Thanks again, Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

In The American Statistician vol 42 (1988) pages 257 - 280, Yoav Benjamini investigates some variations on the box plot, including vaseplots, which maek the width of each box vary proportionally to he estimated density at a particular point. Has anyone implemented these in R ? Thanks as always Peter

...eated with lm..... Relevant code <<<<< model.borrow.logols<- lm(logborrow~age + sex + racgp + yrseduc + needlchg + gallery + totni + inject + poly(year.of.int,3) + druginj + inj.years + HTLV3) fastbw(model.borrow.logols) >>>> Thanks in advance Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research Institutes 71 W. 23rd St www.peterflom.com New York, NY 10010 (212) 845-4485 (voice) (917) 438-0894 (fax)

...rows but variable(s) found have 796 rows >>> However, when I look at 'newdata' (i.e. traintest1) with dim, I find it has 796 rows: <<< > dim(traintest1) [1] 796 612 >>> So, I do not understand the warning. Thanks for any explanation or assistance Peter L. Flom, PhD Brainscope, Inc. 212 263 7863 (MTW) 212 845 4485 (Th) 917 488 7176 (F) [[alternative HTML version deleted]]