Questions tagged «r»

我想完全删除刻面的标签，以创建一种迷你图效果，对于观众来说，这些标签无关紧要，我能想到的最好的方法是： library(MASS) library(ggplot2) qplot(week,y,data=bacteria,group=ID, geom=c('point','line'), xlab='', ylab='') + facet_wrap(~ID) + theme(strip.text.x = element_text(size=0)) 因此，我可以完全摆脱（现在为空白）strip.background以便为“火花线”留出更多空间吗？ 或者，对于像这样的大量二进制值时间序列，是否有更好的方法来获得这种“迷你图”效果？

84 r  graphics  ggplot2 

当我将矩阵子集分配到单个列时，结果是数字类而不是矩阵类（即，将myMatrix [，5]子集分配到第五列）。是否有一种紧凑的方法可以将子集划分为单个列，维护矩阵格式以及维护行/列名称，而无需执行类似以下的复杂操作： matrix( myMatrix[ , 5 ] , dimnames = list( rownames( myMatrix ) , colnames( myMatrix )[ 5 ] )

84 r  matrix  subset  r-faq 

一堆中间文件（称为temp1，temp2等）阻塞了我的内存。是否有可能从内存中删除他们不这样做rm(temp1)，rm(temp2)一次一个？我试过了rm(list(temp1, temp2, etc.))，但这似乎不起作用。

84 r 

如果出于参数考虑，我想要Python中10个长度的向量的最后五个元素，则可以在范围索引中使用“-”运算符，因此： >>> x = range(10) >>> x [0, 1, 2, 3, 4, 5, 6, 7, 8, 9] >>> x[-5:] [5, 6, 7, 8, 9] >>> 在R中执行此操作的最佳方法是什么？有没有一种比我当前使用长度（）函数的方法更简洁的方法？ > x <- 0:9 > x [1] 0 1 2 3 4 5 6 7 8 9 > x[(length(x) - 4):length(x)] [1] 5 6 …

84 r  indexing 

是否有任何简单的方法可以像matlab一样将R脚本编译成独立.exe文件？

84 r 

我想知道如何rmarkdown在同一个文档中生成同时具有纵向和横向布局的pdf。如果有一个纯粹的rmarkdown选择，那将比使用乳胶更好。 这是一个小的可复制示例。首先，.Rmd在RStudio中渲染此图像（按“编织PDF”按钮）将生成pdf，所有页面均采用横向布局： --- title: "All pages landscape" output: pdf_document classoption: landscape --- ```{r} summary(cars) ``` \newpage ```{r} summary(cars) ``` 然后尝试创建混合纵向和横向布局的文档。中的基本设置YAML是根据此处的“包含”部分完成的。该in_header文件“header.tex”只包含\usepackage{lscape}，建议一揽子knitr景观布局在这里。该.tex文件与该文件位于同一目录中.Rmd。 --- title: "Mixing portrait and landscape" output: pdf_document: includes: in_header: header.tex --- Portrait: ```{r} summary(cars) ``` \newpage \begin{landscape} Landscape: ```{r} summary(cars) ``` \end{landscape} \newpage More portrait: ```{r} summary(cars) ``` 但是，此代码导致错误： …

84 r  pdf  rstudio  knitr  r-markdown 

我已经阅读了该问题的答案，它们非常有帮助，但是我特别需要R中的帮助。 我在R中有一个示例数据集，如下所示： x <- c(32,64,96,118,126,144,152.5,158) y <- c(99.5,104.8,108.5,100,86,64,35.3,15) 我想对这些数据拟合模型，以便y = f(x)。我希望它是一个三阶多项式模型。 如何在R中做到这一点？ 另外，R可以帮助我找到最佳的拟合模型吗？

83 r  curve-fitting  data-analysis  polynomial-math 

Stargazer为lm（和其他）物体产生非常漂亮的乳胶表。假设我已经以最大可能性拟合了模型。我希望观星者为我的估算出一张类似lm的桌子。我怎样才能做到这一点？ 尽管有点麻烦，但是一种方法可能是创建一个包含我的估算值的“伪” lm对象-我认为，只要summary（my.fake.lm.object）有效，它就可以工作。这容易做到吗？ 一个例子： library(stargazer) N <- 200 df <- data.frame(x=runif(N, 0, 50)) df$y <- 10 + 2 * df$x + 4 * rt(N, 4) # True params plot(df$x, df$y) model1 <- lm(y ~ x, data=df) stargazer(model1, title="A Model") # I'd like to produce a similar table for the model …

83 r  optimization  lm  stargazer 

我有一个data.frame，其中包含一些具有所有NA值的列，如何从data.frame中删除它们。 我可以使用该功能吗 na.omit(...) 指定一些其他参数？

83 r  dataframe  na 

三年多以前才问过这个问题。给出了一个答案，但是我在解决方案中发现了一个小故障。 下面的代码在R中。我已将其移植到另一种语言，但是已经在R中直接测试了原始代码，以确保问题不在于我的移植。 sunPosition <- function(year, month, day, hour=12, min=0, sec=0, lat=46.5, long=6.5) { twopi <- 2 * pi deg2rad <- pi / 180 # Get day of the year, e.g. Feb 1 = 32, Mar 1 = 61 on leap years month.days <- c(0,31,28,31,30,31,30,31,31,30,31,30) day <- day + cumsum(month.days)[month] leapdays …

83 r  math  geometry  astronomy  azimuth 

我有大约100个igraph对象的列表，其中典型的对象具有大约700个顶点和3500条边。 我想确定更多可能存在联系的顶点组。然后，我的计划是使用混合模型来预测使用顶点和组属性的组内连接顶点数。 某些人可能想对我项目的其他方面做出回应，这可能很棒，但是我最感兴趣的是有关igraph中用于将顶点分组的函数的信息。我遇到过这些社区检测算法，但是我不确定它们的优缺点，或者不确定某些其他功能是否适合我的情况。我在这里也看到了链接，但它们并非特定于igraph。谢谢你的建议。

83 r  igraph 

当我取消列出日期列表时，会将它们重新转换为数字。那是正常的吗？除了重新申请as.Date以外，还有其他解决方法吗？ > dd <- as.Date(c("2013-01-01", "2013-02-01", "2013-03-01")) > class(dd) [1] "Date" > unlist(dd) [1] "2013-01-01" "2013-02-01" "2013-03-01" > list(dd) [[1]] [1] "2013-01-01" "2013-02-01" "2013-03-01" > unlist(list(dd)) [1] 15706 15737 15765 这是一个错误吗？

83 r 

在R中建立模型时，如何保存模型规格，以便可以在新数据上重用它？假设我在历史数据上建立了逻辑回归，但直到下个月才有新的观察结果。最好的方法是什么？ 我考虑过的事情： 保存模型对象并在新会话中加载 我知道某些模型可以使用PMML导出，但是对于导入PMML并没有真正了解 简单来说，当您需要在新的会话中使用模型时，我试图了解您的工作。 提前致谢。

82 r  models 

我在R中有一个大约10,000个元素的列表。假设我只选择元素5、7和9。我不确定如果没有for循环该怎么做。 我想做类似的事情，mylist[[c(5,7,9]]但这不起作用。我也尝试过该lapply功能，但也无法使其正常工作。

82 r  list  subset 

我想使R中的python字典等效。基本上，在python中，我有： visited = {} if atom_count not in visited: Do stuff visited[atom_count] = 1 这个想法是，如果我看到了那个特定的atom_count，我就访问了[atom_count] =1。因此，如果我再次看到那个atom_count，我就不会“做东西”。Atom_Count是一个整数。 谢谢！

82 python  r  dictionary