[R] iplots Java

that the latest iplots R package is correctly installed. I would appreeciate hearing of strategies for solving the problem. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room

Re: [R] Artifacts in pdf() of image() (w/o comments)

Acrobat, so the problems aren't limited to Preview. Best wishes, Mark On 13/08/07, Duncan Murdoch [EMAIL PROTECTED] wrote: Michael Kubovy wrote: On Aug 12, 2007, at 6:24 AM, Duncan Murdoch wrote: Michael Kubovy wrote: Dear r-helpers, In my previous message there were comments

Re: [R] Artifacts in pdf() of image() (w/o comments)

But is it a bug? Can a program anti-alias text and line drawings and not bitmaps? On Aug 13, 2007, at 9:30 AM, Duncan Murdoch wrote: On 8/13/2007 11:07 AM, Michael Kubovy wrote: Dear Friends, Thanks for your input. FYI: Preview doesn't show PDF aliasing in the image I produced if I

Re: [R] Artifacts in pdf() of image() (w/o comments)

things that bother me that are not bugs.) Moreover, I'm not opposed to complaining to Apple, once I have been assured that I'm not reporting a bug where there's none. On Aug 13, 2007, at 9:56 AM, Duncan Murdoch wrote: On 8/13/2007 11:43 AM, Michael Kubovy wrote: But is it a bug? Can

Re: [R] Artifacts in pdf() of image() (w/o comments)

On Aug 12, 2007, at 6:24 AM, Duncan Murdoch wrote: Michael Kubovy wrote: Dear r-helpers, In my previous message there were comments in the code that may have made cutting and pasting awkward. Here it is w/o them. I have two questions: (1) The following produces a pdf with artifacts

[R] Artifacts in pdf() of image()

= '', bty = 'n', axes = FALSE, asp = 1) dev.off() (2) I would like the output to be clipped to a circle, i.e., anything outside the circle tangent to the sides of the square should be transparent. How can I do that? _ Professor Michael Kubovy University

[R] Artifacts in pdf() of image() (w/o comments)

= '', ylab = '', bty = 'n', axes = FALSE, asp = 1) dev.off() (2) I would like the output to be clipped to a circle, i.e., anything outside the circle tangent to the sides of the square should be transparent. How can I do that? _ Professor Michael Kubovy University of Virginia

Re: [R] Standard errors of the predicted values from a lme (or lmer)-object

for the coefficients. Otherwise, you can do some calculations with residuals(lmeresult). Most useful for diagnostic purposes is plot (lmeresult). Perhaps ?estimable in the gmodels package _ Professor Michael Kubovy University of Virginia Department of Psychology USPS

Re: [R] title line

introductions to R available on line? Also, you're not asking your questions well (this comment is not about knowing English well): You would have been better off asking: I have a data frame df - data.frame(x1 = rnorm(3), x2 = rnorm(3)) How do I ? _ Professor Michael

Re: [R] Selecting complementary colours

the posting guide http://www.R-project.org/posting- guide.html and provide commented, minimal, self-contained, reproducible code. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels

[R] Separate x-axis ticks for panels in xyplot()?

'), rot = 45}, {at = v2, labels = c('E', 'F' 'G), rot = 45} ) ) ) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box

[R] Example of mcmcsamp() failing with lmer() output

it? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab

Re: [R] xyplot() and controlling panel.polygon()

', ...) panel.xyplot(x, y, ...) } ) It's pretty clear that panel.superpose is not getting its x and y values after they are split by group and panel. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville

Re: [R] xyplot() and controlling panel.polygon()

