javascript

java

python

c#

android

c++

node.js

php

html

jquery

ios

reactjs

css

.net

git

ruby-on-rails

sql

c

ruby

string

排序-R中的等级和顺序

我无法理解R函数rank和R函数order之间的区别。它们似乎产生相同的输出:

> rank(c(10,30,20,50,40))
[1] 1 3 2 5 4
> order(c(10,30,20,50,40))
[1] 1 3 2 5 4

有人可以帮我一下吗?谢谢

trans by 2020-08-11T22:04:39Z

创建本地R软件包存储库

我想创建一个本地R软件包存储库,以便公司中的用户可以从中安装软件包,系统管理员可以定期更新本地存储库。 当前拒绝访问CRAN镜像。

有没有简单的方法可以做到这一点?

trans by 2020-08-11T12:41:26Z

r-加快write.tab的性能

我有data.frame,我想写出来。 我的data.frame的尺寸是256行乘65536列。 什么是write.csv的更快替代品?

trans by 2020-08-11T10:34:26Z

在R命令行中显示时钟

我想知道是否有一种方法可以在R命令行中显示当前时间,例如在MS DOS中,我们可以使用

Prompt $T $P$G

在每个提示行中都包含时钟。就像是

options(prompt=paste(format(Sys.time(), "%H:%M:%S"),"> "))

会做到这一点,但随后会在设置时固定下来。 我不确定如何使其自动更新。

trans by 2020-08-11T09:54:03Z

R中的NULL和character(0)有什么区别?

NULL和character(0)有什么区别? 整数(0)等?

> identical(NULL, character(0))
[1] FALSE

> is.null(integer(0))
[1] FALSE

> str(character(0))
 chr(0) 

> str(NULL)
 NULL

通常,您似乎可以将is.integer0作为参数传递给函数,并且空向量通常返回为character(0)integer(0)等。

为什么会这样呢? 试想一下,是否有零度测试,la is.integer0

trans by 2020-08-11T08:47:41Z

r-具有xyz坐标的3d表面图

我希望有经验的人可以帮助您如何从xyz数据准备形状文件。 尽管未提供创建形状文件的上述步骤,但在这里可以看到一个很好的数据集,很好地说明了彗星Churyumov–Gerasimenko。

我试图更好地了解如何将曲面应用于给定的XYZ坐标集。 通过R包“ rgl”可以直接使用笛卡尔坐标,但是环绕起来的形状似乎更加困难。 我找到了R包geometry,它提供了QHULL函数的接口。 我尝试使用它来计算Delaunay三角剖分面,然后可以在delaunayn(bunny, options="Qt")中进行绘制。我无法找出与函数delaunayn相关的某些选项,以可能控制这些分面的最大计算距离。 我希望这里的某人可能对从xyz数据改善表面构造有一些想法。

使用“斯坦福兔子”数据集的示例:

library(onion)
library(rgl)
library(geometry)
data(bunny)

#XYZ point plot
open3d()
points3d(bunny, col=8, size=0.1)
#rgl.snapshot("3d_bunny_points.png")

#Facets following Delaunay triangulation
tc.bunny <- delaunayn(bunny)
open3d()
tetramesh(tc.bunny, bunny, alpha=0.25, col=8)
#rgl.snapshot("3d_bunny_facets.png")

enter image description here

这个答案使我相信Qhull的R实现可能存在问题。 另外,我现在尝试了各种设置(例如delaunayn(bunny, options="Qt")),但收效甚微。 Qhull选项在此处概述

编辑:

这是一个球体的附加(更简单)示例。 即使在这里,小平面的计算也不会总是找到最近的相邻顶点(如果旋转球,您会看到一些小平面穿过内部)。

library(rgl)
library(geometry)
set.seed(1)
n <- 10
rho <- 1
theta <- seq(0, 2*pi,, n) # azimuthal coordinate running from 0 to 2*pi 
phi <- seq(0, pi,, n) # polar coordinate running from 0 to pi (colatitude)
grd <- expand.grid(theta=theta, phi=phi)

