小编MYa*_*208的帖子

我想用ggplot2重现以下基本图表.

以下是生成此图的R代码.

set.seed(12345)
Data <- matrix(data = rnorm(n = 30, mean = 0, sd = 1), nrow = 6, ncol = 5)
dimnames(Data) <- list(paste("G", 1:nrow(Data), sep = ""), paste("E", 1:ncol(Data), sep     = ""))
SVD <- svd(Data)
D <- diag(SVD$d[1:min(dim(Data))])
G <- SVD$u%*%sqrt(D)
E <- SVD$v%*%sqrt(D)

dimnames(G) <- list(rownames(Data))
dimnames(E) <- list(colnames(Data))

SVD.Values <- SVD$d
PC.No <- c(1:length(SVD.Values))
PC.SS <- SVD.Values^2
PC.Percent.SS <- PC.SS/sum(PC.SS)*100


library(grDevices)
con.hull.pos <- chull(G)
con.hull <- rbind(G[con.hull.pos, ], G[con.hull.pos[1], ])

getPerpPoints <- function(mat) { …

我想在xtable的每个页面上都有列名.

\documentclass{article}
\usepackage{graphicx}   
\usepackage{longtable}

\begin{document}

<<>>=
x <- matrix(rnorm(1000), ncol = 10)
x.big <- xtable(x,label='tabbig',caption='Example of longtable spanning several pages')
@

<<label = tabSEM.S, echo = FALSE, results = tex >>=
print(x.big,tabular.environment='longtable',floating=FALSE)
@

\end{document}

任何帮助将受到高度赞赏.谢谢

我想在R中为因子实验创建一个无约束的设计矩阵,下面的代码给出了我想要的矩阵.但是代码需要model.matrix针对每个因素以及拦截术语单独命令.我很好奇是否可以通过单个衬垫获得相同的结果.谢谢

y <- c(55, 56, 57, 53, 54, 55, 51, 52, 53, 61, 62, 63)
N <- gl(n = 2, k = 6, length = 2 * 6
        , labels = c("Low", "High")
        , ordered = FALSE)
P <- gl(n = 2, k = 3, length = 2 * 6
        , labels = c("Low", "High")
        , ordered = FALSE)
Data <- data.frame(y, N, P)

X <-
  cbind(
      model.matrix(object = y ~ 1,        data = Data)
    , model.matrix(object = …

我不知道如何更换NA与其他列的值以s列中的R使用dplyr.MWE在下面.

Letters <- LETTERS[1:5]
Char    <- c("a", "b", NA, "d", NA)
df1 <- data.frame(Letters, Char)
df1

library(dplyr]

df1 %>%
  mutate(Char1 = ifelse(Char != NA, Char, Letters))

     Letters Char Char1
1       A    a    NA
2       B    b    NA
3       C <NA>    NA
4       D    d    NA
5       E <NA>    NA

想知道如何将外部文件中的r代码片段包含到四开文档中。

1. 但是，以下代码可以 knitr::read_chunk()正常工作，具体取决于knitr。

   knitr::read_chunk("Ch02.R")
   

   Run Code Online (Sandbox Code Playgroud)

2. 以下使用四开代码的代码运行时没有任何错误，但无法按预期工作。

   #| eval: true
   #| output: false
   #| file: Ch02.R
   

   Run Code Online (Sandbox Code Playgroud)

想知道四开本中是否有任何替代方案knitr::read_chunk()。

我想知道如何在ggplot2中添加facet标题

ggplot(diamonds, aes(cut)) + geom_bar() +   facet_grid(color ~ clarity)

并更改图例标题

ggplot(diamonds, aes(cut, fill=cut)) + geom_bar() +   facet_grid(. ~ clarity)

谢谢你的帮助.

我想分割以下字符串

"ATextIWantToDisplayWithSpaces"

像这样

A Text I Want To Display With Spaces.

我试过这个代码 R

strsplit(x="ATextIWantToDisplayWithSpaces", split=[:upper:])

这会产生这个错误

Error: unexpected '[' in "strsplit(x="ATextIWantToDisplayWithSpaces", split=["

任何帮助将受到高度赞赏.谢谢

我可以使用nlme摘要中提取固定效果summary(fm1).但是如何获得分数却在苦苦挣扎Random effects:.

fm1 <- lme(distance ~ age, Orthodont, random = ~ age | Subject)
summary(fm1)
Linear mixed-effects model fit by REML
 Data: Orthodont 
       AIC      BIC    logLik
  454.6367 470.6173 -221.3183

Random effects:
 Formula: ~age | Subject
 Structure: General positive-definite, Log-Cholesky parametrization
            StdDev    Corr  
(Intercept) 2.3270340 (Intr)
age         0.2264278 -0.609
Residual    1.3100397       

Fixed effects: distance ~ age 
                Value Std.Error DF   t-value p-value
(Intercept) 16.761111 0.7752460 80 21.620377       0
age          0.660185 0.0712533 80  9.265333       0
 Correlation: 
    (Intr)
age -0.848 …

该knitr Chunk选项results = "hold"可以将输出放在后面Chunk Code.我想知道如何将knitr块代码和输出分成两个不同的knitrouts可能与标题Code和Output.在此先感谢您的帮助.

\documentclass{article} 
\begin{document}

<< label=Test, results = "hold" >>=
1:100
args(lm)
@ 
\end{document}

期望的输出

码

1:100
args(lm)

和

产量

 [1]   1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18
 [19]  19  20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36
 [37]  37  38  39  40 …

我想从 Github 中删除自定义域。请提供任何帮助。

实际上，我想要获取https://myaseen208.github.io/DialelAnalysisR/而不是http://www.myaseen208.com/DialelAnalysisR/，如下图所示。