similar to: Multiple comparisons following nlme

Hello, I would like to display two groups of dots in different colors or style and additionnaly a linear regression for all the dots in panel plots of a Trellis graph. Actually to get in each panel the equivalent of: plot(x[mois==3],y[mois==3],col="blue") points(x[mois==9],y[mois==9],col="red") abline(lm(y~x), col="green") "mois" being a grouping

Shalini Raghavan 3M Pharmaceuticals Research Building 270-03-A-10, 3M Center St. Paul, MN 55144 E-mail: sraghavan at mmm.com Tel: 651-736-2575 Fax: 651-733-5096 ----- Forwarded by Shalini Raghavan/US-Corporate/3M/US on 08/16/2004 11:25 AM ----- Shalini

Hi R-people, I am wanting to run Factorial ANOVA followed by Scheffe tests on some spatial subjective data. I'm comparing X-Y independent coordinates against x-y dependent coordinates. There are only four independent spatial coordinates that form a square. I am wondering whether I am doing the right thing, because there doesn't seem to be a simple way of doing this. I have attempted to

Hello R-Users, I have some data sets which change on a daily bases. So far I have imported these sets into R, done all my evaluations resulting in a couple of plots, charts and tables of numbers which I copy&pasted via clipboard into Powerpoint. The procedure is always the same and I wonder, whether there is no easier way for doing so. Is there some type of "report generator"

Hello! My graphs are produced using the "postscript"-option in R (R version 2.10.1 (2009-12-14)). When Greek letters are used on the axis, everything looks fine in the *.ps-file. If included in a LaTeX-file and (on Ubuntu 10.04, fresh install), the Greek letters appear in the DVI- and PS-output, however, if converted with ps2pdf they suddenly disappear. Could anyone suggest a

Dear mailing list readers! Using R and the plot function I stumbled over this little issue: plot(x=NULL, y=NULL, xlim=range(1:10), ylim=range(1e-9:1), log="y") produces an plot empty plot, where the y-axis is in the specified range. Changing ylim to ylim=range(1e-9:1e-3) creates an y-axis in the range of 1e-12 to 1e-6. This appearance of the y-axis is a

Hi, while executing simple sieve script: require ["fileinto", "envelope"]; if header :contains "subject" "test" { redirect "user at domain"; } we got: 2010-09-23T09:25:23+02:00 st1/192.168.10.201 dovecot: [ID 583609 mail.error] master: Error: service(lmtp): child 15524 killed with signal 11 (core not dumped) in case user home directory is

Hello, in the process of setting up a minimal xen host on debian wheezy, the call of xm new resulted in the following output Unexpected error: <type 'exceptions.ImportError'> Please report to xen-devel at lists.xensource.com Traceback (most recent call last): File "/usr/lib/xen-4.1/bin/xm", line 8, in <module> main.main(sys.argv) File

While I was helping a SAS-using friend with an analysis I noticed some differences in the multivariate test statistics, approximate F statistics, and p-values in the manova function using R and proc GLM using SAS. The univariate coefficients are identical. Is there a reason to expect R and SAS to give different results? Thanks, Bill Kristan.

Hello, my data.frame is sort of a collection of process values, i.e. huge run-chart. It consists of a time-stamp in the first column (date as string), factors in the following columns (used for subset-filtering), and some process-data columns. Hereafter, two examples are listed, showing the problems that occour during print: At first the example, that works fine:

Hello, I would like to perform a multivariate regression analysis to model the relationship between m responses Y1, ... Ym and a single set of predictor variables X1, ..., Xr. Each response is assumed to follow its own regression model, and the error terms in each model can be correlated. Based on my readings of the R help archives and R documentation, the function lm() should be able to

Hi, I believe the following is a very minor typo in the documentation for merge() : all logical; all=L is shorthand for all.x=L and all.y=L ^ ^ ^ Looks like 'L' should be 'FALSE'. Someone please alert me if I am wrong. If I am not wrong, is it appropriate to submit a bug for such minor detail? Thanks, Jim James A. Rogers

May I ask to add Return-Path: some meaningful header line to sieve sent mail headers in vacation message? Now the header line isn't generated at all and the effect is as follows: on receiving MTA: 2014-05-09T15:04:32+02:00 host/ip postfix/qmgr[2408]: 41F2F6024E: from=<>, size=900, nrcpt=1 (queue active) in received mail body: Return-Path: <> Would be nice to get: Return-Path:

As a long term user of qpopper with bulletins, I'd like to know if the bulletin feature is planned for dovecot? I really miss bulletins :-), this is very convenient way of providing information to users. Best regards, MU

This is a dovecot sieve log line: 2012-09-26T13:47:52+02:00 host/192.168.1.1 dovecot: [ID 583609 mail.info] lmtp(7565, user at domain): vTtJDTrrYlCNHQAAs1WjhA: sieve: msgid=<MSGID>: stored mail into mailbox 'INBOX' Would it be possible to have mail size logged here? MU -------------- next part -------------- A non-text attachment was scrubbed... Name: smime.p7s Type:

May I ask to add receiving smtp server mail queue id to sieve 'sent vacation response' log line? Now: 2014-05-09T15:04:32+02:00 name/ip dovecot: lmtp(8164, user1 at domain): YvJFLVbSbFPkHwAASpgOhg: sieve: msgid=<000301cf6b87$36be4570$a43ad050$@us.edu.pl> size=46850: sent vacation response to <user2 at domain> Asked for something like this: 2014-05-09T15:04:32+02:00 name/ip

Dear Colleagues, I have scoured the help files and been unable to find an answer to my question. Please forgive me if I have missed something obvious. I have run the following two models, where "category" has 3 levels and "comp" has 8 levels: mod1 <- lmer(x~category+comp+(1|id),data=impchiefsrm) mod2 <- lmer(x~category+comp+category*comp+(1|id),data=impchiefsrm)

Dear R community, I performed a generalized linear mixed model using glmmPQL (MASS library) to analyse my data i.e : y is the response with a poisson distribution, t and Trait are the independent variables which are continuous and categorical (3 categories C, M and F) respectively, ind is the random variable. mydata<-glmmPQL(y~t+Trait,random=~1|ind,family=poisson,data=tab) Do you think it

Hi all, I''m tring to isolate P2P traffic, specifically BitTorrent, for my QoS scripts. I can''t seem to completely isolate ALL BitTorrent traffic. I identify & mark packets and then use tc filters to put them into appropriate classes. My firewall rules (below) do the markings. My VoIP boxes'' and ICMP traffic get highest priority (mark 1). Then comes DNS, SSH,

Macros to use VMI calls from assembly and C languages are introduced. The macros are quite complex, but the end result is rather impressive. The result is that when compiling a VMI kernel, the native code is emitted inline, with no function call overhead, and some wiggle room for register allocation. The hypervisor compatibility code is emitted out of line into a separate section, and patched