x <- rho * cos(grd$theta) * sin(grd$phi)
y <- rho * sin(grd$theta) * sin(grd$phi)
z <- rho * cos(grd$phi)

set.seed(1)
xyz <- cbind(x,y,z)
tbr = t(surf.tri(xyz, delaunayn(xyz)))
open3d()
rgl.triangles(xyz[tbr,1], xyz[tbr,2], xyz[tbr,3], col = 5, alpha=0.5)
rgl.snapshot("ball.png")

enter image description here

trans by 2020-08-11T02:56:43Z

r-在facet_wrap p中添加“浮动”轴标签

我和这个用户有同样的问题-我有一个“锯齿状”的多面图,其中底行的面板少于其他行,并且我希望每列底部的X轴刻度线。

对于该问题,建议的解决方案是设置grid(在ggplot 0.9.2.1中;我认为我要查找的行为在较早的版本中是默认的。)在我的情况下,这是一个较差的解决方案:我的实际轴标签会很长 ,因此将它们放在每行下面会占用太多空间。 结果是这样的:

 x <- gl(3, 1, 15, labels=paste("this is a very long axis label ", letters[1:5]))
 y <- rnorm(length(x))
 l <- gl(5, 3, 15)
 d <- data.frame(x=x, y=y, l=l)
 ggplot(d, aes(x=x, y=y)) + geom_point() + facet_wrap(~l, scales="free_x") + 
   theme(axis.text.x=element_text(angle=90, hjust=1))

enter image description here

Andrie在这里的评论中建议可以在grid中手动完成此操作,但我不知道如何开始。

trans by 2020-08-10T20:03:49Z

r-在ggplot2中的X轴上格式化日期

我在非常非常困难的时间内使x轴看起来适合我的图形。

这是我的数据(通过as.Date生成):

df <- structure(list(Month = structure(1:12, .Label = c("2011-07-31", "2011-08-31", "2011-09-30", "2011-10-31", "2011-11-30", "2011-12-31", "2012-01-31", "2012-02-29", "2012-03-31", "2012-04-30", "2012-05-31", "2012-06-30"), class = "factor"), AvgVisits = c(6.98655104580674,7.66045407330464, 7.69761337479304, 7.54387561322994, 7.24483848458728, 6.32001400498928, 6.66794871794872, 7.207780853854, 7.60281201431308, 6.70113837397123, 6.57634103019538, 6.75321935568936)), .Names = c("Month","AvgVisits"), row.names = c(NA, -12L), class = "data.frame")

这是我尝试绘制的图表:

ggplot(df, aes(x = Month, y = AvgVisits)) + 
  geom_bar() +
  theme_bw() +
  labs(x = "Month", y = "Average Visits per User")

该图表可以正常工作-但是,如果我想调整日期格式,我相信我应该添加以下内容:as.Date

我正在尝试使日期标签为“ MMM-YYYY”

ggplot(df, aes(x = Month, y = AvgVisits)) + 
  geom_bar() +
  theme_bw() +
  labs(x = "Month", y = "Average Visits per User") +
  scale_x_date(labels = date_format("%m-%Y"))

当我绘制它时,我继续出现此错误:

as.Date

尽管对as.Dategeom_bar的格式进行了数小时的研究,但我无法修复它。 谁能解释我在做什么错?

编辑:作为后续思想:您可以使用日期作为因素,还是应该在日期列上使用as.Date

trans by 2020-08-10T19:11:41Z

我应该学习Fortran还是C ++来扩展R?

我使用相当大的数据集进行机器学习(它们仍然适合存储在内存中),并且我在R中编写了一些计算,但发现速度太慢。 因此,我想用我将从R中调用的编译代码替换程序的“关键部分”。我手头的一个示例问题是实现前向后向算法。

