小编Spa*_*man的帖子

我正在尝试使用rCharts包显示R的地图.我开始简单,所以我想在地图上添加一个多边形.但我不知道如何.有任何想法吗？addPolygon不起作用.

map <- Leaflet$new()


map$tileLayer(provider = 'Stamen.TonerLite')

map$setView(c(48.1, 16.7), zoom = 10)
map$addPolygon(
  c(48.99831, 49.08815, 49.08815, 48.99831, 48.99831),
  c(13.42666, 13.42666, 13.56383, 13.56358, 13.42666),
  layerId=c("1"),
  options=opts,
  defaultOptions=opts)
map

当在R中使用hexbin并选择大量的xbins时,我会得到重叠的六边形bin,这是一个例子:

library(hexbin)
ss<-10000
x <- runif(ss,0,1)
y <- runif(ss)
bin <- hexbin(x, y,xbins=100,xbnds=c(0,1),ybnds=c(0,1),IDs=TRUE)
plot(bin)

手动读出箱子表明问题不是由多边形的厚度引起的.这真的是一个错误,或者我只是忽略了什么？如果它是一个错误什么是R中的一个好的替代品(不能在ggplot中使用stat_binhex,因为我需要二进制文件以及数据点包含在哪个bin中供以后导出的完整信息)？

注意:当绘图用pdf(可能不在GUI中)时,重叠变得清晰可见.

我必须为两组大小为N的向量生成随机数.

一组的概率是p,另一组的概率是q = 1-p.(例如,对于p = 0.5的1000个群体,我必须从分布中生成500个随机数,从另一个生成500个随机数).由于这是一个我必须改变'p'的模拟,我编写了我的代码来生成如下:

group1 = rnorm(n = N*p)
group2 = rnorm(n = N*q) # 1st method
group2 = rnorm(n = (N - N*p)) # 2nd method    

使用上述两种方法,R产生的随机数少于它在group2的几行中的应用数量(第一行约占35%,第二种约占12%).

我遇到了与rexp,rpois和runif相同的错误.

以下是两种方法的快照供您参考.

#### EXAMPLE SCRIPT #####

N = 1000
p1 = seq(0.01, 0.99, 0.001)
q1 = 1 - p1


### FIRST METHOD ###

X = data.frame()
for (i in 1:length(p1))
{
X[i, 1] = p1[i]
X[i, 2] = q1[i]
X[i, 3] = length(runif((N * X[i, 1])))
X[i, 4] = length(runif((N …

我需要创建一个 QQ 图来检查我观察到的数据是否符合泊松分布。

这是我的数据框：

df = read.table(text = 'Var1 Freq
 1975   10
 1976   12
 1977    9
 1978   14
 1979   14
 1980   11
 1981    8
 1982    7
 1983   10
 1984    8
 1985   12
 1986    9
 1987   10
 1988    9
 1989   10
 1990    9
 1991   11
 1992   12
 1993    9
 1994   10', header = TRUE)

这df$Freq专栏是我感兴趣的专栏，因为观察结果代表了每年的事件数量。

我知道我必须使用该qqplot函数以及qpois创建理论分位数的函数，但是如何呢？

有没有办法附加两个图形,我用命令获得par(mfrow=c(2,1)),这样两个绘图窗口之间的空间被删除？这可以通过任何图像编辑器轻松完成,但我猜测是否有更公平的方法可以从命令行执行此操作.

提前致谢!

EB

使用版本0.98.11013.假设我在我的脚本中测试以下代码:

for (num in 1:5)
{
    # Print numbers from 1 to the loop variable
    print(1:num)
}

当我跨过每一行(使用CTRL + R)时,我在控制台中获得以下输出:

> for (num in 1:5)
+ {
+ # Print numbers from 1 to the loop variable
+ print(1:num)
+ }
[1] 1
[1] 1 2
[1] 1 2 3
[1] 1 2 3 4
[1] 1 2 3 4 5
> 

现在,说我在我的代码中犯了一个错误并忘记了结束括号}:

for (num in 1:5)
{
    # Print numbers from 1 to the loop variable
    print(1:num) …

我正在使用 RStudio 在 Windows 计算机上运行一些 R 代码，我的代码会定期生成 Excel 文件和 netCDF 文件（最终会生成数十个）。我不希望它们弄乱我的工作目录。有没有办法将文件保存到父目录中名为“输出”的目录（例如：C：/.../原始文件路径/输出）？我想要一种将当前工作目录更改为其他目录的方法。我知道有 getwd() 和 setwd() 但是我如何设置输出目录的路径而不输入整个 Windows 路径（例如：setwd（windows 或 Mac/output 的当前源文件路径）。我的合作者使用一台 Mac，他的输出也存储在那里。

我有这个R代码,我想分配给最大变量A和B的最大值,如下所示:

Max<-(A,B)

我怎么能在R？我搜索但我找不到最大功能.谢谢