search for: ghinea

Dear R Users, There exists a non-central hypergeometric distribution function in the (MCMCpack) package, and a hypergeometric distribution function in the (stats) package. Is there a function for sampling from an extended hypergeometric distribution? Thanks, Narcyz This message is intended for the addressee named and may con...{{dropped}}

Hello fellow R users, Below are two cases using the "mle" operation from the stats4 package. In CASE 1 the code runs fine, in CASE 2 errors occur: CASE 1 x, alpha_current, s, and n are vectors of the same length. ll_beta<-function(b0=0,b1=0) -sum(s*b0+s*b1*x+s*alpha_current-n*log(1+exp(b0+b1*x+alpha_current))) fit_beta<-mle(ll_beta) CASE 2 The error message is as follows

Dear R Users, Any insights into why the following occurs would be helpful.... Firstly: #Evaluating proportions p<-as.list(rep(0,times=length(n))) #creating object of appropriate size. for(j in 1:length(n)){ for(k in 1:length(n[[j]])){ p[[j]][[k]]<- (s[[j]][[k]]/n[[j]][k])}} # 31 x k x num_samples(dim) where k varies The list object s[[j]][[k]] has 3 levels, j, k, and i. j represents the