Tag: r

循环遍历R中的目录中的所有文件,应用多个命令

我需要在R中将一组命令应用到目录中的所有单个.txt文件(大约300)。 我对R不是很熟悉,所以我在网上看到有关循环的所有帮助都是令人困惑的,或者当您需要将多个命令应用于每个文件时,我无法弄清楚如何应用循环。 我需要应用到目录内的每个文件(系统发育树)的命令是(它使用R的ape库): testtree <- read.tree("tree123.txt") unrooted_tr <- unroot(testtree) write.tree(unrooted_tr, file="unrootedtree123.txt") 我如何应用一个循环,将这些命令应用到每个单独的.txt文件(使用R或在Unix命令行中)? 输出(例如unrootedtree123.txt)将需要为每个单独的文件具有不同的名称。 在此先感谢,丹妮。

使用ggplot2在一个canvas中的多个graphics

我正在试图合并两个ggplot2情节到基于这个表的一个: Type RatingA RatingB 1 One 3 36 2 Two 5 53 3 One 5 57 4 One 7 74 5 Three 4 38 6 Three 8 83 我想要制作两个散点图,分别是y轴上的评级平均值和x轴上的types。 这是我如何创build每个graphics: p1 <- ggplot(test, aes(x=reorder(Type, RatingA, mean), y=RatingA)) + stat_summary(fun.y="mean", geom="point") p2 <- ggplot(test, aes(x=reorder(Type, RatingB, mean), y=RatingB)) + stat_summary(fun.y="mean", geom="point") 由于p1和p2具有相同的x轴,所以我想让它们垂直排列。 我看着facet_align,但我找不到能做这项工作的东西。

使用ggplot2在地图上绘制条形图?

我想使用ggplot2为地图上的每个位置生成一个barplot,就像xm根据基础graphics和一些包所做的那样: Let R fly: Visualizing Export Data using R 这与在情节中embedded微型情节有关。 目前,我能做的最好的事情就是在抖动点图上匹配点的大小。 require(ggplot2) require(maps) #Get world map info world_map <- map_data("world") #Creat a base plot p <- ggplot() + coord_fixed() #Add map to base plot base_world <- p + geom_polygon(data=world_map, aes(x=long, y=lat, group=group)) #Create example data geo_data <- data.frame(long=c(20,20,100,100,20,20,100,100), lat=c(0,0,0,0,0,0,0,0), value=c(10,30,40,50,20,20,100,100), Facet=rep(c("Facet_1", "Facet_2"), 4), colour=rep(c("colour_1", "colour_2"), each=4)) […]

data.frame方法最有效的列表?

刚刚和同事谈了这件事情,我们认为值得看看SO土地上的人不得不说什么。 假设我有一个包含N个元素的列表,其中每个元素都是一个长度为X的向量。现在假设我想将其转换为一个data.frame。 与R中的大多数事物一样,有许多方法可以将谚语猫皮肤化,比如as.dataframe ,使用plyr软件包,与cbind结合使用do.call ,预先分配DF并填充它们等等。 提出的问题是当N或X(在我们的例子中是X)变得非常大时会发生什么。 有效的猫(特别是在记忆方面)的本质是否有一种猫皮方法明显优越?

在R / ggplot2中将math符号和下标与常规字母混合使用

我想在ggplot2绘制一个如下所示的标签: Value is $\sigma$, R^{2} = 0.6其中Value is是普通字体, $\sigma$是希腊小写sigma字母, R^{2} = 0.6为带有上标2的R ,接着是等号( = ),接着是0.6 。 这怎么可以在ggplot factor和参数的东西,如xlab,ylab R的xlab,ylab ? 谢谢。

创build一个逗号分隔的vector

我有一个数字向量,其中一个,我试图变成一个字符向量,每个元素之间用逗号分隔。 > one = c(1:5) > paste(as.character(one), collapse=", ") [1] "1, 2, 3, 4, 5" > paste(as.character(one), sep="' '", collapse=", ") [1] "1, 2, 3, 4, 5" 但是,我想要输出如下所示: "1", "2", "3", "4", "5" 我是否缺less粘贴function的一些参数? 帮帮我!?

查找两个数组之间的余弦相似度

我想知道在R中是否有内置函数可以find两个数组之间的余弦相似度(或余弦距离)? 目前,我执行了自己的function,但是我不禁想到R应该已经有了一个。

在`data.table`中使用dynamic列名

我想计算一个data.table中的几个列的每一个的平均值,由另一列分组。 我的问题类似于另外两个关于SO的问题( 一个和两个 ),但我不能将这些问题用于我的问题。 这里是一个例子: library(data.table) dtb <- fread(input = "condition,var1,var2,var3 one,100,1000,10000 one,101,1001,10001 one,102,1002,10002 two,103,1003,10003 two,104,1004,10004 two,105,1005,10005 three,106,1006,10006 three,107,1007,10007 three,108,1008,10008 four,109,1009,10009 four,110,1010,10010") dtb # condition var1 var2 var3 # 1: one 100 1000 10000 # 2: one 101 1001 10001 # 3: one 102 1002 10002 # 4: two 103 1003 10003 # 5: two […]

检查数字是否是整数

我很惊讶得知R没有提供一个方便的函数来检查数字是否是整数。 is.integer(66) # FALSE 帮助文件警告 : is.integer(x)不testingx包含整数! 为此,请使用round ,如例子中的函数is.wholenumber(x)所示。 该示例将此自定义函数作为“解决方法” is.wholenumber <- function(x, tol = .Machine$double.eps^0.5) abs(x – round(x)) < tol is.wholenumber(1) # is TRUE 如果我不得不写一个函数来检查整数,假设我没有阅读上面的评论,我会写一个函数, check.integer <- function(x) { x == round(x) } 我的方法会在哪里失败? 如果你在我的假设鞋子里,你的工作会是什么?

增加y轴上文本和标题之间的距离

Y轴标题看起来离轴文本太近了。 ggplot(mpg, aes(cty, hwy)) + geom_point() 我曾尝试更改theme()的许多参数的值,但似乎没有帮助。