R
Atomの開発が終了した www.itmedia.co.jp2022年12月15日を過ぎちゃったけどvscodeに移行するのでその忘備録。Atomで使っていたのは主にPython, R, LaTeXで、これらをvscodeでも実行できる環境を用意したい PCにはPython、R、LaTeXをインストール済み これら…
write_adv(adv_url_004); Rを使ってboxplotとbeeswarmを重ねて表示する 一度記事を書いているので今回はその2 mecobalamin.hatenablog.com時系列データを曜日でまとめてbeeswarmでプロットする こんなグラフができる boxplot.png元のデータはこれと同じフォ…
write_adv(adv_url_004); boxplotに元データの点を重ねたい 蜂群図（ほうぐんず）というらしいこちらのサイトを参考にさせてもらった Rで箱ひげ図と蜂群図を合体させる | Takuro FUJITA@NIT, Fukui College ビースウォーム | 値を点としてグラフ上にプロット…
write_adv(adv_url_008); 以前edgeRでリードカウントの正規化をした 同じ結果をDESeq2を使って正規化してみる 正規化の方法が違うらしいDESeq2のマニュアル これは読むべき Analyzing RNA-seq data with DESeq2マニュアルはDESeq2をインストール後に browseV…
検索でこのページが引っかかっているようですが 消しちゃったので別の記事を載せます以前pythonとRでやったことをShell scriptでも書いてみた mecobalamin.hatenablog.com mecobalamin.hatenablog.com書いたスクリプトはこれ #!/bin/bash ppap_array=() ppap…
write_adv(adv_url_005); 21 June 2019追記 文章を少し修正した 追記ここまでfeatrureCountの結果をedgeRを使って正規化するこのtutorialのpdfファイルと以下のリンク先を参考にしたStatistical analysis of RNA-Seq data http://www.nathalievialaneix.eu/d…
Rで時系列データの折れ線グラフを作る グラフを書くことよりも日付の処理が面倒だったのでそのメモ write_adv(adv_url_003); 12 December 2020 追記 データは異なるがPythonでも書いてみた mecobalamin.hatenablog.com 追記ここまで 作りたいのは気温と海水…
Atomが起動しなくなったので 再インストールここからダウンロード Atomインストールしたらそのまま使えた パッケージもインストールされてた 設定ファイルは残っているみたいだ起動しなくなった理由だけど Rのスクリプトをscript pluginで実行したら スペル…
17 May 2019追記 メモした内容で実際に正規化してみた mecobalamin.hatenablog.com mecobalamin.hatenablog.com追記ここまで write_adv(adv_url_003); featureCountsでリードカウントしたあとに 複数のデータを比較する場合、 カウントの正規化が必要らしい…
前回作ったbamファイルを使って発現量を調べる mecobalamin.hatenablog.com使うコマンドはfeatureCountこちらのサイトを参考にした featureCount で bam と gtf からリードカウントを得る | Tips for NGS Data Analysis https://ncrna.jp/blog/item/387-feat…
普段使っているエディターはAtom 殆どいじっていないが、 pythonとRを実行できるようにしてある 必要なpackageは scriptとlanguage-rpackageをインストールしたら windowsのシステム環境変数で Rscriptのディレクトリにパスを通す \R-3.5.0\binwindowsの設定…
これまでグラフはpng形式で出力してきた あとから編集したりなにかに載せるなら ベクター形式で出力しておきたい Scalable Vector Graphics - Wikipediaベクター形式を使うことが殆どなく Illsutratorぐらいしか編集ソフトも知らなかったRではsvg形式で出力…
Rで結果を出力するときwrite.table()を使って CSV形式で出力するのが楽だが、 結果をExcelでくれ！ってなる場合もあるそれにクラスタリングした結果を並べ替えたり 色を付けなおしたりするにはExcel便利なので一度Rで計算して結果をExcelファイルに出力するE…
write_adv(adv_url_003); Rでクラスタリングをしてその結果を ヒートマップと樹状図で表示するには pheatmapのライブラリを使うクラスタリングだけならdist()とhclust()を使えばできるが、 ヒートマップも一緒にとなるとpheatmapかな https://www1.doshisha.…
元データが複数あるなら まとめて処理してくれると便利だ例えばirisのデータが種類ごとに分かれていて かつ地域ごとに保存されているとする chibaとsaitamaの2地域で、3種類に分かれている場合がこれ 180910_Iris-data_chiba_setosa_results.txt 180910_Iris…
irisのデータから散布図とヒストグラムを作る 前回も使ったけどirisのデータセットって？irisの正体を調べた方がいた irisの正体 (R Advent Calendar 2012 6日目) - どんな鳥も "アヤメ属の3種 (Iris setosa，Iris versicolor，Iris virginica) の"、 "sepal…
write_adv(adv_url_003); 何かしらデータを取ったらグラフ作ってエラーバーを加えて 統計的有意差があるかを確認をするのだが データ量が増えるとエクセルだとめんどいRで棒グラフを書いてエラーバーを作れるようにしたirisのデータでグラフを書く 結果はこ…
前回pythonで書いたことを今度はRに移植 p <- c("PEN", "APPLE", "PINEAPPLE") word <- "" while(word != "PEN, PINEAPPLE, APPLE, PEN"){ n <- as.integer(runif(4, min = 1, max = 4)) word <- c(p[n[1]]) for(j in 2:4){ word <- paste(word, p[n[j]], se…