. I've been trying to get it from the help pages and from Murrell's book, but neither offers enough of a cookbook for me to figure these things out. Is there something I should have read? Thanks, MIchael On Apr 27, 2007, at 2:25 PM, Deepayan Sarkar wrote: On 4/27/07, Michael Kubovy [EMAIL

Re: [R] xyplot() and controlling panel.polygon()

I can't thank you enough! (I am also grateful to Sundar Dorai-Raj and Bert Gunter for their initial responses.) On Apr 27, 2007, at 4:14 PM, Deepayan Sarkar wrote: On 4/27/07, Michael Kubovy [EMAIL PROTECTED] wrote: Hi Deepayan, Your solution works, anf the polygon are drawn where I wanted

[R] xyplot() and controlling panel.polygon()

, ...) llines(x = c(1, 6), y = c(-24.28, 35.941667), lwd = 2, lty = 3, col = 4) } x[c, g] and y[c, g] describe the polygon I want plotted for group g in panel c. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS

Re: [R] xyplot() and controlling panel.polygon()

Hi Sundar, Thanks for taking a look at this On Apr 26, 2007, at 1:05 PM, Sundar Dorai-Raj wrote: Michael Kubovy said the following on 4/26/2007 7:20 AM: Dear R-helpers, How do I tell panel.polygon what greoup and panel it applies to whithin xyplot(y ~ x | c, groups = g panel

Re: [R] xyplot() and controlling panel.polygon()

Hi Sundar, On Apr 26, 2007, at 10:03 PM, Deepayan Sarkar wrote: On 4/26/07, Michael Kubovy [EMAIL PROTECTED] wrote: Hi Sundar, Thanks for taking a look at this On Apr 26, 2007, at 1:05 PM, Sundar Dorai-Raj wrote: Michael Kubovy said the following on 4/26/2007 7:20 AM: Dear R

[R] Bug in xYplot() with method = 'filled bands'

, Michael Kubovy __ 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.

Re: [R] Positioning in xyplot

, line.col = 'black', side = side, ...) } }) I also wonder why one doesn't need to call panel.xyplot(). Thanks, MK _ Professor Michael Kubovy University of Virginia Department

[R] Changing appearance of band edges in xYplot()

. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab:B019+1-434-982-4751

Re: [R] Changing appearance of band edges in xYplot()

On Apr 21, 2007, at 11:08 AM, Michael Kubovy wrote: Dear R-helpers, I would like the bands in an xYplot() to be gray without dark edges. I've tried playing around with col.bands, but that seems not to affect the edges, although changin col.fill does change the fill as I would expect. I

Re: [R] Positioning in xyplot

Wow! Thanks to Sundar and to Deepayan for this selfless effort.(I rather suspect they had a good time figuring this out.) The outcome is that I will be able to publish just the right figure for the job. I will be sure to aknowledge your contribtution. Michael On Apr 12, 2007, at 12:48 AM,

[R] Positioning in xyplot

= ' = ', shingle.intervals = NULL), data = df, type = 'b', label.curves = F, layout = c(2, 3), ) Currently this puts the panels as follows 3 2 4 1 5 I need: 3 2 4 1 5 How can I do this? Any help will be much appreciated. _ Professor Michael Kubovy

Re: [R] Positioning in xyplot

