similar to: problem with mode of marginal distriubtion of rdirichlet{gtools}

I am a pretty new R user, I am running the latest linux version on xandros, updated with some extra debian packages, and I also run the latest windows version, but prefer linux. I am having trouble downloading "bayesm", it won't do it all from any of the sites on the web, I resorted to this one, http://packages.debian.org/unstable/math/r-cran-bayesm. and got slightly further, but I

Error: could not find function "hier.part" things I have tried: 1. reinstall R (lastest version, on windows XP) 2. install package "gtools" 3. include: library(gtools) require(gtools) 4. how I call this function: hier.part(c$Y, xdata, fam = "gaussian", gof = "Rsqu") 5.when I try to check what's in the package "gtools", I get (hier.part is

Greetings R Community, I am working with the ifelse function and it is returning something unexpected. In the code the line with the MODE1 assignment the output is a vector [1] 4 5 6 but when I put the MODE1 object into the ifelse function [R}'s output for MODE1 is the first number from the string (4). Why is this? Given the supplied vector of x I would assume both the MODE1 and ifelse()

Dear R People: thank you for the many helpful sets of code that I received!!! I combined several of the concepts for the following function: > mode1 function(x) { y <- rle(sort(x)) z <- y$values[y$lengths==max(y$lengths)] return(z) } > xm [1] 22 15 10 30 25 26 2 17 28 2 24 6 26 24 5 22 20 14 > mode1(xm) [1] 2 22 24 26 > This will pick

Is this a common warning? The decoder doesn't return an error on it, but I see it a lot in my test application on windows. It is non existent on my linux box. I haven't tried mingw yet. please note that I'm using visual studio 2008 w/the vcproj that Bjoern Rasmussen made for 0.5.2 (w/some file references removed) at the moment and it is giving a lot of C4554 warnings

I am attempting to use optim to solve a neural network problem. I would like to optimize coefficients that are currently stored in a matrix Y=270 x 1 X= 27- x 14 b1= 10x14 b2= 11x1 V= 10 x 14 set of prior variances. I have the following function: posterior.mode1=function(y,X,b_0,b2,V) { log.like=function(b1) { a_g=compute(b1) z_g=tanh(a_g); z_g=cbind(1,z_g)

Hi, Just thought I'd share some test results of mine in case anyone is interested. Basically the tests are siumulating what our product does with files - although the tests do it a lot quicker (not as many files though). The test is to create 1 million files (each containing the text of the file number) spread over a number of directories. The files are then removed in the same manner as

Dear Greg, Thanks for the new release. The decomposition of the SSQ is just what I need! Regards, Martin. Martin Hoyle, School of Life and Environmental Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK Webpage: http://myprofile.cos.com/martinhoyle >>> gregory_r_warnes at groton.pfizer.com 10/30/02 07:16PM >>> Version 0.7.3 of the gregmisc package

Hola comunidad, Les consulto mi duda. Quisiera generar N números aleatorios tal que su suma esté predeterminada a un valor k. Sé que con la distribución Dirichlet, por ejemplo, se podrían obtener números aleatorios tal que su suma es la unidad. Por ejemplo MCMCpack::rdirichlet(1,rep(10,7)) genera 7 números tal que su suma es 1 (el parámetro rep(10,7) es cualquier cosa). Pero me pregunto si se

>The way the El Torito CD-ROM is going to look like to the program, is >that it looks like an additional hard disk which a very high drive >number (typically 9Fh or so); MS-DOS will normally ignore this since it >asks the BIOS how many drives there are and ignores the rest. >Does the error message come from UDMA2.SYS or...? > -hpa yes, it comes from UDMA2.SYS, it does a

>From spencer.graves at PDF.COM Mon Jun 23 14:36:45 2003 Received: from postal.pdf.com (pdf193.pdf.com [209.128.81.193]) by uhddx01.dt.uh.edu (8.9.2/8.9.2) with ESMTP id OAA09671 for <hodgess at uhddx01.dt.uh.edu>; Mon, 23 Jun 2003 14:36:44 -0500 (CDT) Received: from postage.pdf.com (postage.pdf.com [10.10.8.7]) by postal.pdf.com (Switch-3.0.4/Switch-3.0.0) with ESMTP id

I would like to code the following in R: a1(b1+b2+b3) + a2(b1+b3+b4) + a3(b1+b2+b4) + a4(b1+b2+b3) or in summation notation: sum_{i=1, j\neq i}^{4} a_i * b_i I realise this is the same as: sum_{i=1, j=1}^{4} a_i * b_i - sum_{i=j} a_i * b_i would appreciate some help. Thank you. -- View this message in context: http://r.789695.n4.nabble.com/summation-coding-tp4646678.html Sent from the R

Hello R-users, I am aware that this is a very basic question, but still... Since mode() function is reserved for something else, is there a function to calculate statistical mode for a vector or a matrxi? i.e. mode.function(c(1,2,2,2,3,3))=2 Thanks in advance, Elena Moltchanova Elena MOLTCHANOVA IIASA International Institute for Applied Systems Analysis A-2361 Laxenburg, Austria E-Mail:

Hi, I'd like to take difference for a sequence a between a_i and a_i-2, for instance, a<-c(2,3,4,8,1) I need (2, 5, -3) as a result. If not using a for loop, can anyone help me? Thanks a lot. Dot -- View this message in context: http://www.nabble.com/R-help--how-to-take-difference-in-next-two-elements-tp18927968p18927968.html Sent from the R help mailing list archive at Nabble.com.

Dear R-Users, I need to find a smooth function f() and coefficients a_i that give the best fit to y ~ a_0 + a_1*f(x_1) + a_2*f(x_2) Note that it is the same non-linear transformation f() that is applied to both x_1 and x_2. So my first question is how can I do it in R? A more general question is this: suppose I have a utility function U(a_i, f()), where f() is say a spline. Is there a general

I'm relatively new to R, but I'm attempting to do a non-linear maximum likelihood estimation (mle) in R, with the added problem that I have a non-linear constraint. The basic problem is linear in the parameters (a_i) and has only one non-linear component, b, with the problem being linear when b = 0 and non-linear otherwise. Furthermore, f(a_i) <= b <= g(a_i) for some (simple) f

I am new to generalized linear models and studying McCullagh & Nelder (1989). Especially, I have a problem resembling the \"cheese taste\" example (5.3.1. p. 109) of the book. I tried to analyse the cheese example with R but failed to do so because R allowed me to use logit link function only with binary family that supposes 0 <= y <= 1. Do I need to scale the y\'s or

It obviously would be nice to have such a mode available, for e.g. DVD audio compression. Apparently, the list doesn''t tell me too much about it. My questions are: 1. What is the current status of the 5.1 channel coupling in Vorbis? 2. If I''ll be interested in participation in its development, what is the recommended reading? -------------- next part -------------- An HTML

Hello, I'm very new to using R, but I was told it could do what I want. I'm not sure how best to enter the information but here goes... I'm trying to transfer the following integral into R to solve for ln(gamma_1), on the left, for multiple instances of gamma_i and variable N_i. gamma_i is, for example, (0, 0.03012048, 0.05000000, 0.19200000, 0.44000000, 0.62566845) N_i (N_1 or

Dear useRs (called Frank Harrell, most likely), after having preached for years to my medical colleagues to be cautious with stepwise selection procedures, they chanted back asking for an alternative when using mixed models. There is a half dozen laXXX packages around for all types of linear models, but as far I see there is none for mixed models such as lme. Even boot.stepAIC (which I