标签: pairwise

该包ggsignif对于快速、轻松地指示ggplot图表中的显着比较非常有用。但是，该comparisons调用需要手动键入要比较的每对值。

例如。

library(ggplot2)
library(ggsignif)

data(iris)

ggplot(iris, aes(x=Species, y=Sepal.Length)) + 
  geom_boxplot() +
  geom_signif(comparisons = list(c("versicolor", "virginica"),c('versicolor','setosa')), 
              map_signif_level=TRUE)

我想知道如何通过立即引用所有可能的组合来避免这种情况？例如，expand.grid(x = levels(iris$Species), y = levels(iris$Species))，给出所有组合

           x          y
1     setosa     setosa
2 versicolor     setosa
3  virginica     setosa
4     setosa versicolor
5 versicolor versicolor
6  virginica versicolor
7     setosa  virginica
8 versicolor  virginica
9  virginica  virginica

但如何让这个被接受geom_signif(comparisons=...呢？

软件包信息可在此处获取https://cran.r-project.org/web/packages/ggsignif/index.html

我在 TensorFlow 中有一个 N 维的一维向量，

如何构造成对平方差的总和？

例子

输入向量
[1,2,3]
输出 6
计算为

(1-2)^2+(1-3)^2+(2-3)^2.

如果我将输入作为 N-dim 向量 l，则输出应为 sigma_{i,j}((l_i-l_j)^2)。

添加的问题：如果我有一个 2d 矩阵并且想对矩阵的每一行执行相同的过程，然后对所有行的结果求平均值，我该怎么做？非常感谢！

我正在尝试为垂直和水平滚动方向创建两个可观察的滚动事件。

我尝试使用pairwise()和bufferCount(2,1)运算符从水平滚动事件中过滤垂直滚动事件，但问题是获取prev.scrollTop和 的重复值curr.scrollTop

import { Component, ViewChild, AfterViewInit, ElementRef } from '@angular/core';
import { fromEvent } from 'rxjs';
import { pairwise, tap, filter } from 'rxjs/operators';

@Component({
  selector: 'my-app',
  templateUrl: './app.component.html',
  styleUrls: ['./app.component.css']
})
export class AppComponent implements AfterViewInit {

  @ViewChild('scrollable', {static: false}) scrollable: ElementRef;


  ngAfterViewInit() {

    fromEvent(this.scrollable.nativeElement, 'scroll').pipe(
      pairwise(),
      tap(([prev, curr]) => console.log(prev.target.scrollTop, curr.target.scrollTop)),
      filter(([prev, curr]) => prev.target.scrollTop !== curr.target.scrollTop),
      tap((e) => console.log(e)) // <=    Never reached
    ).subscribe();

  }

}

有任何想法吗？

堆栈闪电复制

我一直在尝试在每次排放中获得新旧价值。我已经看到了使用pairwise或的选项，bufferCount但它们不允许保留第一个值。

目标将来自：

---1---2---3---4---5---

至：

---null,1---1,2---2,3---3,4---4,5---

有任何想法吗？

x = 1:3
y = 1:3
> expand.grid(x = 1:3, y = 1:3)
  x y
1 1 1
2 2 1
3 3 1
4 1 2
5 2 2
6 3 2
7 1 3
8 2 3
9 3 3

使用expand.grid给了我所有的组合。但是，我只想要成对比较，也就是说，我不想要 1 对 1、2 对、2 或 3 对 3 的比较。此外，我只想保留唯一的对，即我想保留1 对 2（而不是 2 对 1）。

总之，对于上述x和y，我想要以下 3 对组合：

  x y
1 1 2
2 1 3
3 2 3

同样，对于x = …

我目前正在尝试各种方法： 1. 相关性。2. 相互信息。3. 距离相关性以找出 X 中的变量与 y 中的因变量之间的关系强度。关联是最快和最简单的（一个样本 1 小时到 300 万条记录和 560 个变量）。相互信息计算大约需要 16 个小时。我也在研究距离相关性，因为它有一个有趣的特性：Xi 和 Y 之间的距离相关性为零，当且仅当它们是独立的。但是我在用 Python 进行计算时遇到了一个问题。

以下是我的数据：

X

prop_tenure prop_12m    prop_6m prop_3m 
0.04        0.04        0.06    0.08
0           0           0       0
0           0           0       0
0.06        0.06        0.1     0
0.38        0.38        0.25    0
0.61        0.61        0.66    0.61
0.01        0.01        0.02    0.02
0.1         0.1         0.12    0.16
0.04        0.04        0.04    0.09
0.22        0.22        0.22    0.22
0.72        0.72        0.73    0.72
0.39        0.39        0.45    0.64

**y**
status
0
0 …

例如,假设我们在Julia中有一个6元素数组Int64[1,1,2,3,3,4].如果我们想要按元素比较两个数组,我们知道我们可以使用".=="; 但我的目标是在上面的数组中进行所有成对比较:如果每对中的元素(i,j)相等,我将其设置为1(或者为真),但如果它们不同,我将其设置为0所有成对比较都存储在6x6矩阵中.没有循环可以在朱莉娅那样做吗？谢谢.

我有 DNA 序列数据。例如，

X="ACGGGT"
Y="ACGGT"

我想知道对齐分数，因此我使用了biopythonpairwise2函数。例如，

from Bio import pairwise2
from Bio.pairwise2 import format_alignment

alignments = pairwise2.align.globalxx(X, Y)
for a in alignments:
    print(format_alignment(*a))

这成功地显示了 DNA 比对，但我只需要如下的分数。有没有办法只显示分数？

我使用了biopython，但如果有更好的方法，我们将不胜感激。