小编Sea*_*ean的帖子

我想完全删除刻面的标签以创造一种迷你效果,因为标签是不相关的观众,我能想到的最好的是:

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以便为"迷你图"留出更多空间吗？

或者,有没有更好的方法来获得像这样的大量二进制值时间序列的" 迷你 "效果？

有没有办法检测R在R内部使用的BLAS版本？我正在使用Ubuntu,我安装了几个BLAS版本 - 我只是不知道从R的角度看哪一个是"活跃的"!

我知道http://r.789695.n4.nabble.com/is-Rs-own-BLAS-td911515.html,Brian Ripley在2006年6月说过这是不可能的 - 但事情有变化吗？

我正在尝试分析使用surveymonkey创建的大型调查,该调查在CSV文件中有数百列,并且输出格式很难使用,因为标题会在两行上运行.

• 有没有人找到一种简单的方法来管理CSV文件中的标题,以便分析是可管理的？
• 其他人如何分析Surveymonkey的结果？

谢谢!

在R中,如何在write.csv(my.data.frame, file='test.csv')不重新读取它的情况下验证是否成功？它似乎没有返回任何东西.我正在考虑做一个file.exists('test.csv')然后在之前和之后抓一个时间戳write.csv()并检查文件上的时间戳是否介于两者之间？有什么建议？

当我尝试下面的任何一个时,轴刻度是1e + 03,1e + 06,1e + 09 - 是否可以获得漂亮的上标10 ^ 3,10 ^ 6,10 ^ 9而不需要手动标记？我似乎记得在过去自动获得这个.

qplot(1:10, 10^(1:10))+scale_y_log10()
qplot(1:10, 10^(1:10), log='y')

继上周我的查询中读取R-mismatched引号中的格式错误的csv后,这些相同的CSV文件也有嵌入的控制字符,例如ASCII 替换字符,即十进制26或0x1A.不幸的是readLines()似乎截断了这个字符的行,所以我很难匹配引号 - 除了丢失这些行中的后面的字段!

我试过readBin()但我无法读到这个文件.我担心我不能干净地把它读成R给你一个例子,我很难在R中创建这些.抱歉不能用一个干净的例子来证明.思考？

更新

现在我很困惑 - 当我使用代码时

 h3 <- paste('1,34,44.4,"', rawToChar(as.raw(c(as.integer(k1), 26, 65))), '",99')
 identical(readLines(textConnection(h3)), h3)

我觉得TRUE我觉得很惊讶!

更新2

 h3
[1] "1,34,44.4,\" HIJK\032A \",99"
> writeLines(h3, 'h3.txt')
> h3a <- readLines('h3.txt')
Warning message:
In readLines("h3.txt") : incomplete final line found on 'h3.txt'
> h3a
[1] "1,34,44.4,\" HIJK"

所以当来自a时,readLines()的反应会有所不同,textConnection()并且它会在SUB字符处静默截断.

如果它有所作为我会感到惊讶,但我在Windows-64上的2.15.2.

更新3

解决这个问题有些模糊的成功......

zb <- file('h3.txt', "rb")
tmp <- readBin(zb, raw(), size=1, n=400) # raw is always of …

在圣诞节之前,我之前曾在多个knitr文档中询问单一样式表.我现在知道我遗漏的关键搜索词是"子文档".

我创建了一个最小的降价示例,但我在创建的文档中获取了编织过程中的标题.我正在使用RStudio来处理它.我见过很多类似的问题,但没有一个答案可行.我已经尝试了各种组块的参数,到目前为止无济于事.

问题:如何摆脱最终针织文件中的标题？

# Master document

```{r master1, comment='', echo=FALSE, message=FALSE, warning=FALSE, results="asis", fig.width= 4., fig.height= 4, fig.cap= ""}
out = NULL
for (i in 1:3) {
  # lets create three different uniform data sets.
  v1 <- sort(rnorm(2))
  uvec <- runif(10, v1[[1]], v1[[2]])
  # now we want to 
  out <- c(out, knit_child('child.Rmd'))
}
cat(paste(out, collapse = '\n'))
```

子文档包含:

## child document `r i`

Here is a summary of the dataset:
```{r echo=FALSE}
summary(uvec)
``` …

我想把一个长长的传说分成两列,但我没有取得任何成功.这是我在其他地方找到的解决方案使用的代码geom='area',虽然它适用于我的其他绘图.我从下面的代码得到的情节如下:

那么如何Q1用两列中的图例进行绘图呢？

NVER <- 10
NGRID <- 20
MAT <- matrix(NA, nrow=NVER, ncol=NGRID)
gsd <- 0.1  # standard deviation of the Gaussians

verlocs <- seq(from=0, to=1, length.out=NVER)
thegrid <- seq(from=0, to=1, length.out=NGRID)
# create a mixture of Gaussians with modes spaced evenly on 0 to 1
# i.e. the first mode is at 0 and the last mode is at 1
for (i in 1:NVER) {
  # add the shape of gaussian i
  MAT[i,] <- dnorm(thegrid, …

我有一百万点和一个大的形状文件-8GB-这太大了,无法加载到我系统上的R内存中.形状文件是单层的,因此给定的x,y最多只能击中一个多边形 - 只要它不完全在边界上!每个多边形都标有一个severity-例如1,2,3.我在64位ubuntu机器上使用R,配备12GB内存.

什么是最简单的方法能够"标签"的数据帧的多边形severity,让我得到了data.frame一个额外的列,即x,y,severity？

我发现在Ubuntu下运行一些计算代码的性能相当差,这是在我用于科学计算的全新无头工作站机器上.我注意到在Ubuntu上运行一些稍微复杂的代码与我用于开发的旧Mac笔记本电脑上的速度有所不同.然而,我已经设法将它提炼成一个非常简单的例子,它仍然表现出比我的旧机器更少的改进:

#include <stdio.h>
#include <math.h>

int main() {
        double res = 0.0;
        for(int i=1; i<200000000; i++) {
                res += exp((double) 100.0/i);
        }
        printf("%lf", res);
        return(0);
}

现在,Mac是运行OS X 10.5的近5年的2.4GHz Core 2 Duo MacBook Pro,它运行此代码大约6.8秒.然而,在运行Ubuntu 11.10的全新3.4GHz Core i7戴尔上大约需要6.1秒!有人可以告诉我这里发生了什么,因为一台近5年的笔记本电脑在全新桌面工作站的10%以内是荒谬的吗？这更加荒谬,因为我可以通过监控工具看到Core i7涡轮增压到接近4GHz!

Mac编译:

gcc -o test test.c -std=gnu99 -arch x86_64 -O2

Ubuntu编译:

gcc -o test test.c -std=gnu99 -m64 -O2 -lm

谢谢,

路易