search for: logprobs

Hello, I have a small problem with developing design matrix X, which I use in estimation the log-likelihood of a multinomial logit model. I have the data: number of observation - 289 number of choice alternative- 3 number of choice specific variables in matrix X -4 matrix X =289x4 I tried to use the function createX, I know that I have to get design matrix 289x12 (am I right?) but

Hello If i want to resampl from the tails of normal distribution , are these commans equivelant??   upper tail:qnorm(runif(n,pnorm(b),1))  if b is an upper tail boundary   or   upper tail:qnorm((1-p)+p(runif(n))  if p is the probability of each interval (the observatins are divided to intervals)   Regards [[alternative HTML version deleted]]

Fellow R users, I'm using the BCE {BCE} function to run a Bayesian sediment mixing model. The aim is to find the optimum % contribution from each of the 4 source areas that can yield the target geochemistry. I have geochemistry for 4 source areas called Rat: Rat<-read.table(text="CaO MgO Na2O Al2O3 Topsoils 2.511250 0.7445500 0.7085500 14.10375 ChannelBanks