[R] plotting predicted curves with log scale in lattice

that the x-variable of llines must be logged to plot the correct values and so the scales argument seems to apply only to the x, y arguments passed to the panel function. Thank you. best, Ken -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives

Re: [R] plotting predicted curves with log scale in lattice

Excuse me for forgetting sessionInfo (below) Ken Knoblauch knoblauch at lyon.inserm.fr writes: I was taken off guard by the following behavior in a lattice plot. I frequently want to add a predicted curve defined at more points than in the formula expression of xyplot. There have been

[R] Sweave rendering of simple character

but I can't seem to make this work. I have also tried the results=verbatim argument. Thanks, in advance for any suggestions. best, ken -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0

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

/2005-March/001650.html https://stat.ethz.ch/pipermail/r-sig-mac/2006-February/002678.html Simon Urbanek suggested one patch on the later thread. My solution, at the time, was to capture the image from the acrobat rendition. best, ken -- Ken Knoblauch Inserm U846 Institut Cellule Souche et

Re: [R] fit a nonlinear model using nlm()

William Simpson william.a.simpson at gmail.com writes: I am trying to fit a nonlinear model using nlm(). The observer is trying to detect a signal corrupted by noise. On each trial, the observer gets stim=signal+rnorm(). In the simulation below I have 500 trials. Each row of stim is a new

Re: [R] A More efficient method?

Keith Alan Chamberlain Keith.Chamberlain at Colorado.EDU writes: Cat=c('a','a','a','b','b','b','a','a','b')# Categorical variable C1=vector(length=length(Cat)) # New vector for numeric values for(i in 1:length(C1)){ if(Cat[i]=='a') C1[i]=-1 else C1[i]=1 } C1 [1] -1 -1 -1 1 1

Re: [R] if statement

- cut(x, quantile(x, c(0, 0.3, 0.7, 1)), c(0, 2, 1), TRUE) ken -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http

Re: [R] Generating artificial datasets with a specific correlationcoefficient.

see mvrnorm in MASS and especially the empirical argument James Milks james.milks at wright.edu writes: I need to create artificial datasets with specific correlation coefficients (i.e. a dataset that returns r = 0.30, etc.) as examples for a lab I am teaching this summer. Is there a

Re: [R] Subset and logical operator error

Sébastien pomchip at free.fr writes: Can you please point to me my syntax mistake or indicate a method to get this type of data.frame subset ? Thank you in advance ID value 1 1 1.2 2 2 1.2 3 3 1.2 4 4 1.2 5 5 A 6 6 A 7 7 A 8 8 A

[R] [R-pkgs] new packages psyphy and MLDS

, bug-reports, etc. are always welcome. Best, Ken Knoblauch -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http

[R] runif with weights

out but the probability of each number from 11 to 100 to come out is 60%. -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10

[R] Selecting complementary colours

+1-434-982-4729 Lab:B019+1-434-982-4751 Fax:+1-434-982-4766 WWW:http://www.people.virginia.edu/~mk9y/ -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91

[R] A function for raising a matrix to a power?

A%*%A%*%A, is there an abbreviation similar to A^3? Atte Tenkanen A=rbind(c(1,1),c(-1,-2)) A [,1] [,2] [1,]11 [2,] -1 -2 A^3 [,1] [,2] [1,]11 [2,] -1 -8 But: A%*%A%*%A [,1] [,2] [1,]12 [2,] -2 -5 -- Ken Knoblauch Inserm U846 Institut Cellule Souche et

[R] How to solve difficult equations?

I don't see the problem, except that you might want to think about what the error message is telling you. A little exploration of your function always helps, too. ss - seq(-2, 2, len = 100) plot(ss, fn(ss), type = l) uniroot(fn, c(-1, 1)) Erreur dans uniroot(fn, c(-1, 1)) : f() values at