我的问题是我应该学习Fortran还是C ++来做到这一点? 我只需要使用数值向量或矩阵。 我主要感兴趣的是哪种语言更容易从R语言学习和使用R进行交互,我也不在乎哪种语言在我的简历上看起来更好。

我已经阅读了R扩展手册,并使用一些简单的Fortran和C ++代码对内联包进行了一些尝试。 我当前的印象是,尽管Rcpp软件包看起来也很有趣,但是Fortran95的学习更简单。 我目前了解R,Python和Matlab。

trans by 2020-08-10T12:24:30Z

r-查找特定日期属于哪个季节

我有一个日期向量,对于每个条目,我都想分配一个季节。 例如,如果日期在21.12之间。 和21.3,我会说是winter。到目前为止,我已经尝试了以下代码,但是无论年份如何,我都无法使它更通用。

my.dates <- as.Date("2011-12-01", format = "%Y-%m-%d") + 0:60
low.date <- as.Date("2011-12-15", format = "%Y-%m-%d")
high.date <- as.Date("2012-01-15", format = "%Y-%m-%d")

my.dates[my.dates <= high.date & my.dates >= low.date] 
 [1] "2011-12-15" "2011-12-16" "2011-12-17" "2011-12-18" "2011-12-19" "2011-12-20" "2011-12-21" "2011-12-22" "2011-12-23" "2011-12-24" "2011-12-25"
[12] "2011-12-26" "2011-12-27" "2011-12-28" "2011-12-29" "2011-12-30" "2011-12-31" "2012-01-01" "2012-01-02" "2012-01-03" "2012-01-04" "2012-01-05"
[23] "2012-01-06" "2012-01-07" "2012-01-08" "2012-01-09" "2012-01-10" "2012-01-11" "2012-01-12" "2012-01-13" "2012-01-14" "2012-01-15"

我试图格式化不带年份的日期,但是它不起作用。

ld <- as.Date("12-15", format = "%m-%d")
hd <- as.Date("01-15", format = "%m-%d")
my.dates[my.dates <= hd & my.dates >= ld] 
trans by 2020-08-10T11:47:48Z

statistics-如何在R中的线性回归模型中计算斜率的95%置信区间

这是R的入门统计练习:

使用rmr数据集,绘制代谢率与体重的关系图。 将线性回归模型拟合到该关系。 根据拟合模型,体重70公斤的预计代谢率是多少? 给出直线斜率的95%置信区间。

rmr数据集位于“ ISwR”包中。 看起来像这样:

> rmr
   body.weight metabolic.rate
1         49.9           1079
2         50.8           1146
3         51.8           1115
4         52.6           1161
5         57.6           1325
6         61.4           1351
7         62.3           1402
8         64.9           1365
9         43.1            870
10        48.1           1372
11        52.2           1132
12        53.5           1172
13        55.0           1034
14        55.0           1155
15        56.0           1392
16        57.8           1090
17        59.0            982
18        59.0           1178
19        59.2           1342
20        59.5           1027
21        60.0           1316
22        62.1           1574
23        64.9           1526
24        66.0           1268
25        66.4           1205
26        72.8           1382
27        74.8           1273
28        77.1           1439
29        82.0           1536
30        82.0           1151
31        83.4           1248
32        86.2           1466
33        88.6           1323
34        89.3           1300
35        91.6           1519
36        99.8           1639
37       103.0           1382
38       104.5           1414
39       107.7           1473
40       110.2           2074
41       122.0           1777
42       123.1           1640
43       125.2           1630
44       143.3           1708

我知道如何计算给定x的预测y,但是如何计算斜率的置信区间?

trans by 2020-08-10T11:39:01Z

在R中创建重复函数调用结果的向量

我有一个函数使用runif计算一些值,因此每次调用它时,结果都会略有不同。 我想计算多次调用该函数的结果的平均值。

为此,最好创建一个具有重复函数调用结果的向量

是否有一种简单的惯用方式来创建重复函数调用的向量? 我尝试

rep(my_function_call(), 10)

