User contributions for Brb
Jump to navigation
Jump to search
11 March 2024
- 21:4121:41, 11 March 2024 diff hist +375 N File:Filtered R.png <syntaxhighlight lang='r'> theta <- seq(0, .80, .01) R_BH <- filtered_R(alpha=.10, S2, p2, theta, method="BH") which.max(R_BH) # 60% <---- so theta=0.6 is the optimal filtering threshold # 61 max(R_BH) # [1] 380 plot(theta, R_BH, type="l", xlab=expression(theta), ylab="Rejections", main="BH cutoff = 0.1") abline(v=.6, lty=2) </syntaxhighlight> current
- 21:3521:35, 11 March 2024 diff hist +3,823 N File:Filtered p.png Note: # x-axis "p cutoff" should be "BH cutoff" or "FDR cutoff". # Each curve represents theta (filtering threshold). For example, theta=.1 means 10% of genes are filtered out before we do multiple testing (or BH adjustment). # It is seen the larger the theta, the more hypotheses are rejected at the same FDR cutoff. For example, #* if theta=0, 251 hypotheses are rejected at FDR=.1 #* if theta=.5, 355 hypotheses are rejected at FDR=.1. <syntaxhighlight lang='r'> BiocManager::install("ALL")... current
- 17:2017:20, 11 March 2024 diff hist +366 Statistics →Correlated data
- 12:3712:37, 11 March 2024 diff hist +555 Genome →Low read count and filtering
- 11:1711:17, 11 March 2024 diff hist +1,447 Genome →DNA methylation, Epigenetics
10 March 2024
- 20:1020:10, 10 March 2024 diff hist +227 Genome →TCGAbiolinks
9 March 2024
- 22:5722:57, 9 March 2024 diff hist +55 NAS →Zimablade
- 22:5322:53, 9 March 2024 diff hist +834 NAS →ODroid xu4
- 19:5719:57, 9 March 2024 diff hist +23 Docker →CasaOS
- 19:5219:52, 9 March 2024 diff hist +137 Docker →CasaOS
8 March 2024
- 16:5416:54, 8 March 2024 diff hist +54 Regression →Example 2: data with outliers
- 16:5116:51, 8 March 2024 diff hist +366 Regression →Example 1: data with outliers
- 16:4916:49, 8 March 2024 diff hist +805 N File:DataOutliers2.png {{Pre}} puree <- read.csv("https://gist.githubusercontent.com/arraytools/e851ed88c7456779557fbf3ed67b157a/raw/9971c61fea1db99acbd9de17ea82679ba9811358/dataOutliers2.csv", header=F) plot(puree[,1], puree[, 2], xlab="X", ylab="Y") abline(lm(V2 ~ V1, data = puree)) # robust regression require(MASS) summary(rlm(V2 ~ V1, data = puree)) abline(rr.huber <- rlm(V2 ~ V1, data = puree), col = "blue") # quantile regression library(quantreg) abline(rq(V2 ~ V1, data=puree, tau = 0.5), col = "red") # theil... current
- 15:1915:19, 8 March 2024 diff hist +155 R packages →Posit Package Manager/RStudio Package Manager
- 12:3912:39, 8 March 2024 diff hist +7 Genome →Low read count and filtering
- 12:1312:13, 8 March 2024 diff hist +1,603 Genome →Low read count and filtering
- 11:1711:17, 8 March 2024 diff hist +1,043 Genome →Low read count and filtering
- 09:1409:14, 8 March 2024 diff hist +95 Online tools →Extract the values of data in a plot
- 08:5908:59, 8 March 2024 diff hist +59 Online tools →Extract the values of data in a plot
- 08:4208:42, 8 March 2024 diff hist +5 Jetson →General information
- 08:3408:34, 8 March 2024 diff hist +242 Video →Video Trimmer
7 March 2024
- 18:2218:22, 7 March 2024 diff hist +130 Recipes →鍋子噴火
- 18:1818:18, 7 March 2024 diff hist +202 R →factor(x , levels = ...) vs levels(x) <-
- 18:1518:15, 7 March 2024 diff hist +255 R →factor(x , levels = ...) vs levels(x) <-
- 18:0418:04, 7 March 2024 diff hist +198 R →stats::relevel()
- 17:5017:50, 7 March 2024 diff hist +321 Docker Applications →Google photo
- 15:4715:47, 7 March 2024 diff hist +143 Ubuntu →Hardware/Device manager
- 15:4115:41, 7 March 2024 diff hist +122 Reproducible →Share your code and data
- 15:3815:38, 7 March 2024 diff hist +113 KVM →Windows VM current
- 15:3615:36, 7 March 2024 diff hist +200 Ubuntu →How to Connect to WiFi using netplan
- 15:3315:33, 7 March 2024 diff hist +502 Regression →Robust regression
- 15:2715:27, 7 March 2024 diff hist +899 N File:DataOutliers.png {{Pre}} puree <- read.csv("https://gist.githubusercontent.com/arraytools/47d3a46ae1f9a9cd47db350ae2bd2338/raw/b5cccc8e566ff3bef81b1b371e8bfa174c98ef38/dataOutliers.csv", header = FALSE) plot(puree[,1], puree[, 2], xlim=c(0,1), ylim=c(0,1), xlab="X", ylab="Y") abline(0,1, lty=2) abline(lm(V2 ~ V1, data = puree)) # robust regression require(MASS) summary(rlm(V2 ~ V1, data = puree)) abline(rr.huber <- rlm(V2 ~ V1, data = puree), col = "blue") # almost overlapped with lm() # quantile regressio... current
- 13:5213:52, 7 March 2024 diff hist +68 Life →Draw 繪畫
- 13:2313:23, 7 March 2024 diff hist +272 Mac →Screenshot utility
- 13:1713:17, 7 March 2024 diff hist +710 Mac →Quicktime/QT player
- 12:2012:20, 7 March 2024 diff hist +171 Regression →Robust regression
- 09:4709:47, 7 March 2024 diff hist +67 Prediction →rpart and decision tree current
6 March 2024
- 22:0722:07, 6 March 2024 diff hist +663 Motioneye →Manually record a video
- 17:0517:05, 6 March 2024 diff hist +454 Motioneye →Motion trigger, Pi zero W
- 13:4913:49, 6 March 2024 diff hist +1,036 Video →Handbrake to reduce file size
- 10:2510:25, 6 March 2024 diff hist +245 Kodi →TV
5 March 2024
- 22:1522:15, 5 March 2024 diff hist +33 Linux →at command: Schedule a task
- 22:1522:15, 5 March 2024 diff hist +301 Vlc →Timer for recording
- 22:0622:06, 5 March 2024 diff hist +266 Linux →at command: Schedule a task
- 21:4521:45, 5 March 2024 diff hist +406 Vlc →Timer for recording
- 19:0319:03, 5 March 2024 diff hist +817 Vlc →Use VLC to record internet radio (suitable if the stream is continuous)
- 15:4415:44, 5 March 2024 diff hist +12 Prediction →rpart and decision tree
- 15:3915:39, 5 March 2024 diff hist +661 Prediction →rpart and decision tree
- 12:4012:40, 5 March 2024 diff hist +236 Prediction →rpart and decision tree
- 12:3012:30, 5 March 2024 diff hist +483 Prediction →rpart and decision tree