runif function and extreme values:

[Aldi Kraja <aldi@wubios.wustl.edu>]

Thanks to all who responded: Rossini, Connor, Alzola, Plate, Lucke, McKinney
The problem was with the range I was creating by using the round function.
0.5-6.5

Also an interesting idea was given by Lucke:

You are not simulating a six-sided die.  You are sampling a continuous distribution.  Use sample(1:6,size=n,replace=T). Your extremes are diminished because of distribution truncation at 1 and 6. 

McKinney provided a detailed solution:
"rollingADie"<-
function(n, max1)
{
	# n nr. of obs. and max represent max number selected of a rolling die
	y <- sample(1:max1, size = n, replace = T)
	hist(y, main = "The rolling die", probability = T)
	return(NULL)
}

set.seed(123)
rollingADie(100000,6)
Thanks again,
Aldi
--