小编SHR*_*ram的帖子

# data 
set.seed (123)
xvar <- c(rnorm (1000, 50, 30), rnorm (1000, 40, 10), rnorm (1000, 70, 10))
yvar <-   xvar + rnorm (length (xvar), 0, 20)
myd <- data.frame (xvar, yvar)


# density plot for xvar
            upperp = 80   # upper cutoff
            lowerp = 30   # lower cutoff
            x <- myd$xvar
            plot(density(x))
            dens <- density(x)
            x11 <- min(which(dens$x <= lowerp))
            x12 <- max(which(dens$x <= lowerp))
            x21 <- min(which(dens$x > upperp))
            x22 <- max(which(dens$x > upperp))
            with(dens, polygon(x = c(x[c(x11, …

我正在使用Rmarkdown和 Knitr使用Rstudio.

以下是打印脚本和输出到html.

```{r}
summary(cars)
```

但是以下只会打印输出,即嵌入图.

```{r, echo=FALSE}
plot(cars)
```

我的情况与上面不同,我想呈现脚本但不应该在html中运行,因为这将需要很长时间(小时,如果不是几天)来运行.所以我刚刚做了评论标志.

```{r}
#summary(cars)
```

但我需要一个更好的方法来做到这一点 - 有没有更好的方式呈现脚本而不运行它.

我正在编写一个函数,需要检查(和哪个!)列(变量)是否具有所有缺失值(NA,<NA>).以下是该功能的片段:

test1 <- data.frame (matrix(c(1,2,3,NA,2,3,NA,NA,2), 3,3))
test2 <- data.frame (matrix(c(1,2,3,NA,NA,NA,NA,NA,2), 3,3))

na.test <-  function (data) {
  if (colSums(!is.na(data) == 0)){
      stop ("The some variable in the dataset has all missing value,
     remove the column to proceed")
      }
      }
na.test (test1)

Warning message:
In if (colSums(!is.na(data) == 0)) { :
  the condition has length > 1 and only the first element will be used

Q1:为什么上述错误和任何修复？

Q2:有没有办法找到哪一列都有NA,例如输出列表(变量名或列号)？

我有两个文件夹,folder1并folder2与周围的200个文件每个要么*rda或*R.我想从两个目录中读取所有文件和数据集.我怎样才能做到这一点？

路径:

folder1:  C:\folder1
folder2:  C:\folder2 

我的审判

setwd("C:/folder1")
myls <- ls() # do work as this will only list that are already loaded in the system 
setwd("C:/folder2")
myls2 <- ls()
myls  # do work as this will only list that are already loaded in the system 

我知道这是一个简单的问题,但我没有任何答案.

我想从X轴移除标签和轴,但是添加新的刻度.

plot(1:10, ylab = "")
at1 <- seq(1, 10, 0.1)
axis(side = 1, at = at1, labels = FALSE)

我无法摆脱y标签.

我是外行人,unix和sofar我在Windows中使用R. 例如,我在R会话中输入以下内容(在R gui中).

# this is a my funny example script 
X <- 1:10
Y <- 21:30
plot(X, Y)
myfun <- function (x){
              x1 <- x^0.2
              return (x1)
             }
myfun(X)

在两种情况下,如何在unix shell中实现这一点 -

(1)直接在命令行中通过interpeter(2)创建脚本并运行脚本.

请考虑我是unix的外行人.

我正在尝试从以下数据创建结构图:

mydf <- data.frame ( group = rep (1:5, each = 20), z = rnorm (20, 10, 1),
                  x = c(rnorm (20, 2, 0.5), rnorm (20, 2, 0.5),
          rnorm (20, 9, 0.5), rnorm (20, 9, 0.5),rnorm (20, 5, 0.5)),
       y = c(rnorm (20, 2, 0.5), rnorm (20, 9, 0.5), rnorm (20, 2, 0.5),
      rnorm (20, 9, 0.5), rnorm (20, 2, 0.5)))

means <- aggregate(. ~ group, data = mydf, mean)
gmx <-mean (mydf$x)
gmy <- mean (mydf$y)

library(ggplot2)
ggplot(mydf, aes(x, …

我试图将标记放在条形内以创建条形图(水平或垂直)以提供像刻度的视图.

只是一个小例子:

myd <- data.frame (names = c("A", "B", "C", "D"), height = c(2.1, 3.5, 3.5,1.5))
require(ggplot2)
c <- ggplot(myd, aes(factor(names), height, fill = names))
c + geom_bar()
c + geom_bar() + coord_flip()

有没有一种简单的方法或包可以实现这一目标？

我有一个数据集,重复次数不等.我希望通过删除那些不完整的条目(即复制小于最大值)来对数据进行子集化.只是一个小例子:

set.seed(123)
mydt <- data.frame (name= rep ( c("A", "B", "C", "D", "E"), c(1,2,4,4, 3)), 
                   var1 = rnorm (14, 3,1), var2 = rnorm (14, 4,1))
 mydt
       name     var1     var2
1     A 2.439524 3.444159
2     B 2.769823 5.786913
3     B 4.558708 4.497850
4     C 3.070508 2.033383
5     C 3.129288 4.701356
6     C 4.715065 3.527209
7     C 3.460916 2.932176
8     D 1.734939 3.782025
9     D 2.313147 2.973996
10    D 2.554338 3.271109
11    D 4.224082 3.374961
12    E 3.359814 2.313307
13    E 3.400771 …

我想用多边形重绘R中的任何草图.但我需要一项艰苦的工作才能找到坐标(每个点的X或Y值).

是否有任何R包或其他软件可以从图像中找到坐标？因此,输出将是具有相当高的点数(x,y)的数据,以便可以重新创建该图.

编辑:........

以下是示例:

(1)地图轮廓例如:

第二个例子:

(2)对象形状示例: