小编Prr*_*dep的帖子

所以我对lm() - 命令感到有些困惑.我试了一下lm(x~y, mydata),并lm(y~x, mydata)和我有不同的输出.那么,哪个变量用作x,哪个变量用作y？我很抱歉问这样一个noob问题,但我不确定,我找不到任何解释该命令参数的东西!

我有一个圆柱数据框架结构.基于条件,我在运行时期间对数据帧进行了子集化.我观察到数据帧在子集化后被转换为向量.我已经使用该as.data.frame()函数实现了数据帧结构.

# random generation of the values    
df <- data.frame(a=sample(1:1000,100))
#checking the class of the object
class(df)
#dimensions
dim(df)
#[1] 100   1
#subsetting the data with a random value present in the df, here 547
df_sub <- df[-df$a==547,]
# checking the subset dataset class
class(df_sub)
#[1] "integer"

我想知道如何在不明确使用该as.data.frame()函数的情况下保留数据框架结构.

我一直在尝试使用 ggpairs（来自 GGally）作为探索性工具。我遇到了轴刻度标签相互覆盖的问题。如果这是直的 ggplot2，我会把它们放在 45 或 90 度角上，一切都会好起来，但我不知道如何为 ggpairs 做到这一点。在互联网上浏览让我感觉这可能是可以用包装来做的事情，但不知道如何继续。有人可以帮忙吗？

一些代码的结果很丑陋：

library(dplyr)
library(ggplot2)
library(GGally)
library(MASS)
cars <- dplyr::select(Cars93, Type, Price, MPG.city, MPG.highway, 
                  EngineSize, Origin)
ggpairs(cars)

我得到了丑陋的 ggpairs，左下角的标签非常混乱

谢谢您的关注。

如何可视化两个序列的完全比对？

library(Biostrings)
s1 <-DNAString("ACTTCACCAGCTCCCTGGCGGTAAGTTGATCAAAGGAAACGCAAAGTTTTCACTTCACCAGCTCCCTGGCGGTAAGTTGATCAAAGGAAACGCAAAGTTTTCAAGAAGACTTCACCAGCTCCCTGGCGGTAAGTTGATCAAAGGAAACGCAAAGTTTTCAAG")
s2 <-DNAString("GTTTCACTACTTCCTTTCGGGTAAGTAAATATATGTTTCACTACTTCCTTTCGGGTAAGTGTTTCACTACTTCCTTTCGGGTAAGTAAATATATAAATATATAAAAATATAATTTTCATCAAATATATAAATATATAAAAATATAATTTTCATCAAATATATAAAAATATAATTTTCATC")
pairwiseAlignment(s1,s2)

输出：

Global PairwiseAlignmentsSingleSubject (1 of 1)
pattern: [1] ACTTCACCAGCTCCCTGGCGGTAAGTTGATCAAAGGAAACGCAAAGT--TTTCAC---...CTTCACCAGCTCCCTGGCGGTAAGTTG-ATCAAAGG---AAACGCAAAGTTTTCAAG 
subject: [1] GTTTCACTACTTCCTTTCGGGTAAGTAAAT-ATATGTTTCACTACTTCCTTTCGGGTA...TATATAAATATATAAAAATATAATTTTCATCAAATATATAAAAATATAATTTTCATC 
score: -394.7115 

在这里，仅显示了部分对齐？您知道绘制或打印路线的任何现有功能吗？

我有两个数据:list包含271个值,listfull包含355个值(其中271个与列表中的相同).我想要的是在这两个文件中打印出84个不同的值.当我编写代码以打印"相同"值时,它们可以工作,但不适用于"不同情况".我试图找出原因,但不能.你可以帮我解决这个问题.谢谢

list<-read.table("C:\\Data\\list.txt", header=T)
listfull<-read.table("C:\\Data\\listfull.txt", header=T)



for (i in 1:271)
{
  for (j in 1:355)
  {
    if(list$Cow_ID[i]==listfull$Cow_ID[j])
    bo<-data.frame(listfull[j,])

   }
    write.table(bo,"C:\\Data\\store.txt",row.names = FALSE, dec = ".", na = "NA", sep = " ", append = TRUE, col.names = FALSE)
}

//以上代码有效!

//但不是以下内容:

for (i in 1:355)
{
for (j in 1:271)
{
if(listfull$Cow_ID[i]!=listfull$Cow_ID[j])
  bo<-data.frame(listfull[i,])

  }

  write.table(bo,"C:\\Data\\store.txt",row.names = FALSE, dec = ".", na = "NA", sep = " ", append = TRUE, col.names = FALSE)
}

这里list和listfull的样子如下:

Lisfull List
Cow_ID  Cow_ID
26       26 …