On Apr 10, 2007, at 7:07 PM, Deepayan Sarkar wrote: On 4/10/07, Michael Kubovy [EMAIL PROTECTED] wrote: Dear R-helpers, I have an xyplot of the following kind: xYplot(y ~ x | p, groups = factor(gg, levels = c('1', '5', '2', '4', '3')), strip = strip.custom(var.name = 'g

[R] Failure of mcsamp() but not mcmcsamp()

gridBase MASS JGR iplots JavaGDrJava 0.4-3 7.2-33 1.4-15 1.0-5 0.3-6 0.4-14 Thanks, MK _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

[R] Why does lmList() fail when lm() doesn't?

6.8269 3.4134 1.9976 0.2059 # Residuals 7 11.9612 1.7087 _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick Road

Re: [R] partial R

to proceed is install.packages('car') ?cr.plots (You may have to restart R to get the help on a newly-installed package.) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels

Re: [R] Doing partial-f test for stepwise regression

not matter? You can tell which is which by looking at the degrees of freedom. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

Re: [R] Newbie: Combn and scripting

() for each line into another file. How do I do this? data.csv: a,b,c,d 1,2,3.4 g,3,6,t etc x=data.csv, m=3 _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room

Re: [R] distribution graph

?violinplot (You need to install the UsingR package first.) On Mar 23, 2007, at 4:06 AM, Plessen, Christian von wrote: I am looking for a way to produce a distribution graph as in the example: (http://cecsweb.dartmouth.edu/release1.1/datatools/dgraph.php?

Re: [R] package:AlgDesign and .Random.seed

On Mar 21, 2007, at 4:16 AM, Uwe Ligges wrote: Michael Kubovy wrote: Dear r-helpers, Could you please help me solve the following problem: When I run require(AlgDesign) trt - LETTERS[1:5] blk - 10 trtblk - 3 BIB - optBlock(~., withinData = trt, blocksizes = rep(trtblk, blk)) In response

[R] package:AlgDesign and .Random.seed

0.1-40.14-16 0.4-3 7.2-32 1.4-15 1.0-5 JavaGDrJava 0.3-6 0.4-14 _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels

[R] Help with xYplot() in package Hmisc

- y + e ll - y - 1 ul - y + 1 g - c(rep(1, 4), rep(2, 4), rep(3, 4)) df - data.frame(x, y, ll, ul, g) xYplot(Cbind(y, ll, ul) ~ x, groups = g, type = 'p', data = df) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

Re: [R] logging mouse clicks

Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab:B019+1-434-982-4751 Fax:+1-434

[R] Problem with ci.lmer() in package:gmodels

7.056447 0 Days 10.52773 13.63372 7.389946 1.646900 0 _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

[R] Is the gmodels package being maintained?

and/or address with a current maintainer and/or provide a current email address? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick

Re: [R] Is the gmodels package being maintained?

Hi, Finding his email address was not immediate, but I finally did, and did bring the problem to Greg's attention @ rochester, and the message didn't bounce this time. On Mar 9, 2007, at 7:43 AM, Peter Dalgaard wrote: Michael Kubovy wrote: Dear r-helpers, I sent a cc of a recent

Re: [R] Extracting the p of F statistics from lm

. try Sfitting[4]$coefficients[,4] I'm not sure that this is the best way, but it works with the example for lm() summary(lm.D9)[4]$coefficients[,4] # (Intercept) groupTrt # 9.547128e-15 2.490232e-01 _ Professor Michael Kubovy University of Virginia Department

Re: [R] Table Construction from calculations

or www.ms.washington.edu/stat390/winter07/R_primer.pdf) which allows you do avoid writing loops. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

[R] Failure to run mcsamp() in package arm

I would appreciate a pointer to what the problem might be. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick Road

Re: [R] Failure to run mcsamp() in package arm

More problems. If I run sim(fm1 - lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) from the lmer() help page. I get the error Error in mvrnorm(n.sims, bhat[j, ], V.beta) : 'Sigma' is not positive definite On Mar 7, 2007, at 1:30 PM, Michael Kubovy wrote: Dear r-helpers, I can run

Re: [R] Failure to run mcsamp() in package arm

arm R2WinBUGSlme4 Matrix lattice 0.8-18 1.2-11.0-13 2.0-4 0.9975-13 0.9975-11 0.14-16 MASS JGR iplots JavaGD rJava 7.2-321.4-15 1.0-5 0.3-50.4-14 On Mar 7, 2007, at 4:30 PM, Michael Kubovy wrote: More

Re: [R] Question

, c(Var1 * 100, Var2 * 200)) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011

Re: [R] How to override ordering of panels in xyplot()

On Mar 3, 2007, at 3:19 PM, Deepayan Sarkar wrote: On 3/3/07, Michael Kubovy [EMAIL PROTECTED] wrote: Dear r-helpers, I'm conditioning an xyplot on a variable whose levels are'low', 'med', 'high'. How do I override the alphabetical ordering for the panels of the plot? This has less

[R] Help with paste()

Dear r-helpers, Could you please tell me what's missing: rbind(paste('txt.est',1:24, sep = '')) txt.est1, ... txt.est24 are vectors that I wish to rbind. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

[R] How to override ordering of panels in xyplot()

Dear r-helpers, I'm conditioning an xyplot on a variable whose levels are'low', 'med', 'high'. How do I override the alphabetical ordering for the panels of the plot? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box

Re: [R] mtext bold font problem

( Difference))) ) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982

