小编emu*_*rak的帖子

我刚刚发现了这个fitdistrplus软件包,我把它和Poisson发行版一起运行等等.但是在尝试使用二项式时我遇到了问题:

set.seed(20)
#Binomial distributed, mean score of 2
scorebinom <- rbinom(n=40,size=8,prob=.25)


fitBinom=fitdist(data=scorebinom, dist="binom", start=list(size=8, prob=mean(scorebinom)/8))

我收到错误:

Error in fitdist(data = scorebinom, dist = "binom", start = list(size = 8,  : 
  the function mle failed to estimate the parameters, 
                with the error code 100
In addition: There were 50 or more warnings (use warnings() to see the first 50)

我从这个包中看到了很多关于负二项分布的文档,但关于二项式并不多.这个功能确实出现了支持这种分布(虽然fitdistr在MASS没有).

有什么想法吗？

我正在寻找通过属性对SpatialPolygonsDataFrame进行子集设置，但是我想允许它返回空的SpatialPolygonsDataFrame。

如果像这里讨论的那样将SpatialPolygonsDataFrame类型的对象像data.frames一样对待，我们应该能够获取并使用空对象。

我很感兴趣，因为我想将其合并到一个函数中，该函数可能会尝试通过本质上不选择任何功能的属性进行子集化。

owd <- getwd() 
setwd(system.file("shapes", package = "maptools")) 
library(maptools) 
nc90 <- readShapeSpatial("co37_d90") 
setwd(owd)

nc90@data[nc90@data$AREA>0.15,]               # returns data.frame
bigctys <- nc90[nc90@data$AREA>0.15,]         # SpatialPolygonsDataFrame 
nc90@data[nc90@data$AREA>0.25,]               # returns empty data.frame
bigestctys <- nc90[nc90@data$AREA>0.25,]      # ERROR

有没有办法使这项工作？如果没有，有没有办法初始化一个空的SpatialPolygonsDataFrame对象？我想在这样的对象上执行的将来动作涉及在现有地图上进行绘图，因此我希望无论如何都可以生成图像，即使是空白图像也是如此。

我想更改默认的xyplot设置，其中轴标签按面板交替显示。

xyplot(yield~N | P+K, data=npk)

我意识到这是有意避免在相邻面板上重叠极端轴的方法，但是对于如上所示的分类x轴标签，这是不必要的，对于发布来说看起来很奇怪。

有没有一种方法可以将所有轴标签放在如下所示的同一侧（我在图形程序中编辑过）？