search for: airoldi

...e CODE above estimates the likelihood of of the Fisher scores for (example | class) and then implements the Bayes rule to return the maximum a-posteriori class. Is that correct? Any pointer towards that direction is appreciated. Please cc to edo at stat.cmu.edu your reply. Thanks Edoardo M. Airoldi http://www.stat.cmu.edu/~eairoldi BH 232L (412) 268.7829 PC Lab (412) 268.8719

...25)) > hist(z,breaks=seq(-1,1,by=.25),probability=TRUE) I think the values on the Y axis are messed up, e.g., it should top at 0.3 (relative frequency) for the bin [0.25 0.50). How is 'Density' computed? best regards, Edo ------------------------------------------------- Edo Airoldi, Ph.D. Department of Computer Science & Lewis-Sigler Institute for Integrative Genomics Princeton University, NJ 08544 609-258-8326 (lab phone) 609-258-8004 (fax)

...ing reply message sequence 1 on thread 0x15b5e250 2007-08-15 20:52:09.414 R[711] tossing reply message sequence 2 on thread 0x15b5e250 2007-08-15 20:52:09.454 R[711] tossing reply message sequence 3 on thread 0x15b5e250 > ------------------ however nothing obvious seems to go wrong. edo airoldi

> hi all, > I am using optim. I am getting the following error message: > > Error in optim(par = start.vals[, h], fn = post.func.pois, gr = post.grad. > pois, : > L-BFGS-B needs finite values of fn > > If I look at optim typing '> optim' it seems that the error comes from > inside .Internal(optim), so I wonder how can I see the code for .Internal(

Hello all, I am using libraries (mva,cluster) to produce dendrograms. With 1000 examples the dendrogram gets too crowded, and i am wondering whether there is an option (which i cannot find) to set the number of leaf nodes, like in matlab, and return the plot and the assignment map examples -> leaf nodes. Any suggestion is appreciated. Thanks Edo

hello all, I am trying to print a ps file as part of a function as in: func <- function (..., filename="temp.ps") { # some stuff [...] # plot eval( cat("postscript(\"",filename,"\")\n", sep="") ) plot(...) abline(...) dev.off() # more stuff [...] } but it does not work. Nor it does with 'paste' instead

hi there, i have one more question about LDA. just to make surei understand, suppose we have two classes, then if i specify a prior=c(.3,.7) in lda(...) this will affect my between classes covariance matrix as in: SB = (.3*m1 - .7*m2) %*% inv(Sigma) %*% t(.3*m1 - .7*m2) [is Sigma affected ?] and the threshold to decide which class to assign 'test' data = log(.3/.7) if i specify a

...umley) 17. Re: Smooth of a time serie ( Henrique Patr?cio Sant'Anna Branco ) 18. degrees of freedom in a LME model (Federico Calboli) 19. RE: within group variance of the coeficients in LME (J.R. Lockwood) 20. lm diagnostics and qr (fwd) (Jean Eid) 21. Re: dendrograms (Edoardo M Airoldi) 22. Re: lm diagnostics and qr (fwd) (John Fox) 23. assignment in lists (Philippe Grosjean) 24. problems with library in 1.7.1 (R. Heberto Ghezzo) 25. Re: Can't save a graph to pdf in R for MacOS (p.b.pynsent) 26. Re: Plots using POSIX (Prof Brian Ripley) 27. Re: problems with libra...

hi all, I am fitting a logistic regression model on binary data. I care about the fitted probabilities, so I am not worried about infinite (or non-existent) MLEs. I use: > glm(Y~., data=X, weights=wgt, family=binomial(link=logit), maxit=250) I understand the three ways to fit model, and in my case Y is a factor, one column > Y <- c(rep("A",679), rep("B",38))

hi all, I am fitting a logistic regression model on binary data. I care about the fitted probabilities, so I am not worried about infinite (or non-existent) MLEs. I use: > glm(Y~., data=X, weights=wgt, family=binomial(link=logit), maxit=250) I understand the three ways to fit model, and in my case Y is a factor, one column > Y <- c(rep("A",679), rep("B",38))