[R] problem with mfg argument of par

)) for (ix in 4:1) { par(mfg = pos[ix, ]) plot(1:5) } Thank you in advance. R version 2.4.1 Patched (2007-01-23 r40561) i386-apple-darwin8.8.1 locale: C attached base packages: [1] stats graphics grDevices utils datasets methods [7] base -- Ken Knoblauch

[R] using scan to record user's input

: MASS 7.2-31 -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http://www.pizzerialesgemeaux.com/u846

Re: [R] using scan to record user's input

are. On 3/11/07, Ken Knoblauch [EMAIL PROTECTED] wrote: I'm using scan in a script to record a series of responses of the user as a function of some graphs that I put up on the screen. A toy version would be, y - rep(NA, 3) for (ix in seq( length(y) ) ) { y[ix] - scan( n = 1

Re: [R] using scan to record user's input

7.2-30 2.22-19 Otherwise, I'm as stumped as you are. On 3/11/07, Ken Knoblauch [EMAIL PROTECTED] wrote: I'm using scan in a script to record a series of responses of the user as a function of some graphs that I put up on the screen. A toy version would be, y - rep(NA, 3) for (ix

Re: [R] using scan to record user's input

starts running after the last brace and gets as its input the next line sent from the console. Either one of the solutions you propose (using source or putting it in a function) should fix your problem. Best, Pat On 3/11/07, Ken Knoblauch [EMAIL PROTECTED] wrote: I get it on a powerpc Mac

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

how? Must be as special Virginia Brand of lmList. That Test does not turn up in my output, and the only place I found it in the sources was as a comment ## scatter plot matrix plots, generally based on coef or ranef Dieter -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau

[R] error using user-defined link function with mixed models (LMER)

difference appears to be the valideta parameter (it's NULL in the logexposure family). -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33

Re: [R] error using user-defined link function with mixed models (LMER)

of the binomial and logexposure families, the big difference appears to be the valideta parameter (it's NULL in the logexposure family). -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel

[R] enumerates all possible combinations

does expand.grid do what you want? expand.grid(c(0, 1), c(0, 1), c(0, 1)) Var1 Var2 Var3 1000 2100 3010 4110 5001 6101 7011 8111 Hi all R users, I want to create a matrix having n columns and

[R] Regression lines

Try this version of your function and then think about it tst - function () { attach (attitude) x - rating y - learning detach (attitude) plot (x, y) abline(v=mean(x)) abline(h=mean(y)) abline (lm(y~x)) cc - coef(lm(x ~ y)) abline (-cc[1]/cc[2], 1/cc[2]) } My simpleminded understanding of

[R] Regression lines

This should do the trick: mind_reader - function() { ll - letters[round(runif(6, 1, 26))] ff - ll[1] for (ix in 2:length(ll)) { ff - paste(ff, ll[ix], sep = ) } if (exists(ff)) { cat(The function that you were

[R] how to adjust link function in logistic regression to predict the proportion of correct responses in 2AFC task?

272595 -- Ken Knoblauch Inserm U371 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http://www.lyon.inserm.fr/371/ [[alternative HTML version

[R] Vertical line in densityplot?

I think that you are mixing lattice and base graphics. This works for me: library(lattice) x-rnorm(100) densityplot(x, panel = function(x, ...) { panel.densityplot(x, ...)

Re: [R] deriv when one term is indexed

-1.612492e-01 2.988391e-02 beta1 1.269000e-05 1.792914e-05 beta2 3.844041e-05 5.388546e-05 beta3 6.108281e-06 8.762678e-06 gamm 2.481102e+00 2.542966e+00 R.version.string # XP [1] R version 2.4.0 Patched (2006-10-24 r39722) On 11/18/06, Ken Knoblauch [EMAIL PROTECTED] wrote

Re: [R] deriv when one term is indexed

sqrt(F) = 0.0068 Apologies for the error, to the authors of MASS, who wrote this function and within whose package it can be found. ken Ken Knoblauch wrote: I have found a way to approach the question that I posed, that is having deriv and deriv3 deal with an indexed term in the formula

[R] deriv when one term is indexed

R.version.string # XP [1] R version 2.4.0 Patched (2006-10-24 r39722) On 11/18/06, Ken Knoblauch [EMAIL PROTECTED] wrote: Thank you for your rapid response. This is reproducible on my system. Here it is again, with, I hope, sufficient detail to properly document what does not work and what

[R] deriv when one term is indexed

packages: boot MASS lattice 1.2-26 7.2-29 0.14-13 Thanks for any suggestions. best, Ken -- Ken Knoblauch Inserm U371 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91

Re: [R] deriv when one term is indexed

instead, there). Thanks for any suggestions. best, Ken -- Ken Knoblauch Inserm U371 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10

[R] deriv when one term is indexed

) Waiting for profiling to be done... 2.5%97.5% Blev -1.612492e-01 2.988387e-02 beta1 1.269000e-05 1.792914e-05 beta2 3.844042e-05 5.388546e-05 beta3 6.108282e-06 8.762679e-06 gamm 2.481102e+00 2.542966e+00 Ken Knoblauch wrote: Thank you for your rapid response

[R] Fw: nested linear model; with common intercept

(?) that my Site-coefficients are now calculated relative to my reference Site (treatment.contrasts), *but* that my TideCode levels now relate to their reference level within Site. Is that correct? Thank you in advance for help. Regards, Mark Difford. -- Ken Knoblauch Inserm U371 Institut

[R] Fitting a cumulative gaussian

If the data asymptote at 0 and 1, then you can use glm with the binomial family with either the logistic or probit links. If the data are from an n-alternative forced choice procedure or if the data do not asymptote at 0 and 1 for some reason or other, then you need to try other procedures.

[R] (no subject)

Yes, you can, as documented in ?connections, but on the Mac you would have to use: read.delim(pipe(pbpaste)) You can also use pbcopy to copy to the clipboard. Dear R users, I am trying to get data from the clipboard into R on MacOSX. I tried the following, but got an error message:

[R] Importing data from clipboard on Mac OSX

Oops, forgot to include the Subject. sorry for that sloppiness. Yes, you can, as documented in ?connections, but on the Mac you would have to use: read.delim(pipe(pbpaste)) You can also use pbcopy to copy to the clipboard. Dear R users, I am trying to get data from the clipboard into R on

[R] Aligning ragged text columns

Hi Hadley, I find that things line up better in data.frames data.frame(c1 = c(Heading 1, , Heading 2, ), + c2 = c(This is some info, about heading 1, This is some info, about heading )) c1c2 1 Heading 1 This is some info 2 about heading 1 3 Heading 2 This

[R] Aligning ragged text columns

On Jul 19, 2006, at 11:40 AM, ken knoblauch wrote: Hi Hadley, I find that things line up better in data.frames data.frame(c1 = c(Heading 1, , Heading 2, ), + c2 = c(This is some info, about heading 1, This is some info, about heading )) c1c2 1 Heading 1 This is some

[R] R on MAC OS X

Looks like you forgot the leading / as in Data-read.table(/Users/SaraMM/PhD/Analises-LitterBags/Dados- Litter.txt,head=T) Dear all, I have been usig R for some time, but now I have a MAC instead of a PC, am I am having problems in reading files... I have tried:

[R] adding grid lines to an xyplot when there is only 1 panel - corrected

, groups = AREA, type = b, auto.key = list(space = right)) Ken Knoblauch -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive Neuroscience 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10

[R] adding grid lines to an xyplot when there is only 1 panel

-2.070V4 V3 -0.042V4 V4 xyplot(est ~ ISITE, data = all2, groups = AREA, type = b, auto.key = list(space = right)) Ken Knoblauch -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive Neuroscience 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0

[R] simple if statement

How about the following, if you really want characters or just leave as factor i - round(runif(10, 1, 4)) years - as.character(factor(i, labels = c(Freshman, Sophomore, Junior, Senior))) HTH, ken __ R-help@stat.math.ethz.ch mailing list

[R] trace of matrix product

of crossprod() or trcrossprod() -- Robin Hankin Uncertainty Analyst National Oceanography Centre, Southampton European Way, Southampton SO14 3ZH, UK tel 023-8059-7743 * -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive Neuroscience 18 avenue du

[R] Plotting a segmented function

') But this seems very cumbersome. Any help is much appreciated. Thanks Jacob -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive Neuroscience 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http

[R] R] highlight an area below a line

the graph in a certain interval, say x=1,2? I understand that I should use polygon, but I have not found the right way to give the result that I want. Thanks a lot! Georg -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive Neuroscience 18 avenue du Doyen Lépine 69500 Bron France tel

[R] R] Sweave scientific real display format (e.g. 5e-12)

figure and couldn't google out how to make construct a pair of \Sexpr s or a LaTeX macro that would include 5\cdot 10^{-12} into the LaTeX output istead of 5e-12 . Any ideas? Thank you Gábor -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive Neuroscience 18 avenue du

[R] Sweave scientific real display format (e.g. 5e-12)

of it. * Dear All, I couldn't figure and couldn't google out how to make construct a pair of \Sexpr s or a LaTeX macro that would include 5\cdot 10^{-12} into the LaTeX output istead of 5e-12 . Any ideas? Thank you Gábor -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive

[R] Tobit estimation?

of freedom, p= 0 n= 100 -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive Neuroscience 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http://www.lyon.inserm.fr/371

[R] transform

I have a=c(1,1,4,3,5); then a [1] 1 1 4 3 5 is there any function to transform a to b: b a [1] 1 3 4 5 How about: unique(sort(a)) or sort(unique(a)) Thanks! -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive Neuroscience 18 avenue du Doyen Lépine 69500 Bron

[R] FW: Re: Doubt about nested aov output

that do this. -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive Neuroscience 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http://www.lyon.inserm.fr/371

Re: [R] FW: Re: Doubt about nested aov output

it there are three different rats labelled 1. There is a rat 1 on treatment 1 and a rat 1 on treatment 2 and a rat 1 on treatment 3. Thus the levels of Rat do not designate the experimental unit, it is the levels of Treatment:Rat that do this. -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept

[R] Generating correlated data from uniform distribution

~~ - (p.dalgaard at biostat.ku.dk) FAX: (+45) 35327907 Ken Knoblauch Inserm U371, Cerveau et Vision Department of Cognitive Neurosciences 18 avenue du Doyen Lepine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: 06 84 10 64

[R] x*x*x*... vs x^n

there! anyone got any comments? -- Robin Hankin Uncertainty Analyst National Oceanography Centre, Southampton European Way, Southampton SO14 3ZH, UK tel 023-8059-7743 Ken Knoblauch Inserm U371, Cerveau et Vision Department of Cognitive Neurosciences 18 avenue du Doyen Lepine 69500

[R] Is it possible to get the first letter of a word?

or What about; strsplit(Training, split=)[[1]][1] [1] T Ken Knoblauch Inserm U371, Cerveau et Vision Department of Cognitive Neurosciences 18 avenue du Doyen Lepine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: 06 84 10 64 10 http

[R] Fitting Theoretical Distributions to Daily Rainfall Data

. fitdistr from the MASS package works quite well, too. Dear List Members, I need a bit help about fitting some theoretical distributions (such as geometric, exponential, lognormal or weibull distribution) Ken Knoblauch Inserm U371, Cerveau et Vision Department of Cognitive

[R] help with kolmogorov smirnov test

What does 'with ties in' mean? with some identical elements (par ex., au moins une paire ex-equo) HTH Ken Knoblauch Inserm U371, Cerveau et Vision Department of Cognitive Neurosciences 18 avenue du Doyen Lepine 69675 Bron cedex France tel: +33 (0)4 72 91 34 77 fax: +33 (0

[R] French Curve

to have used one, could think of a name, but we looked it up in a universal French dictionary on the web, and it came up with ``un pistolet''. Ken Knoblauch Inserm U371, Cerveau et Vision Department of Cognitive Neurosciences 18 avenue du Doyen Lepine 69675 Bron cedex France tel

RE: [R] French Curve

PROTECTED]: From: Ken Knoblauch I remember that my father had a French curve: it was a plastic template used for drawing which had several smooth edges of varying curvature. You could use it to draw a wide variety of curved shapes. No doubt the French called it something else

Re: [R] problem using uniroot with integrate

example, but I'm not sure what I could have added, if I understand what you mean by working example, because my function wasn't working. best, ken Quoting Sundar Dorai-Raj [EMAIL PROTECTED]: Ken Knoblauch wrote on 3/9/2005 10:27 AM: Hi, I'm trying to calculate the value

[R] problem using uniroot with integrate

status major2 minor0.1 year 2004 month11 day 15 language R ken Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen

[R] apply a function to a rolling subset of a vector

this on the list a year or so ago, but I couldn't find anything in the archives. Thanks in advance, Whit [[alternative HTML version deleted]] Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France tel: +33 (0)4 72 91 34 77 fax

Re: [R] Easy cut paste from Excel to R?

minor0.1 year 2004 month11 day 15 language R Mac OS 10.3.8 Excel 2004, version 11.1 (040909) -Don At 11:09 PM +0100 2/17/05, Ken Knoblauch wrote: Here

Re: [R] Easy cut paste from Excel to R?

Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: 06 84 10 64 10 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

[R] Easy cut paste from Excel to R?

Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: 06 84 10 64 10 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] Easy cut paste from Excel to R?

. Quoting Prof Brian Ripley [EMAIL PROTECTED]: On Thu, 17 Feb 2005, Uwe Ligges wrote: Ken Knoblauch wrote: I tried the interesting suggestion below, discussed in several postings yesterday on the help-list, on my Mac (0S 10.3.7) but could not get it to work, as shown in the tests indicated

[R] Visual stimulus presentation using R?

, suggestion? Cheers, Christoph. Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: 06 84 10 64 10 __ [EMAIL PROTECTED] mailing

[R] AIC in glm.nb and glm(...family=negative.binomial(.))

in the two instances. Thank you, in advance. Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: 06 84 10 64 10 __ [EMAIL PROTECTED

[R] adding a line to a single panel of a lattice plot

}, layout=c(3,1),aspect=x, scales=list(cex=1.2), xlab=list(label=x), ylab=list(label=y) Thanks for any help or directions thereto, in advance. Ken Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France tel: +33 (0

[R] Re :adding a line to a single panel of a lattice plot

precociously. Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: 06 84 10 64 10 Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen Lepine

RE: [R] Plot symbols for more than 25 groups

-project.org/posting-guide.html Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 92 34 61 portable: 06 84 10 64 10 __ [EMAIL PROTECTED] mailing

Re: [R] Difficulties in interaction between R and latex (prosper)

] Department of Economic Affairs http://www.mayin.org/ajayshah Ministry of Finance, New Delhi Ken Knoblauch Inserm U 371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 92 34 61 portable: 06 84 10 64

RE: [R] Copula calculation in R?

}} __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- Ken Knoblauch Inserm U371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France Tel: +33 (0)4 72 91 34 77 Fax: +33 (0)4 72 91 34 61 Portable: 06 84 10 64 10 email: [EMAIL

[R] Multiple comparisons with a glm

I've never seen anything written about multiple comparisons, as in the multcomp package or with TukeyHSD, but using a glm. Do such procedures exist? Are they sensible? Are there any packages in R that implement such comparisons? Thank you. -- Ken Knoblauch Inserm U371 Cerveau et Vision 18