但是它只调用一次函数,并将结果重复10次。 我希望函数求值10次,并求结果的向量。

trans by 2020-08-10T05:20:41Z

r-将数据帧一起添加到for循环中

我有一个for loop,它在每次迭代后都会产生一个数据帧。 我想将所有数据帧附加在一起,但发现困难。 以下是我的正在尝试,请提出解决方法:

d = NULL
for (i in 1:7) {

  # vector output
  model <- #some processing

  # add vector to a dataframe
  df <- data.frame(model)

}

df_total <- rbind(d,df)
trans by 2020-08-10T04:35:23Z

无法在R上加载rJava

我希望在R x64 3.1.2中加载rJava。操作系统-Windows 8.1 64位

尽管安装似乎可以正常工作:

  > install.packages("rJava")
    Installing package into ‘C:/Users/sony/Documents/R/win-library/3.1’
    (as ‘lib’ is unspecified)
    --- Please select a CRAN mirror for use in this session ---
    trying URL 'http://cran.utstat.utoronto.ca/bin/windows/contrib/3.1/rJava_0.9-6.zip'
    Content type 'application/zip' length 758898 bytes (741 Kb)
    opened URL
    downloaded 741 Kb

package ‘rJava’ successfully unpacked and MD5 sums checked

The downloaded binary packages are in
        C:\Users\sony\AppData\Local\Temp\RtmpamYUH7\downloaded_packages

加载程序包时出现错误:

library(rJava)
Error in get(Info[i, 1], envir = env) : 
  lazy-load database 'C:/Users/sony/Documents/R/win-library/3.1/rJava/R/rJava.rdb' is corrupt
In addition: Warning message:
In get(Info[i, 1], envir = env) : internal error -3 in R_decompress1
Error: package or namespace load failed for ‘rJava’

我知道我正在使用R-64位,同时将环境变量Path设置为:

C:\Program Files\Java\jdk1.7.0_51\bin;C:\Program Files\Java\jre7\bin\server;

同样,Java版本也是64位。 我不了解错误背后的原因。

通过R手动设置时:

options(java.home="C:\\Program Files\\Java\\jdk1.7.0_51\\bin")
library(rJava)

我收到此错误:

enter image description here

Error : .onLoad failed in loadNamespace() for 'rJava', details:
  call: inDL(x, as.logical(local), as.logical(now), ...)
  error: unable to load shared object 'C:/Users/sony/Documents/R/win-library/3.1/rJava/libs/x64/rJava.dll':
  LoadLibrary failure:  The specified module could not be found.

Error: package or namespace load failed for ‘rJava’

有什么建议吗? 我卸载了某些旧版本的Java后,似乎出现了该错误。

其他:另外,我没有将JAVA_HOME设置为环境变量。

卸载并重新安装R似乎也不能解决问题。

删除较旧的Java版本也无济于事。

trans by 2020-08-10T02:48:32Z

R Shiny-动态将tabPanel添加到tabsetPanel(使用renderUI)

我正在开发一个闪亮的应用程序,其中使用了conditionPanel784,该应用程序是在用户输入某些特定输入时生成的。 因此,我想使用if()...785函数使tabsetPanel出现/消失。

我现在的难题是,conditionPanel元素的数量(参数3004122342172918918785)也取决于用户输入,方式是:有时我需要1个,其他时候我需要更多3004122342172918918786。

怎么做? 我尝试在3004122342172172918786参数中包含conditionPanel784或包含简单的if()...条件,但它(但并不令人惊讶...)不起作用。

trans by 2020-08-09T18:22:39Z

r-Rmarkdown yam中的多个作者和字幕

我正在尝试按照此pandoc示例将多个作者添加到yaml元数据块中的Rmarkdown文件中。 pdf将在RStudio(0.98.932版)中生成,但没有作者信息。

---
title:  'This is the title: it contains a colon'
author:
- name: Author One
  affiliation: University of Somewhere
- name: Author Two
  affiliation: University of Nowhere
