[R] help with aggregate(): tables of means for terms in an mlm

(term, :) # browser() means - aggregate(resp, factors, mean) # rownames(means) - ... # means - means[, -(1:length(factors)] } termMeans.mlm(soils.mod, Contour) Error in FUN(X[[1L]], ...) : arguments must have same length thanks, -Michael -- Michael Friendly Email: friendly

Re: [R] mahalanobis

misunderstood something in the use of the function mahalanobis? Thanks. -- Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3

[R] R implementations of scatterplot/map labeling algorithims?

-09-27 Author - Elliot Noma DOI - 10.1007/BF02294366 Link - http://www.springerlink.com/content/c4k6205r83156165 Does anyone know of a related R (or other) public implementation of a solution to this problem? thanks, -Michael -- Michael Friendly Email: friendly AT yorku DOT

[R] example() and example(package=) like data() and data(package=) ?

there should not be parallel uses, example() and example(package=foo) not requiring a topic to be specified? That would be very handy, especially for exploring a new package or developing tools for documenting them, with output. -Michael -- Michael Friendly Email: friendly AT yorku DOT ca

Re: [R] means over factors in mlm terms

and provide commented, minimal, self-contained, reproducible code. -- Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3

Re: [R] Repeated Measures MANOVA in R

-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Michael Friendly Email: friendly AT yorku DOT ca Professor

[R] plot image matrix with row/col labels

, 35 Ouvriers, 36 Journaliers, 37 Domestiques, 38 Chevaux, 39 Chiens, 40 Moralite ) dimnames(mat) - list(rnames, cnames) # how to plot the image matrix with row/col names??? # show the image matrix plot(imagematrix(mat)) -- Michael Friendly Email: friendly AT yorku DOT ca

Re: [R] making matrix monotonous

@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax

[R] calculating similarity/distance among hierarchically classified items

to approach this. I'd appreciate any suggestions. -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3 CANADA

Re: [R] comparative density estimates

with the same bandwidth. Best wishes, Z On Thu, 23 Mar 2006 14:25:53 -0500 Michael Friendly wrote: I have two series of events over time and I want to construct a graph of the relative frequency/density of these events that allows their distributions to be sensibly compared. The events

[R] comparative density estimates

= argument to select a BW/grayscale plot. thanks, -Michael [1] http://www.math.yorku.ca/SCS/Gallery/milestone/ -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp

Re: [R] R newbie example code question

-- __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- Michael Friendly Email: friendly AT yorku DOT ca Professor

Re: [R] attributes of a data.frame

one right now. So that's another reason why it's not in the base: it doesn't need to be, you can just go find and install that contributed package! -- Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416

[R] loess: choose span to minimize AIC?

, gcv=gcv) return(result) } cars.lo - loess(dist ~ speed, cars) (values - loess.aic(cars.lo)) $span [1] 0.75 $aicc [1] 6.93678 $aicc1 [1] 167.7267 $gcv [1] 5.275487 -- Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. York University Voice

Re: [R] loess: choose span to minimize AIC?

-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Michael Friendly Sent: Thursday, November 17, 2005 9:58 AM To: R-help@stat.math.ethz.ch Subject: [R] loess: choose span

[R] visualizing/summarizing a large, sparse logistic regression

, or a more complex model involving Agency and Office plus potential interactions among A-H. Can someone point me in some useful directions or to some similar examples? thanks, -Michael -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416

Re: [R] plotting on a reverse log scale

-- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3 CANADA __ R-help

[R] plotting on a reverse log scale

the major ones, say at 100 year intervals. thanks, -Michael -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3

Re: [R] plotting on a reverse log scale

: On 7/6/2005 3:36 PM, Michael Friendly wrote: I'd like to do some plots of historical event data on a reverse log scale, started, say at the year 2000 and going backwards in time, with tick marks spaced according to log(2000-year). For example, see: http://euclid.psych.yorku.ca/SCS/Gallery

[R] abbreviate or wrap dimname labels

or abbreviate.or.wrap()? thanks, -Michael -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3 CANADA

Re: [R] abbreviate or wrap dimname labels

12:30 PM To: Michael Friendly Cc: R-Help Subject: Re: [R] abbreviate or wrap dimname labels On Fri, 2005-04-15 at 12:12 -0400, Michael Friendly wrote: For a variety of displays (mosaicplots, barplots, ...) one often wants to either abbreviate or wrap long labels, particularly when

Re: [R] Two-factorial Huynh-Feldt-Test

