Tag: r

条件长度> 1,只有第一个元素将被使用

我有一个dataframe,旅程: > head(trip.mutations) Ref.y Variant.y 1 TC 2 GC 3 AC 4 TC 5 CA 6 GA 我想添加第三列,mutType,遵循这些规则: for (i in 1:nrow(trip)) { if(trip$Ref.y=='G' & trip$Variant.y=='T'|trip$Ref.y=='C' & trip$Variant.y=='A') { trip[i, 'mutType'] <- "G:C to T:A" } else if(trip$Ref.y=='G' & trip$Variant.y=='C'|trip$Ref.y=='C' & trip$Variant.y=='G') { trip[i, 'mutType'] <- "G:C to C:G" } else if(trip$Ref.y=='G' & trip$Variant.y=='A'|trip$Ref.y=='C' & […]

在R中使用ggplot2,我如何在不同区域使graphics的背景变成不同的颜色?

我使用ggplot2软件包在R中制作一个简单的条形图。 我不想将灰色的默认值分成五个区域,每个区域都是不同的(但是同样是低调的)颜色。 我如何做到这一点? 更具体地说,我想要五个颜色区域从0-25,25-45,45-65,65-85和85-100运行,其中颜色代表劣于青铜,青铜,银,金和铂分别。 一个配色scheme的build议也非常欢迎。 感谢您的帮助。

R + ggplot:绘制多个页面

我试图写一个函数,绘制一个ggplot facet_wrap图在多个页面。 这只是一个黑客,因为这个function似乎是在ggplot2function待办事项列表。 我做了一些小的计算,以find我需要的页数,我的data.frame的行数, data.frame需要我等。我非常有信心,这一切工作。 pdf(filename) for (i in seq(num_pages)){ slice = seq(((i-1)*num_rows)+1,(i*num_rows)) slice = slice[!(slice > nrow(df.merged))] df.segment=df.merged[slice,] p <- ggplot(df.segment, aes(y=mean,x=phenotype)) p <- p + geom_bar(stat="identity",fill="white",colour="black") p + facet_wrap("ID",scales="free_y",ncol=n_facets,nrow=n_facets) } dev.off() 我的问题是,通过在pdf()和dev.off()函数之间将它们全部封装在for循环中,for循环似乎并不等待ggplot完成它的工作,非常快速地通过其循环并输出无效的PDF。 如果我设置i = 1 ,启动pdf() ,在for循环中运行上面的代码,然后设置i=2 ,然后运行代码,等等,直到我感到无聊( i=3 ),然后closures设备由此产生的PDF是辉煌的。 有没有一种方法,我可以让循环等待最后一行完成绘图,然后再进入下一个迭代?

Rdataframe中的行数基于组

我有这样的数据框在R : ID MONTH-YEAR VALUE 110 JAN. 2012 1000 111 JAN. 2012 2000 . . . . 121 FEB. 2012 3000 131 FEB. 2012 4000 . . . . 因此,每年每个月有n行,他们可以以任何顺序(意思是他们都不是连续的,并在rest)。 我想计算每个MONTH-YEAR有多less行,即JAN有多less行。 2012年,FEB有多less。 2012年等。 像这样的东西: MONTH-YEAR NUMBER OF ROWS JAN. 2012 10 FEB. 2012 13 MAR. 2012 6 APR. 2012 9 我试图做到这一点: n_row <- nrow(dat1_frame %.% […]

ggmap错误:GeomRasterAnn是用不兼容的ggproto版本构build的

我正在使用ggmap,并得到下面的错误: Error: GeomRasterAnn was built with an incompatible version of ggproto. Please reinstall the package that provides this extension. 我已经安装了ggmap(2.6.1)和ggplot2(2.2.0)的最新版本,但仍然有同样的错误。

分割文件名称

我有这样的文件的名称: name1.csv ,我想提取这个string的两个子string。 一个将name1存储在一个variables中,另一个存储扩展名csv ,而不在另一个variables中存储点。 我一直在寻找是否有像Java的indexOffunction,允许做这种操作,但我什么也没有find。 任何帮助?

如何在R中读取utf-8格式的数据?

我的系统:win7 + R-3.0.2。 > Sys.getlocale() [1] "LC_COLLATE=Chinese (Simplified)_People's Republic of China.936;LC_CTYPE=Chinese (Simplified)_People's Republic of China.936;LC_MONETARY=Chinese (Simplified)_People's republic of China.936;LC_NUMERIC=C;LC_TIME=Chinese (Simplified)_People's Republic of China.936" 在微软记事本中保存了两个相同内容的文件: 一个保存为ansi格式, 另一个保存为utf8格式,数据为马来西亚航空公司M370的死亡名称。 或者你可以这样创build文件。 1)将数据复制到Microsoft记事本。 乘客姓名,性别,出生dateHuangTianhui,男,1948/05/28姜翠云,女,1952/03/27李红晶,女,1994/12/09 2)在记事本中以ansi格式保存为test.ansi。 3)在记事本中以utf-8格式保存为test.utf8。 read.table("test.ansi",sep=",",header=TRUE) #can work fine read.table("test.utf8",sep=",",header=TRUE) #can't work 然后,我把编码设置为utf-8。 options(encoding="utf-8") read.table("test.utf8",sep=",",header=TRUE,encoding="utf-8") In read.table("test.utf8", sep = ",",header=TRUE,encoding = "utf-8") : invalid input found on input connection […]

在向量中查找多个元素

假设我有以下vector: > x <- sample(1:10,20,replace=TRUE) > x [1] 8 6 9 9 7 3 2 5 5 1 6 8 5 2 9 3 5 10 8 2 我怎样才能find哪些元素是8或9?

一次重塑多个值

我有一个很长的数据集,我想扩大,我很好奇,如果有一种方法,在R一步使用reshape2或tidyr包完成这一切。 dataframedf如下所示: id type transactions amount 20 income 20 100 20 expense 25 95 30 income 50 300 30 expense 45 250 我想要得到这个: id income_transactions expense_transactions income_amount expense_amount 20 20 25 100 95 30 50 45 300 250 我知道我可以通过例如reshape2获得部分path: dcast(df, id ~ type, value.var="transactions") 但是,是否有办法一次性重塑整个DF,同时解决“交易”和“金额”两个variables? 理想情况下,新的更适合的列名称?

在一张图中绘制多条线

尝试使用ggplot将多行图表绘制成一个图表,但不知道如何使用我的数据集。 不确定是否需要更改数据结构(转置?) 数据看起来像这样: Company 2011 2013 Company1 300 350 Company2 320 430 Company3 310 420 我也试过换位: Year Company1 Company2 Company3 2011 300 320 310 2013 350 430 420 为此,我可以绘制一个值的使用; ggplot(data=df, aes(x=Year, y=Company1)) + geom_line(colour="red") + geom_point(colour="red", size=4, shape=21, fill="white") 但是我不知道如何组合所有的公司,因为我再也没有一个对象“公司”了。 有什么build议么?