問題タブ [neuroscience]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
r - Rで、それらのヒストグラムにスパイク時間の指数分布をプロットする方法は?
したがって、私の質問は、前回の質問の後の展開に続きます。スパイクトレインのラスタープロットとしてスパイク時間を取得しようとしています。発火率を 100 に設定し、20 回の試行でスパイク トレインを取得しました。そのコードは次のとおりです。
これにより、次の結果が得られます。
これがすべて完了した後の次のタスクは、Inter-Spike 間隔のベクトルを取得し、それらのヒストグラムを取得することでした。ISI の分布は指数分布に従うため、同じデータで ISI の指数分布をプロットすると、ヒストグラムの高さによって作成される曲線と一致します。したがって、スパイク間のタイミングを最初に取得するために、次を使用しました。
次に、スパイク間の間隔とそのヒストグラムのベクトルを取得するために、次のコマンド ラインを使用しました。
そして、それは私にこのプロットを与えました：
今、私が欲しいのは、スパイク間の間隔の指数分布の性質を正当化するヒストグラム内の指数分布をプロットすることです。使用するパラメーターと使用するレートについて混乱しています。誰かが Interspike インターバル プロットで作業したことがある場合は、助けてください。また、データが不完全なようで申し訳ありませんが、何か不足している場合はお知らせください。
r - コードを遅くすることなく、Rのリストの各ベクトル(長さの異なる)に最初と最後の要素を追加する方法は?
私はRとStackOverFlowが初めてです。ですから、何か足りないものがある場合は、お知らせください。神経の振る舞いをよりよく理解するために、不均一なプロセスをシミュレートすることに取り組んでいます。私のコードは、私が持っている方法で動作します。たとえば、20試行はそれぞれ数秒続き2ます (各試行はスパイクトレインを表します)。次に、リストSpikeTimesには20ベクトルのリストが表示されます。各ベクトルは、その特定の試行でスパイクが発生した場所のタイムスタンプに対応しています。[例えば。SpikeTimes[1]これは0.002250802 0.053934034...1.971574170 2.023090384、最初のスパイク トレインで でスパイクが発生したことを意味0.002250802し0.053934034ます。2制限時間の秒数を超えてタイムスタンプが表示される理由はわかりませんが、後で作業します]。私のコードは次のようになります--
私の問題は、リスト内の各ベクトルについてSpikeTimesです。これはスパイクのタイム スタンプを提供しますが、スパイク列の開始 (つまり0) と終了 (つまり )も含めたいと考えています。2したがって、このリストを追加して、各ベクトルに最初のエントリが として含まれ0、最後のエントリが として含まれるようにし2ます。
私の場合は次のSpikeTimes[1]ように見え0 0.002250802 0.053934034...1.971574170 2、他の人SpikeTimes[i]は同様に見えます。SpikeTimes <- c(0, SpikeTimes)最初に入力しようとし0ましたが、最初の要素としてベクトルを使用して、リストに20個ではなく21個のベクトルしかありませ0んでした（なぜそれが起こったのかがわかります）。コードが遅くならないようにするにはどうすればよいですか? 私はRの初心者であり、インターネットで読んでもこの特定の問題の解決にはなりません。あらゆる種類の入力をいただければ幸いです。