Bela Bauer wrote: I'm still looking for an efficient way to print the new summary. Is there any easy way to tell the summary or print functions about the corrected degrees of freedom? The standard way is not to print the adjusted df, but rather just the adjusted p-values (along with the

[R] debian /etc/apt/sources.list for CRAN?

ones used instead. -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3 CANADA

[R] re: patterns of missing data: determining monotonicity

== == == === === == == == == == 0 0 0 002 2 3 3 4 -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS

[R] patterns of missing data: determining monotonicity

0 0 00111 3 8 0000 0 10101 3 == == == === === == == == == == 0000 0 22334 -- Michael Friendly Email: [EMAIL PROTECTED] Professor

Re: [R] SAS or R software

, the active hard work of many R developers has led to R graphics for which the *default* results for many graphs avoid many of the egregious design errors introduced in SAS in the days of pen-plotters (+ signs for points, cross-hatching for area fill). -Michael -- Michael Friendly Email: [EMAIL

Re: [R] adjusting the map of France to 1830

/11/2004, Michael Friendly wrote: I'm doing some analyses of historical data from France in 1830 on 'moral statistics' that I'd like to show on a map. I've done most of my analyses in SAS, but a few things would work better in R. To do this, I have to adjust the modern map, library(maps) map

[R] adjusting the map of France to 1830

; proc gremove data=gfrtemp out=gfrance; by dept; id id; run; -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3

[R] read.delim problem with trailing spaces

. -Michael -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3 CANADA __ [EMAIL

[R] Lattice .ps graphic is rotated in LaTeX slides

, or how could I do it differently so it would work? thanks -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3 CANADA

Re: [R] Lattice .ps graphic is rotated in LaTeX slides

way: \rotatebox{180}{\includegraphics[angle=90,height=.6\textheight]{fig/barley2x3.ps}} thanks, -Michael Marc Schwartz wrote: On Fri, 2004-10-01 at 10:09, Michael Friendly wrote: I've generated a version of the classic dotplot of the barley data with library(lattice) data(barley) trellis.device

re: [R] boxplot notches

are equal. I believe this is explained in @Article{McGill-etal:78, author = R. McGill and J. W. Tukey and W. Larsen, year = 1978, title =Variations of Box Plots, journal = TAS, volume = 32, pages =12--16, } -- Michael Friendly Email: [EMAIL

[R] problem with abline for x.y

(x.y, col=gray, lwd=2) library(car) data.ellipse(child, parent, plot.points=FALSE, levels=c(0.40, 0.68), lty=2) The resulting figure may be seen at http://euclid.psych.yorku.ca/SCS/Gallery/Private/galton.jpg -Michael -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York

RE: [R] ANOVA procedure on the sufficient statistics

means and MSE, however. Larsen's procedure generates a weighted data set that gives the same ANOVA table as the raw data, but requires that the weight= argument be used in aov(). (weights are not supported in all related functions) -Michael -- Michael Friendly Email: [EMAIL PROTECTED

re: [R] ANOVA procedure on the sufficient statistics

. -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3 CANADA __ [EMAIL

[R] manova() with a data frame

also help if model.frame returned something to indicate which variable had an invalid type.) -Michael -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS

Re: [R] AGREP

be calculated using reduction operators, or even a correlation, based on the positions of the 1s. Thinking in arrays always helps avoid explicit loops. -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele

[R] nedit syntax highlighting patterns for R?

: == NEdit: style name required in style specification: Note:darkRed:Italic == NEdit: expecting quoted string in language mode specification: ==q .r .R:: -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814

re: [R] how to overlap plots

Try the distplot() function in package vcd. -Michael -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3 CANADA

Re: [R] correlation and causality examples

] Dear All, I'am looking for examples showing that correlation does not imply causality, the targeted audience consists of undergraduate students (their first year at the university but in the BioMathStat track). All practicals are under R. -- Michael Friendly Email: [EMAIL PROTECTED] Professor

Re: [R] Beginner: Homogenity of Variances

summary. This is implemented and described in my STAT2DAT macro, http://www.math.yorku.ca/SCS/sasmac/stat2dat.html It should be easy to translate into R. -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814

[R] ungrouping grouped data

this in R, but surely there must be an easier way. (In APL2 it would be frequency\[1]galton) thanks, -Michael -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca

[R] Minard's Challenge: Re-Visioning Minard Contest

with reactions to this challenge and what they deem useful to share with all readers. Other ways to reply include posting a web URL where readers can view the details or a direct email reply to me. -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice

Re: [R] Graphing

How about :-) system(sas multitot) where multitot.sas is goptions device=gif; proc gchart; hbar devices; by year; -Michael -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Street