From: "G. Sawitzki" <gs at statlab.uni-heidelberg.de>
Date: 28 agosto 2005 14:11:18 GMT+02:00
To: jago at mclink.it
Dear Stefano,
this small exaple leads to a crash of R. I did not try it on
versions other than the Mac version. So I am sending it to you
directly. If it is a littice problem, could you pass it to Deepayan
Sarkar? Thank you.
g.
==
#pbinom
library(grid)
library(lattice)
n<-20
psteps<-50
binomtable<- function (n,psteps){
x<-(0:(10*n))/10
p<- (0:psteps)/psteps
dd<-expand.grid(x=x,p=p)
dd$F<- pbinom(dd$x,n,dd$p)
dd$x0<-trunc(dd$x)
dd
}
bt<-binomtable(n=5,psteps=100)
bt[bt$x-bt$x0>=0.9,]$F<-NA
wireframe(bt$F~bt$x*bt$p,bt,groups=bt$x0,shade=TRUE) # leads to R
crash
#wireframe(bt$F~bt$x*bt$p,bt,shade=TRUE) #ok