similar to: gibbs sampling for mixture of normals

Hello All, I would like to optimize a (log-)likelihood function subject to a number of linear constraints between parameters. These constraints are equality constraints of the form A%*%theta=c, ie (1,1) %*% 0.8,0.2)^t = 1 meaning that these parameters should sum to one. Moreover, there are bounds on the individual parameters, in most cases that I am considering parameters are bound between zero

Dear all Currently, I'm implementing the generalized hyperbolic distribution into Splus. Unfortunately the Bessel function is not implemented in Splus. In R the Bessel function does exist but it is an internal function and I'm not able to look at the code. Is there any possibility to see the code of the Bessel function in R or does anybody has an implementation of the Bessel function in

I have a monthly price index series x, the related return series y = diff(log(x)) and a POSIXlt date-time variable dp. I would like to apply annual blocks to compute for example annual block maxima and mean of y. When studying the POSIX classes, in the first stage of the learning curve, I computed the maximum drawdown of x: > mdd <- maxdrawdown(x) > max.dd <- mdd$maxdrawdown > from

Hello, I also try to fit a skewed distribution (like skewed student t) to data points. Do you have an idee howto do this??? thank you fabrice dusonchet *********************************************************************************** This email and any files transmitted with it are confidentia...{{dropped}}

Sigfried, I am not a statistician, but I have learned that according to the Central Limit Theorem (CLT) sums of random variables, regardless of their form, will tend to be normally distributed. CLT does not require the variables in the sum to come from the same underlying distribution. Ciao, Hannu Kahra Progetti Speciali Monte Paschi Asset Management SGR S.p.A. Via San Vittore, 37 IT-20123

my rails version is 2.2.0 i have created a rss feed It is working fine in IE but in fire fox only xml content is showned like "This XML file does not appear to have any style information associated with it. The document tree is shown below." any idea ? -- Posted via http://www.ruby-forum.com/. --~--~---------~--~----~------------~-------~--~----~ You received this message because

Dear all, in the coming Winter Semester, I will be a teaching assistant for a course in Survival Analysis. My job will be to do the lab sessions. The software used for these lab sessions will be R. Most of the students have a background in social sciences and the only stats package they used so far is most likely SPSS. So I assume they might be quite surprised the first time they see R

Hello: would you please help me with the following glm question? for the R function glm, what I understand is: once you specify the "family", then the link function is fixed. My question is: is it possible I use, for example, "log" link function, but the estimation approach for the guassian family? Thanks, Shuangge Ma, Ph.D. ******************************************** *

I have R-0.62.3 and RedHat 5.0 with fort77. At the office R is running but when I install a package, eg. 'akima' I get a message make: f77: Command not found 'R INSTALL library' does not accept fort77, but R installation checks f77/g77/fort77. Is there an easy way to edit some file and replace f77/g77 with fort77? At home I have the same R and Linux versions, but when

Hi, I am trying to maximize a utility function using optim. I have a problem with optim, since some of the values in the caw, mom, me and btm matrixes in the code bellow are missing. Is there a handy way just to skip the missing values in the loop? g <- 5 retp <- NULL object <- function (x) { b1 <- x[1] b2 <- x[2] b3 <- x[3] for(i in 1:nrow(ret)){ for(j

Hi, when using a laptop and trying to install packages for R for Windows 2.2.1 I get the following error > utils:::menuInstallPkgs() --- Please select a CRAN mirror for use in this session --- Warning: unable to access index for repository http://cran.dk.r-project.org/bin/windows/contrib/2.2 Warning: unable to access index for repository

Hi, I have a graph with seven parallel horizontal lines. Is it possible to shade the area between two adjacent lines? Thank you in advance, Hannu

Hi, I have adviced my students to install R and an editor on a USB drive for working in the computer class. With R everything works fine following these instructions: http://personal.bgsu.edu/~mrizzo/Rmisc/usbR.htm. But several editors (e.g., Tinn-R and WinEdt) require administrator rights. I have found out that Emacs (Vincent Goulet's Emacs for Windows Modified

Full_Name: Hannu Kahra Version: 1.0.0 OS: Windows Submission from: (NULL) (153.1.50.127) x1 <- lag(x, k=1) returns the original series x. The time series is not shiftet at all. Only the time is shiftet: for x1 new start = old start -1 and new end = old end -1. When I run, e.g. lm(x~x1) the result is not what I expect. Both x and x1 are time series. Is this a bug or am I forgetting something?

Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file 'a:/project.txt' please help --------------------------------- [[alternative HTML version deleted]]

Hi, I just found out that Tinn-R_2.3.6.0 is now available on http://sourceforge.net/projects/tinn-r/. When configuring Tinn-R: R > Configure > Permanent (Rprofile.site) I get error C:\Program Files\R\R-2.12.0\bin\etc\RProfile.site The above file was not found. Please try to repeat the procedure! The correct path for RProfile.site is C:\Program Files\R\R-2.12.0\etc\Rprofile.site I was

Hi, Can you tell me if there are an R package or code for STAR model estimation and test misspecification. If no, how i could do this. Thanks in advance Best regards AJMI Noomen Phd student TUNISIA --------------------------------- [[alternative HTML version deleted]]

Hi R Users: This package provides two generic functions for performing Markov chain sampling in a naive way for a user-defined target distribution, which involves only continuous variables. The function "gibbs_met" performs Gibbs sampling with each 1-dimensional distribution sampled with Metropolis update using Gaussian proposal distribution centered at the previous state. The function

I would be very grateful for any help from members of this list for what might be a simple problem... We are trying to simulate the behaviour of a clinical measurement in a series of computer experiments. This is simple enough to do in R if we assume the measurements to be Gaussian, but their empirical distribution has a much higher peak at the mean and the distribution has much longer tails.

Hey everyone, I have read the paper by Ritter and Tanner(1992) on Griddy-Gibbs sampler and I am trying to implement it in R without much luck. I was wondering if anyone had used this or could point me to any example code. Thanks, Liz --------------------------------- [[alternative HTML version deleted]]