Questions tagged «r»

R是一种免费的开源编程语言和软件环境,用于统计计算,生物信息学,可视化和通用计算。请提供最少且可复制的示例以及所需的输出。使用`dput()`作为数据,并通过`library()`调用指定所有非基本包。不要为数据或代码嵌入图片,而应使用缩进的代码块。对于统计相关的问题,请使用https://stats.stackexchange.com。

3
向ggplot2线图添加图例
我对ggplot2中的图例有疑问。我设法在同一张图中绘制三条线,并希望使用所用的三种颜色添加图例。这是使用的代码 library(ggplot2) require(RCurl) link<-getURL("https://dl.dropbox.com/s/ds5zp9jonznpuwb/dat.txt") datos<- read.csv(textConnection(link),header=TRUE,sep=";") datos$fecha <- as.POSIXct(datos[,1], format="%d/%m/%Y") temp = ggplot(data=datos,aes(x=fecha, y=TempMax,colour="1")) + geom_line(colour="red") + opts(title="TITULO") + ylab("Temperatura (C)") + xlab(" ") + scale_y_continuous(limits = c(-10,40)) + geom_line(aes(x=fecha, y=TempMedia,colour="2"),colour="green") + geom_line(aes(x=fecha, y=TempMin,colour="2"),colour="blue") + scale_colour_manual(values=c("red","green","blue")) temp 和输出 我想添加一个图例,其中使用了三种颜色和变量名称(TempMax,TempMedia和TempMin)。我努力了 scale_colour_manual 但找不到确切的方法。 不幸的是,原始数据已从链接站点中删除,无法恢复。但是它们来自具有这种格式的气象数据文件 "date","Tmax","Tmin","Tmed","Precip.diaria","Wmax","Wmed" 2000-07-31 00:00:00,-1.7,-1.7,-1.7,-99.9,20.4,20.4 2000-08-01 00:00:00,22.9,19,21.11,-99.9,6.3,2.83 2000-08-03 00:00:00,24.8,12.3,19.23,-99.9,6.8,3.87 2000-08-04 …
143 r  ggplot2  legend  r-faq 


14
确定列中NA值的数量
我想计算NA数据框列中的值数量。假设我的数据框被调用df,而我正在考虑的列的名称为col。我提出的方法如下: sapply(df$col, function(x) sum(length(which(is.na(x))))) 这是一种很好/最有效的方法吗?
142 r  dataframe 