date: "`r format(Sys.time(), '%d %B %Y')`"
tags: [nothing, nothingness]
abstract: |
  This is the abstract.

  It consists of two paragraphs.
output: pdf_document
---

我还想自定义标题并添加字幕。 可能?

trans by 2020-08-09T04:44:21Z

R中缩写的州名

我有一个带有状态全称的可变状态的大文件。 我想用州缩写(即“纽约”代表“纽约”)代替它。 是否有一种简单的方法(除了使用多个if-else命令之外)? 可能使用“替换”语句? 谢谢。

trans by 2020-08-08T20:40:52Z

r-如何很好地注释ggplot2(手动)

使用ggplot2我通常使用geom_text和类似position=jitter的注解。

但是-对于一个好的情节,我经常发现值得手动注释。 如下所示:

data2 <- structure(list(type = structure(c(5L, 1L, 2L, 4L, 3L, 5L, 1L, 
2L, 4L, 3L, 5L, 1L, 2L, 4L, 3L, 5L, 1L, 2L, 4L, 3L), .Label = c("EDS", 
"KIU", "LAK", "MVH", "NA*"), class = "factor"), value = c(0.9, 
0.01, 0.01, 0.09, 0, 0.8, 0.05, 0, 0.15, 0, 0.41, 0.04, 0.03, 
0.52, 0, 0.23, 0.11, 0.02, 0.64, 0.01), time = c(3L, 3L, 3L, 
3L, 3L, 6L, 6L, 6L, 6L, 6L, 15L, 15L, 15L, 15L, 15L, 27L, 27L, 
27L, 27L, 27L), year = c(2008L, 2008L, 2008L, 2008L, 2008L, 2007L, 
2007L, 2007L, 2007L, 2007L, 2007L, 2007L, 2007L, 2007L, 2007L, 
2006L, 2006L, 2006L, 2006L, 2006L)), .Names = c("type", "value", 
"time", "year"), row.names = c(1L, 3L, 4L, 5L, 6L, 7L, 9L, 10L, 
11L, 12L, 13L, 15L, 16L, 17L, 18L, 19L, 21L, 22L, 23L, 24L), class = "data.frame")
ggplot(data2, aes(x=time, y=value, group=type, col=type))+
geom_line()+
geom_point()+
theme_bw()+
annotate("text", x=6, y=0.9, label="this is a wrong color")+
annotate("text", x=15, y=0.6, label="this is a second annotation with a wrong color")

问题是,我无法获得文本注释颜色来匹配线条颜色。 我想我可以用手动比例尺解决此问题,但是我希望有更好的方法吗?

trans by 2020-08-08T20:39:27Z

包-R CMD ch中没有全局变量的可见绑定

当我使用类似subset之类的功能(使用列表元素的逐字名称作为参数)时,在检查软件包时注意到,我获得了注释“全局变量无可见绑定”。

例如,使用数据框:

foo <-data.frame(a = c(TRUE,FALSE,TRUE),b = 1:3)

我可以做一些愚蠢的事情,例如:

subset(foo,a)
transform(foo,a=b)

哪个工作按预期进行。 但是,R CMD中的R代码检查不理解它们是指元素,并且抱怨全局变量没有任何可见的绑定。

尽管这样做可以,但我真的不喜欢在包装中包含注释,而是希望它通过检查时完全没有错误,警告和注释。 我也不是很想重做我的代码。 有没有一种方法可以编写这些代码,从而使参数不指向全局变量?

trans by 2020-08-08T16:46:14Z

机器学习-计算R中的AUC?

给定分数向量和实际类别标签向量,如何为R语言或简单英语中的二进制分类器计算单数AUC度量?

“ AUC:更好的措施...”的第9页似乎需要了解类标签,而这是MATLAB中的一个我不理解的示例

R(Actual == 1))

因为R(不要与R语言混淆)被定义为向量而是用作函数?

trans by 2020-08-08T15:17:10Z

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 下一页 共53页