小编Spa*_*man的帖子

我有一个庞大的IP地址列表.有没有办法使用R来获得这些地址的相应纬度和经度？

谢谢.

拉维

我有2个节点类型,其中TypeA将始终指向TypeB,TypeB没有出站边.

如何使用igraph将其指示为有向图？

我试图使一个textAreaInput框闪亮,跨越100%的网页,并在浏览器最小化/最大化时调整大小.我可以通过提供参数width = 100%来使用此行为创建一个简单的textInput.即使width在textInput和textAreaInput手册页中具有相同的描述,向textAreaInput提供相同的参数也不会产生相同的行为.这是理想的行为还是错误？

最小的工作示例 -

library(shiny)

shinyApp(
    #UI
    ui = fluidPage(
        fluidRow(
            column(12,
                textAreaInput("big_box", "Big box", value = "", width = '100%', rows = 5, resize = "both")
            )
        ), 
        fluidRow(
            column(12,
                textInput("long_box", "Long box", value = "", width = '100%')
            )
        )
    ),
    #Server
    server = function(input, output) {
    }
)

输出示例 -

干杯

我有一个175行6列的.csv文件。我想追加第176行。我的代码如下：

x <- data.frame('1', 'ab', 'username', '<some.sentence>', '2017-05-04T00:51:35Z', '24')

write.table(x, file = "Tweets.csv", append = T)

我希望看到的是： 在此处输入图片描述

相反，我的结果是： 在此处输入图片描述

我应该如何更改我的代码？

我已经下载了一些气候再分析数据,花了很多时间从grib转换为cdf,最后我已经设法让它在R中工作.我的下一个问题是时间是一个字符向量,类似于下列:

tt =c(
    "20090101","20090101.25","20090101.5","20090101.75","20090102",
    "20090102.25","20090102.5","20090102.75","20090103","20090103.25",
    "20090103.5","20090103.75","20090104","20090104.25","20090104.5", 
    "20090104.75")

前四位数是年份,后面是月份,然后是白天,最后是一天中的一小部分:yyyymmdd ...

如何用格式yyyy-mm-dd HH:MM将其转换为posixct时间？

我会尝试

as.POSIXct(TT, '%Y%米%d')

但由于日间分数,这在这里不会有所帮助.

最终结果应该是:

    tt
       [1] "2009-01-01 00:00"    "2009-01-01 06:00" "2009-01-01 12:00"  
"2009-01-01 18:00" "2009-01-02 00:00"   

谁有人建议解决方案？

注意:我可以用一半的方式

as.Date(tt, "%Y%m%d")

但不知道如何处理分数

我试图通过使用 doParallel 包和 foreach 包来让我的代码在 R 上并行运行。我还使用 sf 包来操作 shp 文件。我确保所有代码都在 foreach 循环中使用 %do% 工作，因此如果出现错误，我可以更好地追踪它。我的代码使用 foreach 和 %do% 工作得很好，但是当我更改它时 do%dopar% R 会继续给我以下错误：

{ 中的错误：任务 1 失败 - “找不到函数“st_geometry_type””

尽管我在 R 脚本的顶部明确使用了 require(sf) 。我制作了一个小函数，如果语句为真，则仅打印“检查”以复制错误。

require(sf)
require(doParallel)
doParallel::registerDoParallel(cores = 2)

testforeach <- function(sfObject)
{
  foreach(i=1:10) %dopar% {
    if (st_geometry_type(sfObject[i,]) == "LINESTRING")
    {
      print("check")
    }
  }
}

当我运行这段代码时，它会抛出同样的错误：

{ 中的错误：任务 1 失败 - “找不到函数“st_geometry_type””

但是，当我用 %do% 替换 %dopar% 时，它会打印出所有预期的“检查”消息。这是 R 中的错​​误还是我遗漏了什么？我尝试重新安装我的软件包，但这似乎没有任何影响，因为我继续遇到相同的错误。任何帮助将不胜感激。

当n增加时,R中双循环的速度非常慢.有没有办法提高for循环的速度？

    set.seed(1)
    n=1000

    y=rnorm(n)
    x1=rnorm(n)
    x2=rnorm(n)

    lm.ft=function(y,x1,x2)
      lm.fit(cbind(1,x1.bar,x2.bar), y)$coef

    res=array(,dim=c(1,3,n,n))
    for(i in 1:n)
      for(j in 1:n){
       x1.bar=x1-x1[i]
       x2.bar=x2-x2[j]
       res[,,i,j]=lm.ft(y,x1.bar,x2.bar)
      }

当在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创建理论分位数的函数，但是如何呢？