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

Excuse me for forgetting sessionInfo (below) Ken Knoblauch 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 >

[R] plotting predicted curves with log scale in lattice

panel.xyplot(x, y, ...) } ) Note 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 Cerv

[R] Sweave rendering of simple character

; ?HSP 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 Fran

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

h/pipermail/r-sig-mac/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 Kn

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

William Simpson 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 trial. > On e

Re: [R] A More efficient method?

Keith Alan Chamberlain 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 1 -1 -1 1

Re: [R] if statement

e=1 and the bottom > > 30% as Case=0. How do I do that? > > Or, how about, x <- rnorm(100) Case <- 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

Re: [R] Subset and logical operator error

Sébastien 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 > subset

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

see mvrnorm in MASS and especially the empirical argument James Milks 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 way to do tha

[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] [R-pkgs] new packages psyphy and MLDS

cisms, 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

[R] Selecting complementary colours

at 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. _ Prof

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

%*%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

[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

[R] problem with mfg argument of par

grid(1:2, 1:2)) 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

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

a time, so the for-loop 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 Kn

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

; locale: > en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8 > > attached base packages: > [1] "stats" "graphics" "grDevices" "utils" "datasets" > "methods" "base" > > other attached package

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

tils" "datasets" > "methods" "base" > > other attached packages: > MASS KernSmooth > "7.2-30" "2.22-19" > > Otherwise, I'm as stumped as you are. > > On 3/11/07, Ken Knoblauch <[EMAIL PROTECTED]> wr

[R] using scan to record user's input

t; "grDevices" "utils" "datasets" "methods" [7] "base" other attached packages: MASS "7.2-31" -- Ken Knoblauch Inserm U846 Institut Cellule Souche et Cerveau Département Neurosciences Intégratives 18 a

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

ermail/r-help/2006-December/122353.html ken Douglas Bates a écrit : > Was make.link() used in the example code? > > On 2/11/07, Ken Knoblauch <[EMAIL PROTECTED]> wrote: >> >> Isn't it the case, that since R 2.40 that all one ought to need do is >> define one

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

ogexposure 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: +33 (

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

rs > under the plot. Could someone please suggest 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 Knob

[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 2

[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 t

[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] how to adjust link function in logistic regression to predict the proportion of correct responses in 2AFC task?

ripley at stats.ox.ac.uk > Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ > University of Oxford, Tel: +44 1865 272861 (self) > 1 South Parks Road, +44 1865 272866 (PA) > Oxford OX1 3TG, UKFax: +44 1865 2725

[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, ...) panel.abli

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 f

Re: [R] deriv when one term is indexed

un] * GL^gamm, >> +data = dd, start = st >> +) >> > >> > confint(dd.nls) >> Waiting for profiling to be done... >> 2.5%97.5% >> Blev -1.612492e-01 2.988386e-02 >> beta1 6.10

[R] deriv when one term is indexed

beta.rouge", "beta.vert", "beta.bleu", "gamm"))) >> +.grad[, "Blev"] <- 1 >> +.grad[1:17, "beta.rouge"] <- .expr1[1:17] >> +.grad[18:34, "beta.vert"] <- .expr1[1:17] >> +.grad[35:51, &qu

[R] deriv when one term is indexed

.d2) 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 re

Re: [R] deriv when one term is indexed

t; +data = dd, start = st) >> >> confint(dd.nls) > Waiting for profiling to be done... >2.5%97.5% > Blev -1.612492e-01 2.988388e-02 > beta1 6.108282e-06 8.762679e-06 > beta2 1.269000e-05 1.792914e-05 > beta3 3.844042

[R] deriv when one term is indexed

version 2.4.0 Patched (2006-11-10 r39843) i386-apple-darwin8.8.1 locale: C attached base packages: [1] "stats" "graphics" "grDevices" "utils" "datasets" [6] "methods" "base" other attached packages: boot MASS lat

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

de > 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 Cellule Souche et Cerveau Département Neurosciences Intégratives 18 avenue du Doyen Lépine

