Displaying 2 results from an estimated 2 matches for "sdplus".
Did you mean:
dplus
2010 Aug 05
2
try-error within for loop
...(1,s)))
taxdiv<-taxondiveO(abd,taxa)
if(class(result)=="try-error") result[[i]] <- NA
else result[[i]] <-taxdiv #it doesnt work
result[[i]]<-as.data.frame(result[[i]])[1,]},silent=T)
taxind<-do.call(rbind,result)
taxind[1:5,]
number D Dstar Lambda Dplus SDplus
1 589 94.72456 94.93875 142.4468 94.88140 55885.15
2 428 80.38656 80.63660 328.7728 80.90432 34627.05
3 282 88.03613 88.45139 332.8963 88.53554 24967.02
4 1320 80.43822 80.51943 328.1416 80.63355 106436.29
5 1077 80.43121 80.53073 287.2267 80.62954 86838.01
and instead I want...
2001 Oct 02
4
plot of Bernoulli data
I have some Bernoulli data something like this:
x<-sort(runif(100,1,20))
p<-pnorm(x,10,3)
y<-as.numeric(runif(x)<p)
plot(x,y)
lines(x,p)
This plot is not very satisfactory because the ogive does not visually
fit the (0,1) points very well, and also because the points tend to fall
on top of one another. The second problem can be eliminated by adding
vertical jitter. However I was