Re: [R] question about boxplot

, and the character rotation. Would you consider an easy way out---an alternative with reasonable defaults? data(ToothGrowth) bwplot(dose ~ len | supp, ToothGrowth) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

[R] Help with xlab and ylab distance from axes

Dear r-helpers, In basic graphics, I have a figure with x and y axes suppresed. I would like to move the xlab and the ylab closer to the axes. Do I have to use mtext()? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box

[R] Suppressing \newcolumntype{} declaration in latex.default() in Hmisc

Dear r-helpers, When use latex() on a matrix, I set the option dcolumn = T. As a result, in front of each tabular I get \newcolumntype{.}{D{.}{.}{-1}} The LaTeX compiler complains about these multiple redeclarations. _ Professor Michael Kubovy University of Virginia

Re: [R] Suppressing \newcolumntype{} declaration in latex.default() in Hmisc

these multiple redeclarations. Is there a way to suppress this? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick Road

Re: [R] Suppresing default text in pairs.lmList() in package = nlme

Thanks Ken and Dieter, I added xlab = '' and the text 'Scatter Plot Matrix' produced by pairs () applied to an lmList object (in nlme) went away. Skip the rest --- which is mainly autobiographical --- unless you're curious. On Feb 11, 2007, at 4:25 AM, Dieter Menne wrote: Michael Kubovy

[R] OT: question about transformation

If x is a r.v. 1/log(x) ~ N (approximately) what does that tell me about the distribution of x (i.e., is it a standard distribution?) Can you point me to a list where I would have appropriately posted this? _ Professor Michael Kubovy University of Virginia

[R] Suppresing default text in pairs.lmList() in package = nlme

I would like to suppress the text 'Scatter Plot Matrix' that appears under the plot. Could someone please suggest how? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels

Re: [R] How to count the number of NAs in each column of a df?

wonder if we have a place to keep R programming idioms (which probably get unnecessarily reinvented). Is the R-Wiki suitable? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels

Re: [R] Replace individual values in a data frame with NA

? a-c(1,2,3,4) b-c(5,6,7,8) data-data.frame(a,b) data[3, 2] - NA ?Extract _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick Road

[R] How to add the variable name to a qqplot or densityplot in the diagonal of an splom?

'), gp = gpar(fontsize = 10)) }, varname.col = 'blue', varname.cex = 1 ) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

[R] smartpred depends on fitted() in flexmix?

loading of the flexmix package. newx = seq(0, 1, len=n) points(newx, predict(fit, data.frame(x=newx)), type=b, col=2, err=-1) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels

[R] Strange behavior of abc.ci() in boot

JavaGD 1.2-27 1.4-3 0.1-40.14-16 0.4-3 7.2-31 1.4-15 1.0-5 0.3-5 rJava 0.4-13 _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400

[R] How to count the number of NAs in each column of a df?

I would like to remove columns of a df which have too many NAs. I think that summary() should give me the information, I just don't know how to access it. Advice? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

Re: [R] lmer2 error under Mac OS X on PowerPC G5 but not on Dual-Core Intel Xeon

)) Error in as.double(start) : Calloc could not allocate (903190944 of 4) memory On Jan 28, 2007, at 8:35 PM, Michael Kubovy wrote: (fm1 - lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)) Error in as.double(start) : Calloc could not allocate (888475968 of 4) memory

