similar to: Fisher LDA and prior=c(...) argument

> 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 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, I am using the following R version: > R version 2.4.1 (2006-12-18) > Copyright (C) 2006 The R Foundation for Statistical Computing > ISBN 3-900051-07-0 I believe I found a bug in the 'hist' function, when 'probability=TRUE'. I looked in the archives and I came across problems with the 'hist' functions (e.g., bug PR# 944, posted in 2001), however, a

installed from the binary image (http://cran.r-project.org/bin/ macosx/) on an intel mac 10.4.10 throws three warnings: ----------------- R version 2.5.1 (2007-06-27) Copyright (C) 2007 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or

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

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

Hi all How can I change the limits (xlim or ylim) in a plot that has been already created? For example, consider this naive example curve(dbeta(x,2,4)) curve(dbeta(x,8,13),add=T,col=2) When adding the second curve, it goes off the original limits computed by R for the first graph, which are roughly, c(0,2.1) I know two obvious solutions for this, which are: 1) passing a sufficiently large

Version 1.1 of the MiscPsycho package had been uploaded to CRAN. The package has been updated to include the following: 1) The irt.ability() function that estimates examinee ability given a set of item parameters. The function is very general and can be used to estimate ability when there are only dichotomous items (1-, 2-, or 3PL), only polytomous items (generalized partial credit model), or a

Version 1.1 of the MiscPsycho package had been uploaded to CRAN. The package has been updated to include the following: 1) The irt.ability() function that estimates examinee ability given a set of item parameters. The function is very general and can be used to estimate ability when there are only dichotomous items (1-, 2-, or 3PL), only polytomous items (generalized partial credit model), or a

Dear R-help list, I've got a data set in long format - each subject can have several (varying in number) measurements, with each record representing one measurement. I want to assign a sequence number to each measurement, starting at 1 for a person's first measurement. I can do this with the by function, but there must be an easier way. Here's my code - id is id number, age is

Hello, Is it possible to build a function, a posteriori, by getting the names of the arguments and the formula in a list or a vector ? When I use the function as.function(), I don't know how to define explicitly the names of my variables as formal arguments. Thanks _______________________________________ Florent BATY CNRS UMR 5558 Dynamique des Populations Bact?riennes Facult? de

Hola. Me interesa meter una serie de regresiones lineales en una lista a la que voy a denominar L. El problema es que a posteriori necesito sacar valores de esas regresiones. Sin embargo, al llamar a L[1] este no es un objeto del tipo "lm" sino que es del tipo "list" y no puedo aplicarle funciones relacionadas con la regresión (tal y como puede ser hacer un summary para

Jean Marc, I have been fiddling around with preprocess to try to get it sounding less metallic/underwater- so based on some things said recently, I have been first fiddling with mean_post and mean_prior to try to get them to do change the update rate, so that update_noise gets called at the right times. Haven't made it too deep into the code, but I have a few questions in speex_preprocess:

Hi all, If I run LDA on the same data (2-class classification) with default(no priors specified in the lda function) vs. "prior=c(0.5, 0.5)", the results are different. The (0.5, 0.5) priors give better 1-classify-to-1 rate, and the proportional priors(default, no priors specified in the lda function) give better 0-classify-to-0 rate, for both training and testing data sets. However,

Dear all, I am studying a dataset using the aov() function. The independant variable 'cds' is a factor() with 8 levels and here is the result in studying the dependant variable 'rta' with aov() : > summary(aov(rta ~ cds)) Df Sum Sq Mean Sq F value Pr(>F) cds 7 0.34713 0.04959 2.3807 0.02777 Residuals 92 1.91635 0.02083 The dependant variable

Hi All, I have data about seed predation (SP) in fruits of three differents colors (yellow, motted, dark) and in two fruiting seasons (2007, 2008). I performed a GLMM (lmer function, lme4 package) and the outcome showed that the interaction term (color:season) was significant, and some combinations of this interaction have significant Pr(>|z|), but I don't think they are the right

[Mod: Sent to linux-security instead of linux alert -- alex] Dave G. <daveg@ESCAPE.COM> wrote: > /* vixie crontab buffer overflow for RedHat Linux > * > * I dont think too many people know that redhat uses vixie crontab. > * I didn''t find this, just exploited it. The vulnerability involves an unguarded sscanf call in env.c. Enlarging the buffer to the largest

Hi, everyone. GAD package analyses complex ANOVA models with any combination of orthogonal/nested and fixed/random factors, as described by Underwood (1997). There are two restrictions: (i) data must be balanced; (ii) fixed nested factors are not allowed. Homogeneity of variances is checked using Cochran's C test and "a posteriori" comparisons of means are done using

Good afternoon, gentlemen! After several days studying and researching on categorical data (various forums with answers from the owner of the library - all incipient) how to interpret the output the function MCMCglmm, come to enlist the help of you, if someone has already worked with MCMCglmm function in the case of variables ordinal dependent. I've read and reread all the pdf's of the

Hola doblett. No sé si te lo han comentado ya, pero yo le echaría un ojo al paquete ReporteRs. Un saludo, Miguel. El 06/06/2017 a las 16:00, doblett escribió: > Muchas gracias a todos por las sugerencias y por la rápida respuesta. > Mi impresión es que no tenemos una solución cómoda para este tema y al > menos en mi día a día es muy útil. > > - Tables, funciona muy bien