kernel-density

    3熱度

    1回答

    我有一個數據框中的蛋白質 - 蛋白質相互作用數據,標題爲s1m。每個DB和AD對做出的互動,我可以繪製它還有: > head(s1m) DB_num AD_num [1,] 2 8153 [2,] 7 3553 [3,] 8 4812 [4,] 13 7838 [5,] 24 3315 [6,] 24 6012 圖中的數據的樣子: 然後我用這個網站的代碼,我發

    0熱度

    3回答

    我對R非常非常新,所以請原諒我的問題的基本性質。總之,我做了很多谷歌搜索試圖回答這個問題,但是我發現即使是基本的指南,以及論壇上的簡單討論,都會比我擁有更多的先前知識,特別是在概述所有編碼術語是什麼,改變它們意味着什麼。 總之,我有一個標籤格式化表格,其中包含三列數據,我希望在單個圖表上繪製密度圖。我希望線條可以是不同的圖案(虛線,虛線等等,因爲任何東西都可以很容易地將它們區分開來,因爲我的主管是

    5熱度

    2回答

    我有一個x,y點分佈,我得到KDE到scipy.stats.gaussian_kde。這是我的代碼,以及如何輸出如下(在x,y數據可以從here獲得): import numpy as np from scipy import stats # Obtain data from file. data = np.loadtxt('data.dat', unpack=True) m1, m2

    0熱度

    1回答

    我正在製作一個像這樣的密度圖,How to create a density plot in matplotlib?即時將axvlines設置爲圖上的少量點,我的問題是我需要知道什麼是最高峯的精確x值。 我可以用一個循環找到它:這款BB有峯值的x值,但要做到這一點循環的時間 density = gaussian_kde(data) aa = 0 bb = 0 for

    8熱度

    3回答

    我需要儘可能精確地找到核密度估計的峯值(連續隨機變量的模態值)。我能找到的近似值: x<-rlnorm(100) d<-density(x) plot(d) i<-which.max(d$y) d$y[i] d$x[i] 但是,當計算d$y精確的功能是已知的。我怎樣才能找到模式的確切值?

    1熱度

    1回答

    我具有由大小1×5的向量的數據,每個代表一個派克:[x,y,r,g,b],x和y是位置:0 <= x <= M,0 <= y <= N。 r,g,b是像素的顏色:0 <= r,g,b <= 255。 我想使用多變量Epanechnikov核估計密度估計。我讀了基本上有兩種方法: 乘法方法 - 計算每個維度的內核,然後將它們相乘。 計算向量的範數並計算該值的內核。 這兩種方法究竟是如何處理我的數據?

    0熱度

    1回答

    我需要一張不同密度線的圖,每張圖都有另一種顏色。這是一個示例代碼(但要小得多),使用內置的data.fame USArrests。我希望可以使用它? colors <- heat.colors(3) plot(density(USArrests[,2], bw=1, kernel="epanechnikov", na.rm=TRUE),col=colors[1]) lines1E <- fu

    4熱度

    1回答

    我有一個關於kde2d (Kernel density estimator).的問題我正在爲變量的相同空間中的兩組不同的數據計算兩個不同的kde2d。當我用filled.contour2或contour比較時,我發現散點圖中具有較低密度點的集合(其中總點數少於10個)的集合具有較高的輪廓值密度。我期望具有更高點密度的組具有更高的密度輪廓值,但正如我上面所說的那樣並非如此。必須選擇帶寬(h)?我使用

    2熱度

    1回答

    我正在使用beanplot包的beanplot函數,並且我找不到一種方法來使用wd=參數並獲得好的結果。 我想要什麼 顯示beanplots的寬度取決於樣本的大小。 (現在)明白'wd ='的用法以及如何使用它。 到目前爲止,當我嘗試使用WD參數,作爲列表它給了我一個錯誤,作爲一個載體,它給了我一個很奇怪(WD好像乘到密度估計值) 例 library(beanplot) set.seed(200

    1熱度

    2回答

    我有點通過在邊緣的情況下將R密度()函數的行爲不解... 假設我添加越來越多的具有x點= 0進入模擬數據集。我期望的是,密度估計值將在x = 0時快速收斂(我對這意味着......的意圖模糊)。在實踐中,配合肯定會越走越窄,但很慢,見下面的序列圖中: plot(density(c(0,0)), xlim=c(-2,2)) plot(density(c(0,0,0,0)), xlim=c(-2,2