File: dmnom.R

package info (click to toggle)
r-cran-combinat 0.0-7-1
  • links: PTS, VCS
  • area: main
  • in suites: squeeze
  • size: 124 kB
  • sloc: makefile: 1
file content (15 lines) | stat: -rw-r--r-- 373 bytes parent folder | download | duplicates (5) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
"dmnom"<-
function(x, size = sum(x), prob = stop("no prob arg"))
{
#       DATE WRITTEN: 22 May 1995           LAST REVISED:  22 May 1995
#       AUTHOR:  Scott Chasalow
#
	p <- max(length(x), length(prob))
	x <- rep(x, length = p)
	prob <- rep(prob, length = p)
	prob <- prob/sum(prob)
	if(sum(x) != size)
		0
	else exp(logfact(size) + sum(x * log(prob) - logfact(x)))
}