17
用最新的非NA值替换NA
在data.frame(或data.table)中,我想用最接近的先前非NA值“填充” NA。data.frame以下是使用向量(而不是)的简单示例: > y <- c(NA, 2, 2, NA, NA, 3, NA, 4, NA, NA) 我想要一个fill.NAs()允许我构造的函数yy: > yy [1] NA NA NA 2 2 2 2 3 3 3 4 4 我需要对许多(总计〜1 Tb)小型data.frames(〜30-50 Mb)重复此操作,其中所有条目都是NA。解决这个问题的好方法是什么? 我煮的丑陋的解决方案使用此功能: last <- function (x){ x[length(x)] } fill.NAs <- function(isNA){ if (isNA[1] == 1) { isNA[1:max({which(isNA==0)[1]-1},1)] …
141 r  data.table  zoo  r-faq 

6
在R中如何使用“ <<-”(作用域分配)?
我刚读完R前言中有关范围定义的内容,并对&lt;&lt;-分配感到非常好奇。 手册显示了一个(非常有趣的)示例&lt;&lt;-,我认为我理解了。我仍然缺少的是什么时候有用的上下文。 因此,我很乐意从您那里读到一些示例(或示例的链接),这些示例说明何时使用&lt;&lt;-有趣/有用。使用它可能会有什么危险(看起来很容易松动),以及您可能想共享的任何提示。

12
计算每个唯一值的出现次数
假设我有: v = rep(c(1,2, 2, 2), 25) 现在,我想计算每个唯一值出现的次数。unique(v) 返回唯一值是多少,但不是多少。 &gt; unique(v) [1] 1 2 我想要能给我的东西 length(v[v==1]) [1] 25 length(v[v==2]) [1] 75 但作为更一般的单行代码:)接近(但不太完全)的东西是这样的: #&lt;doesn't work right&gt; length(v[v==unique(v)])
140 r  count  unique 

5
R中的“ throw”的等效项
如何“抛出” R中的错误?我有一个函数,它接受一个数据框和一些列名,并对其进行填充。如果列不存在,我希望该函数停止并停止所有依赖于此的函数。 我已经看过recover并browse和traceback,但同时,他们似乎是接近,但不是我所期待的。
139 r  error-handling 

6
没有轴,图例等的ggplot2图
我想使用生物导体的六边形(我可以这样做)来生成填充整个(png)显示区域的图-无轴,无标签,无背景,无螺母。
139 r  ggplot2 



4
“ * apply”家庭真的没有向量化吗?
因此,我们习惯对每个R的新用户说“ apply不是矢量化的,请查看Patrick Burns R Inferno Circle 4 ”,其中说(我引用): 常见的反射是在apply系列中使用功能。这不是 向量化,而是循环隐藏。apply函数的定义中包含一个for循环。lapply函数掩盖了循环,但是执行时间往往近似等于显式的for循环。 的确,快速查看apply源代码可以发现循环: grep("for", capture.output(getAnywhere("apply")), value = TRUE) ## [1] " for (i in 1L:d2) {" " else for (i in 1L:d2) {" 到目前为止还可以,但是看看lapply还是vapply实际上可以看到完全不同的图片: lapply ## function (X, FUN, ...) ## { ## FUN &lt;- match.fun(FUN) ## if (!is.vector(X) || is.object(X)) ## …

9
为组合的ggplots添加通用图例
我有两个ggplots,它们与水平对齐grid.arrange。我浏览了很多论坛帖子,但是我尝试的所有内容似乎都是现在已更新并命名为其他名称的命令。 我的数据看起来像这样; # Data plot 1 axis1 axis2 group1 -0.212201 0.358867 group2 -0.279756 -0.126194 group3 0.186860 -0.203273 group4 0.417117 -0.002592 group1 -0.212201 0.358867 group2 -0.279756 -0.126194 group3 0.186860 -0.203273 group4 0.186860 -0.203273 # Data plot 2 axis1 axis2 group1 0.211826 -0.306214 group2 -0.072626 0.104988 group3 -0.072626 0.104988 group4 -0.072626 0.104988 …
138 r  ggplot2  legend  gridextra 

11
R无法做到的MATLAB能做什么?[关闭]
按照目前的情况,这个问题不适合我们的问答形式。我们希望答案会得到事实,参考或专业知识的支持,但是这个问题可能会引起辩论,争论,民意调查或扩展讨论。如果您认为此问题可以解决并且可以重新提出,请访问帮助中心以获取指导。 8年前关闭。 已锁定。该问题及其答案被锁定,因为该问题是题外话,但具有历史意义。它当前不接受新的答案或互动。 我经常听到人们抱怨MATLAB许可证的价格昂贵。然后,我不知道为什么他们不只是用八度或[R 。但是后者对吗?可以使用R取代MATLAB吗?

3
如何在R中保存data.frame?
我在R中创建了一个不是很大的data.frame,但是要花很多时间才能构建。我要将其另存为文件,然后可以再次在R中打开该文件吗?
137 r  dataframe 

14
ggplot2中具有边缘直方图的散点图
是否有办法像下面的示例中那样用边际直方图创建散点图ggplot2?在Matlab中,它是scatterhist()函数,R也存在等效项。但是,我还没有在ggplot2上看到它。 我通过创建单个图形开始尝试,但是不知道如何正确排列它们。 require(ggplot2) x&lt;-rnorm(300) y&lt;-rt(300,df=2) xy&lt;-data.frame(x,y) xhist &lt;- qplot(x, geom="histogram") + scale_x_continuous(limits=c(min(x),max(x))) + opts(axis.text.x = theme_blank(), axis.title.x=theme_blank(), axis.ticks = theme_blank(), aspect.ratio = 5/16, axis.text.y = theme_blank(), axis.title.y=theme_blank(), background.colour="white") yhist &lt;- qplot(y, geom="histogram") + coord_flip() + opts(background.fill = "white", background.color ="black") yhist &lt;- yhist + scale_x_continuous(limits=c(min(x),max(x))) + opts(axis.text.x = theme_blank(), axis.title.x=theme_blank(), …

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.