[R] lmer2 error under Mac OS X on PowerPC G5 but not on Dual-Core Intel Xeon

rJava 1.4-15 1.0-5 0.3-5 0.4-13 * lmer runs the example w/o a problem I just tried to run it on on Intel-based MacPro, and lmer2 ran without a hitch. _ Professor Michael Kubovy University of Virginia Department

[R] Unbalanced design help

discussion of the analysis of such data, and tools in R to use? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick Road

[R] Adding lines to xYplot

I am using xYplot to plot data and CIs. How do I add several lines to the figure? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

Re: [R] Finding the effect of Box-Cox transformation using vis.boxcoxu

any option to put user's own dataset. Can anyone tell me how to put my own dataset here i.e. data? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer

Re: [R] Newbie question: Statistical functions (e.g., mean, sd) in a transform statement?

NA 52.41755 NA 83.46864 55.64576 21.52078 6 54.01291 52.48922 53.97689 NA 160.47902 80.23951 46.33038 _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room

Re: [R] ANCOVA

On Jan 6, 2007, at 8:34 AM, John Cardinale wrote: Are there any R function which can do analysis of covariance? ?lm RSiteSearch('ancova') _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904

Re: [R] Help with filled.contour()

To Jeff: a tip of the hat. I have another question after Jeff's solution: On Jan 3, 2007, at 5:09 PM, Jeffrey Horner wrote: Michael Kubovy kubovy at virginia.edu writes: I tried and it gave a strange result. See http://people.virginia.edu/~mk9y/mySite/twoGaussian.R and http

Re: [R] New to R

by the more general http://cran.r-project.org/doc/manuals/R-intro.pdf _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick Road

[R] Help with filled.contour()

) *** It has four drawbacks. I would like: (1) to remove the white grid; (2) to remove the white contours; (3) its appearance to be smooth; (4) to remove the key. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

Re: [R] Help with filled.contour()

On Jan 1, 2007, at 2:13 PM, Dieter Menne wrote: Michael Kubovy kubovy at virginia.edu writes: The following plot is a first approximation to what I need: .. see below *** It has four drawbacks. I would like: (1) to remove the white grid; (2) to remove

Re: [R] slightly extended lm class

() in the nlme package? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1

Re: [R] Difference of array and vector

of attributes: http://www.burns-stat.com/pages/Spoetry/essentials.pdf There Burns writes (p. 3): An S array is merely a vector that has a dim attribute, and optionally a dimnames attribute. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS

Re: [R] sequential row selection in dataframe

- GAGurine[-seq(nr, len, 4),] # drop every 4th, giving 235 rows nr - 3 # or 4 will give 236 rows. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

Re: [R] Rating competitors

often. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab

Re: [R] complex barplot enquiry

unfit 3110 0 0 0 0 0 0 0 0 unsound 44185 3 1 0 0 0 0 0 1 Could you please post commands to define the df 'offend'? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

Re: [R] Newbie data organisation/structures question...

(p1t1 = subset(dat, p == 1 t ==1)$iv, p1t2 = subset(dat, p == 1 t ==2)$iv) # correcting your requirement to get your result There are many other ways of getting the results you need as Marc Schwartz pointed out in his reply. _ Professor Michael Kubovy University

Re: [R] how to replace some objects?

) num - recode(let, 'a' = 1; 'b' = 2; else = 3 ) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA

[R] Problem with glmmADMB

0.4-12 7.2-30 0.14-16 _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1

Re: [R] digital length

?Round On Dec 14, 2006, at 2:52 AM, XinMeng wrote: How can I control the digital length of data? e.g: 0.1234 is the output of an algorithm. What I want is 0.12 instead. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box

Re: [R] A smal fitting problem...

) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab:B019