[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. Tw

[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] (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] Aligning ragged text columns

" ") dd Heading 1 This is some info about heading 1 Heading 2 This is some info about heading On Jul 19, 2006, at 11:40 AM, ken knoblauch wrote: > Hi Hadley, > > I find that things line up better in d

[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

[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

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

V4 >-4.126V4 V1 >-1.481V4 V2 >-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

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

-1.481V4 V2 -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

[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

t also uses %*% instead 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 Neuroscien

[R] Plotting a segmented function

3,f4) plot(x,F,type='l') 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

[R] R] highlight an area below a line

how can I highlight the area below 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 D

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

ot;$5cdot ^{-12}$" untried, of course, but maybe you could squeeze something out 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

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

squeeze something out 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

[R] Tobit estimation?

x2 1.5945012.825190 -2.985390 Scale= 1.493047 Loglik(model)= -66.1 Loglik(intercept only)= -109.7 Chisq= 87.28 on 2 degrees of freedom, p= 0 n= 100 > -- Ken Knoblauch Inserm U371 Cerveau et Vision Dept. of Cognitive Neuroscience 18 avenue du Doyen Lé

[R] transform

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 B

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

. 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

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

not designate the "experimental unit", it is the levels of Treatment:Rat 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

[R] Generating correlated data from uniform distribution

(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~ - (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)

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

about the same. I was not expecting a factor of 3 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 De

[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

[R] Fitting Theoretical Distributions to Daily Rainfall Data

ogram. 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 Dep

[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 f

RE: [R] French Curve

, Andy" <[EMAIL 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 t

[R] French Curve

here, of age 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 6

Re: [R] problem using uniroot with integrate

rking 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, > &

[R] problem using uniroot with integrate

system powerpc, darwin6.8 status major2 minor0.1 year 2004 month11 day 15 language R ken ____ Ken Knoblauch Inserm U 371 Cer

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

ments 1 to 20, the second element of y would contain the sum of elements 2:21, and so on. I thought I had seen 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

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

nhagen Denmark Ph: (+45) 35327918 > ~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 > 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

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

t" "opened" "no" "yes" > > > > version > _ > platform powerpc-apple-darwin6.8.5 > arch powerpc > os darwin6.8.5 > system powerpc, darwin6.8.5 > status

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

ack down readTableHeader yet. The warning occurs even without headers in the data. Quoting Prof Brian Ripley <[EMAIL PROTECTED]>: > On Thu, 17 Feb 2005, Uwe Ligges wrote: > > > Ken Knoblauch wrote: > > > >> I tried the interesting suggestion below, discussed i

[R] Easy cut & paste from Excel to R?

t;clipboard", "r") : unable to open connection In addition: Warning message: cannot open file `clipboard' platform powerpc-apple-darwin6.8 arch powerpc os darwin6.8 system powerpc, darwin6.8 s

[R] Visual stimulus presentation using R?

gt; >So what I need is: > > 1. put images onto the (full!)screen (qick) > 2. read keyboard input > 3. write results (to an R structure, presumably) > >Any idea, suggestion? > > >Cheers, > Christoph. > > > _

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

e in what is calculated 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 __

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

he solution. Sorry to have bothered the help 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 3

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

,lty=2) #would like this only on panel 1 }, 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 C

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

h/mailman/listinfo/r-help > PLEASE do read the posting guide! > http://www.R-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

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

oing wrong? How do I get the graph made using R > to sit horizantally (i.e. landscape), and fill the screen? I tried to > say "horizontal=T" and that doesn't work. > > -- > Ajay Shah Consultant > [EMAIL PROTECTED] Department of Economic Af

RE: [R] Copula calculation in R?

mley Assoc. Professor, Biostatistics > > [EMAIL PROTECTED] University of Washington, Seattle > > > > __ > > [EMAIL PROTECTED] mailing list > > https://www.stat.math.ethz.ch/mailman/listinfo/r-help > >

[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 Visi