Search Archives
Search tips
from:Name
Search by author name, e.g. from:Duncan Murdoch
"exact phrase"
Match an exact phrase
word1 word2
Match messages containing both words
Date range
Use the date pickers to filter results to a time period
Use the list dropdown to narrow results to a specific mailing list. Combine from: with other terms to filter by author and content.
Package mvtnorm provides dmvnorm, pmvnorm that can be used to compute Pr(X=x,Y=y) and Pr(X<x,Y<y) for a bivariate normal. Are there functions that would compute Pr(X<x,Y=y)? I'm currently...
Berton Gunter <gunter.berton at gene.com> wrote: > plot(myscores[,1:2],pch=c('s','c','v')[groups]) Thanks, this works. How to understand the result of the expression > c(1,2,'3)[groups] [1] 3 3 3 3 3...
Sasha Pustota <popgen at gmail.com> wrote: > Berton Gunter <gunter.berton at gene.com> wrote: > > plot(myscores[,1:2],pch=c('s','c','v')[groups]) > > Thanks, this works. How to understand the result of > the expression > > > c(1,2,'3...
On p. 48 of "Statistics Complements" to the 3rd MASS edition, http://www.stats.ox.ac.uk/pub/MASS3/VR3stat.pdf I read that the orthogonal rotations of Z Lambda^-1 remain uncorrelated, where Z is the PC and Lambda...
Ok, that was just my wishful thinking. Is there a way to plot repeated labels that identify groups, e.g. factor(c(rep("s",50),rep("c",50),rep("v",50))) instead of 1--150 row indices, using something like...
I hope this time I'm using the "iris" dataset correctly: ir <- rbind(iris3[,,1], iris3[,,2], iris3[,,3]) lir <- data.frame(log(ir)) names(lir) <- c("a","b","c","d") I'm trying to understand the meaning of expressions...
Suppose I compute a correlation matrix ... r <- cor(log(rbind(iris3[,,1], iris3[,,2], iris3[,,3]))) and then want to plot the values in color, e.g. image(r, col = terrain.colors(200)) Why the matrix in the plot is...
jim holtman <jholtman at gmail.com> wrote: > does this explain it? > > > groups <- factor(c(rep("Z",5),rep("X",5),rep("Y",5))) > > > > groups > [1] Z Z Z Z Z X X X X X Y Y Y Y Y...
I'm trying to integrate one dimension of a bivariate normal. It works when additional parameters are passed explicitly: library(mvtnorm) bivInt <- function(x,y,mx,my,r) { dmvnorm(c(x, y), mean=c(mx, my), sigma=rbind(c(1...
It's in the Venables & Ripley MASS (ed 3) book in the section on principal components. The context is as follows > ir <- rbind(iris[,,1], iris[,,2], iris[,,3]) > ir.species <- factor(c(rep("s",50),rep("c",50),rep...
I'm trying to follow an example from ?clara, which plots two clusters: > x <- rbind(cbind(rnorm(10,0,8), rnorm(10,0,8)), cbind(rnorm(20,50,8), rnorm(20,50,8))) > clarax <- clara(x, 2) > clusplot(clarax) Suppose...
Gabor Grothendieck <ggrothendieck at gmail.com> wrote: > Sasha Pustota <popgen at gmail.com> wrote: > > ir <- rbind(iris3[,,1], iris3[,,2], iris3[,,3]) > > lir <- data.frame(log(ir)) > > names(lir) <- c("a","b","c","d") > > > > I'm trying to understand the...
I wrote: > On p. 48 of "Statistics Complements" to the 3rd MASS edition, > http://www.stats.ox.ac.uk/pub/MASS3/VR3stat.pdf > I read that the orthogonal rotations of Z Lambda^-1 remain > uncorrelated, where Z is the PC...
Brian Ripley <ripley at stats.ox.ac.uk> wrote: > 'iris' in S-PLUS is not the same as 'iris' in R, rather similar to 'iris3' > in R. Thank you! This answered all my questions. > You need the fourth (2002) edition...
Thanks Jay. I realized that I was doing it a silly way shortly after I posted and that the answer i was looking for is simply condXY(y, x, my, mx, r) * dnorm(y, my) condXY <- function(y, x, my...
P.S. just "pnorm(x, mean=m, sd=s)", not "1-pnorm(x, mean=m, sd=s) + pnorm(-x, mean=m, sd=s)" On Wed, Oct 1, 2008 at 2:04 PM, Sasha Pustota <popgen at gmail.com> wrote: > Thanks...
Can't find what you're looking for? Try searching with Google .