[R] Sweave and warning messages

How does one tell Sweave() to include analysis warning messages in the verbatim output? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

Re: [R] stat question - not R question so ignore if not interested

a different approach: require(sm) xy.sm - sm.regression(xy$x, xy$y) Whenever there's a big discrepancy between an OLS fit and a robust one, we should not pursue the OLS one w/o reinterpretation, which others have discussed in their replies. _ Professor Michael Kubovy

Re: [R] Help with response CIs for lme

question is --- if the CIs I obtained were correct --- is there a less laborious way to obtain them? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

[R] overplot() examples fail (package:gplots)

_ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab:B019+1-434-982-4751 Fax

[R] Problem with pairs() in nlme

- subset(vfr, slant %in% c(10, 20, 30)) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903

[R] Slight discrepancy between predict.lm() and all.effects()

[[1]]$fit, lwr = hsb.eff[[1]]$lower, upr = hsb.eff[[1]]$upper ) Is this due to rounding or algorithm? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102

[R] Help with response CIs for lme

]] [2,], hsb.int[[1]][1,]+hsb.int[[1]][3,], hsb.int[[1]][1,]+hsb.int[[1]] [2,] + hsb.int[[1]][3,] + hsb.int[[1]][4,])) If so, is there an easier way to write it? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

[R] Help with plot(augPred()) plot(comparePred()) in nlme

_ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab:B019+1-434-982

Re: [R] have I an actual matrix?

, 558, 134, 53, 9) names(nums) - c(2002, 2002, 2002, 2002, 2002, 2003, 2003, 2003, 2003, 2004, 2004, 2004, 2004) barplot(nums) or require(gplots) barplot2(nums, plot.grid = TRUE, las = 1) # offers many nice options _ Professor Michael Kubovy University of Virginia

Re: [R] Equivalent to Stata command

Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab:B019+1-434-982-4751 Fax:+1-434

[R] CI -- fixed effects -- response -- mixed models

to you for guidance on how to proceed. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office

Re: [R] no possible to load a package correctly dowloaded by the R prompt

('partitions') You could have said install.packages('untb', dependencies = TRUE), and this would have done the job for you. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102

[R] Help with unlist

Dear r-helpers, I have a list whose elements are str(durCut[[1]]) Ord.factor w/ 5 levels vLowlowmed..: 3 2 5 2 2 2 4 4 3 5 ... How do I unlist durCut into an ordered factor? str(unlist(durCut)) int [1:3024] 3 2 5 2 2 2 4 4 3 5 ... _ Professor Michael Kubovy

[R] plot(all.effects)

= coln, nrow = 2, ncol = 2, more = last[i]) } Error in plot.effect(x[[(i - 1) * cols + j]], row = i, col = j, nrow = rows, : formal argument row matched by multiple actual arguments _ Professor Michael Kubovy University of Virginia Department of Psychology

[R] Constructing predictions from HPDinterval() after lmer()

on 5 and 129 DF, p-value: 2.2e-16 rd.eff - all.effects(rd) rd.ci - data.frame(y = rd.eff[[1]]$fit, Lower = rd.eff[[1]] $lower, Upper = rd.eff[[1]]$upper) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

Re: [R] plotting 95% confidence bands on a simple linear regression model from lm()

) ?summary.effect _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729

[R] Hide line ends behind unfilled circles?

(2, 3, 3, 4, 6, 4, 5, 5, 6, 7, 7, 7, 8, 8)] y1 - yy[c(2, 3, 3, 4, 6, 4, 5, 5, 6, 7, 7, 7, 8, 8)] plot(yy ~ xx, pch = aa, cex = 3) segments(x0, y0, x1, y1) Can anyone suggest a way of insuring that the lines are hidden behind the unfilled circles? _ Professor Michael

[R] Compact presentation of multiple figures

), mar = c(1.1, 2.1, 1.1, 0.1)) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office

  1   2   >