cockroachALData: Spike Trains of several Cockroach Antennal Lobe Neurons...
In STAR: Spike Train Analysis with R
Description
Usage
Format
Details
Source
References
Examples
Description
Four (CAL1S and CAL1V), three (CAL2S and
CAL2C), three (e060517spont and e060517ionon),
three (e060817spont, e060817terpi, e060817citron
and e060817mix), two (e060824spont and
e060824citral) and four (e070528spont and
e070528citronellal) Cockroach (Periplaneta americana) antennal lobe neurons
(putative projection neurons) were recorded simultaneously
and extracellularly during spontaneous activity and odors (vanilin,
citral, citronellal, terpineol, beta-ionon)
responses from six different animals. The data sets contain the sorted
spike trains of the neurons.
Usage
1
2
3
4
5
6
7
8
9
10
11
12
13
14
Format
CAL1S is a named list with 4 components
("neuron 1", "neuron 2", "neuron 3",
"neuron 4"). Each
component contains the spike train (ie, action potentials occurrence
times) of one neuron recorded during 30 s of spontaneous
activity. Times are expressed in seconds.
CAL1V is a named list with 4 components
("neuron 1", "neuron 2", "neuron 3",
"neuron 4").
Each component is a named list
with 20 components: "stim. 1", ..., "stim. 20". Each
sub-list contains the spike train of one neuron during 1 stimulation
(odor puff) with vanillin
(http://en.wikipedia.org/wiki/Vanillin). Each acquisition was 10
s long. The command
to the odor delivery valve was on between sec 4.49 and sec 4.99.
CAL2S is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3"). Each
component contains the spike train (ie, action potentials occurrence
times) of one neuron recorded during 1 mn of spontaneous
activity. Times are expressed in seconds.
CAL2C is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3").
Each component is a named list
with 20 components: "stim. 1", ..., "stim. 20". Each
sub-list contains the spike train of one neuron during 1 stimulation
(odor puff) with citral
(http://en.wikipedia.org/wiki/Citral). Each acquisition was 14 s long. The command
to the odor delivery valve was on between sec 5.87 and sec 6.37.
e060517spont is a named list of with 3 components
("neuron 1", "neuron 2", "neuron 3"). Each
component is a spikeTrain object (ie, action potentials occurrence
times) of one neuron recorded during 61 s of spontaneous
activity. Times are expressed in seconds.
e060517ionon is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3").
Each component is a repeatedTrain object
with 19 spikeTrain objects: "stim. 1", ..., "stim. 19". Each
spikeTrain contains the spike train of one neuron during 1 stimulation
(odor puff) with beta-ionon
(http://commons.wikimedia.org/wiki/Image:Beta-Ionon.svg). Each acquisition was 15 s long. The command
to the odor delivery valve was on between sec 6.07 and sec 6.57.
e060817spont is a named list of with 3 components
("neuron 1", "neuron 2", "neuron 3"). Each
component is a spikeTrain object (ie, action potentials occurrence
times) of one neuron recorded during 60 s of spontaneous
activity. Times are expressed in seconds.
e060817terpi is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3").
Each component is a repeatedTrain object
with 20 spikeTrain objects: "stim. 1", ..., "stim. 20". Each
spikeTrain contains the spike train of one neuron during 1 stimulation
(odor puff) with terpineol
(http://en.wikipedia.org/wiki/Terpineol). Each acquisition was 15 s long. The command
to the odor delivery valve was on between sec 6.03 and sec 6.53.
e060817citron is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3").
Each component is a repeatedTrain object
with 20 spikeTrain objects: "stim. 1", ..., "stim. 20". Each
spikeTrain contains the spike train of one neuron during 1 stimulation
(odor puff) with citronellal
(http://en.wikipedia.org/wiki/Citronellal). Each acquisition was 15 s long. The command
to the odor delivery valve was on between sec 5.99 and sec 6.49.
e060817mix is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3").
Each component is a repeatedTrain object
with 20 spikeTrain objects: "stim. 1", ..., "stim. 20". Each
spikeTrain contains the spike train of one neuron during 1 stimulation
(odor puff) with a mixture of terpinaol and citronellal
(the sum of the two previous stim.). Each acquisition was 15 s long. The command
to the odor delivery valve was on between sec 6.01 and sec 6.51.
e060824spont is a named list of with 2 components
("neuron 1", "neuron 2"). Each
component is a spikeTrain object (ie, action potentials occurrence
times) of one neuron recorded during 59 s of spontaneous
activity. Times are expressed in seconds.
e060824citral is a named list with 2 components
("neuron 1", "neuron 2").
Each component is a named list
with 20 components: "stim. 1", ..., "stim. 20". Each
sub-list contains the spike train of one neuron during 1 stimulation
(odor puff) with citral
(http://en.wikipedia.org/wiki/Citral). Each acquisition was 15 s long. The command
to the odor delivery valve was on between sec 6.01 and sec 6.51.
e070528spont is a named list of with 4 components
("neuron 1", "neuron 2", "neuron 3",
"neuron 4"). Each
component is a spikeTrain object (ie, action potentials occurrence
times) of one neuron recorded during 60 s of spontaneous
activity. Times are expressed in seconds.
e070528citronellal is a named list with 4 components
("neuron 1", "neuron 2", "neuron 3",
"neuron 4").
Each component is a repeatedTrain object
with 15 spikeTrain objects: "stim. 1", ..., "stim. 15". Each
spikeTrain contains the spike train of one neuron during 1 stimulation
(odor puff) with citronellal
(http://en.wikipedia.org/wiki/Citronellal). Each acquisition was 13 s long. The command
to the odor delivery valve was on between sec 6.14 and sec 6.64.
Details
Every repeatedTrain object of these data sets has an attribute named
stimTimeCourse containing the openng and closing times of the
odor delivery valve.
The data were recorded from neighboring sites on a NeuroNexus
(http://neuronexustech.com/) silicon probe. Sorting was done
with SpikeOMatic with superposition resolution which can AND
DOES lead to artifcats on cross-correlograms.
Source
Recording and spike sorting performed by Antoine Chaffiol
antoine.chaffiol@univ-paris5.fr at the Cerebral Physiology
Lab, CNRS UMR 8118:
http://www.biomedicale.univ-paris5.fr/physcerv/physiologie_cerebrale.htm.
References
http://www.biomedicale.univ-paris5.fr/physcerv/C_Pouzat/Doc/ChaffiolEtAl_FENS2006.pdf
Examples
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
## load CAL1S data
data(CAL1S)
## convert the data into spikeTrain objects
CAL1S <- lapply(CAL1S,as.spikeTrain)
## look at the train of the 1sd neuron
CAL1S[["neuron 1"]]
## fit the 6 different renewal models to the 1st neuron spike train
compModels(CAL1S[["neuron 1"]])
## look at the ISI distribution with the fitted invgauss dist for
## this 1st neuron
isiHistFit(CAL1S[["neuron 1"]],model="invgauss")
## load CAL1V data
data(CAL1V)
## convert them to repeatedTrain objects
CAL1V <- lapply(CAL1V, as.repeatedTrain)
## look at the raster of the 1st neuron
CAL1V[["neuron 1"]]
## load e070528spont data
data(e070528spont)
## look at the spike train of the 1st neuron
e070528spont[["neuron 1"]]
## load e070528citronellal data
data(e070528citronellal)
## Get the stimulus time course
attr(e070528citronellal[["neuron 1"]],"stimTimeCourse")
## look at the raster of the 1st neuron
plot(e070528citronellal[["neuron 1"]],stim=c(6.14,6.64))
## Not run:
## A "detailed" analysis of e060817 were 2 odors as well as there mixtures
## were used.
## Load the terpineol, citronellal and mixture response data
data(e060817terpi)
data(e060817citron)
data(e060817mix)
## get smooth psths with gsspsth0
e060817terpiN1PSTH <- gsspsth0(e060817terpi[["neuron 1"]])
e060817terpiN2PSTH <- gsspsth0(e060817terpi[["neuron 2"]])
e060817terpiN3PSTH <- gsspsth0(e060817terpi[["neuron 3"]])
e060817citronN1PSTH <- gsspsth0(e060817citron[["neuron 1"]])
e060817citronN2PSTH <- gsspsth0(e060817citron[["neuron 2"]])
e060817citronN3PSTH <- gsspsth0(e060817citron[["neuron 3"]])
e060817mixN1PSTH <- gsspsth0(e060817mix[["neuron 1"]])
e060817mixN2PSTH <- gsspsth0(e060817mix[["neuron 2"]])
e060817mixN3PSTH <- gsspsth0(e060817mix[["neuron 3"]])
## look at them
## Neuron 1
plot(e060817terpiN1PSTH,stimTimeCourse=attr(e060817terpi[["neuron 1"]],"stimTimeCourse"),colCI=2)
plot(e060817citronN1PSTH,stimTimeCourse=attr(e060817citron[["neuron 1"]],"stimTimeCourse"),colCI=2)
plot(e060817mixN1PSTH,stimTimeCourse=attr(e060817mix[["neuron 1"]],"stimTimeCourse"),colCI=2)
## Neuron 2
plot(e060817terpiN2PSTH,stimTimeCourse=attr(e060817terpi[["neuron 2"]],"stimTimeCourse"),colCI=2)
plot(e060817citronN2PSTH,stimTimeCourse=attr(e060817citron[["neuron 2"]],"stimTimeCourse"),colCI=2)
plot(e060817mixN2PSTH,stimTimeCourse=attr(e060817mix[["neuron 2"]],"stimTimeCourse"),colCI=2)
## Neuron 3
plot(e060817terpiN3PSTH,stimTimeCourse=attr(e060817terpi[["neuron 3"]],"stimTimeCourse"),colCI=2)
plot(e060817citronN3PSTH,stimTimeCourse=attr(e060817citron[["neuron 3"]],"stimTimeCourse"),colCI=2)
plot(e060817mixN3PSTH,stimTimeCourse=attr(e060817mix[["neuron 3"]],"stimTimeCourse"),colCI=2)
## Make now fancier plots with superposed psths ####
## Take into account the fact that the stimuli onsets are not identical
## Neuron 1
plot(e060817mixN1PSTH$mids-0.02,e060817mixN1PSTH$ciUp,type="n",ylim=c(0,max(e060817mixN1PSTH$ciUp)),xlim=c(5,14),xlab="Time (s)",ylab="Firing rate (Hz)",main="Neuron 1 e060817")
rect(5.99,0,6.49,max(e060817mixN1PSTH$ciUp),col="grey80",border=NA)
abline(h=0)
polygon(c(e060817mixN1PSTH$mids-0.02,rev(e060817mixN1PSTH$mids-0.02)),c(e060817mixN1PSTH$ciLow,rev(e060817mixN1PSTH$ciUp)),col=rgb(1,0,1,0.5),border=NA)
polygon(c(e060817citronN1PSTH$mids,rev(e060817citronN1PSTH$mids)),c(e060817citronN1PSTH$ciLow,rev(e060817citronN1PSTH$ciUp)),col=rgb(1,0,0,0.5),border=NA)
polygon(c(e060817terpiN1PSTH$mids-0.04,rev(e060817terpiN1PSTH$mids-0.04)),c(e060817terpiN1PSTH$ciLow,rev(e060817terpiN1PSTH$ciUp)),col=rgb(0,0,1,0.5),border=NA)
lines(e060817terpiN1PSTH$mids-0.04,e060817terpiN1PSTH$freq,col=rgb(0,0,1),lwd=2)
lines(e060817citronN1PSTH$mids,e060817citronN1PSTH$freq,col=rgb(1,0,0),lwd=2)
lines(e060817mixN1PSTH$mids-0.02,e060817mixN1PSTH$freq,col=rgb(0,0,0),lwd=2)
legend(8,0.9*max(e060817mixN1PSTH$ciUp),c("Terpineol","Citronellal","Mixture"),col=c(4,2,1),lwd=2)
## Neuron 2
plot(e060817mixN2PSTH$mids-0.02,e060817mixN2PSTH$ciUp,type="n",ylim=c(0,max(e060817mixN2PSTH$ciUp)),xlim=c(5,14),xlab="Time (s)",ylab="Firing rate (Hz)",main="Neuron 2 e060817")
rect(5.99,0,6.49,max(e060817mixN2PSTH$ciUp),col="grey80",border=NA)
abline(h=0)
polygon(c(e060817mixN2PSTH$mids-0.02,rev(e060817mixN2PSTH$mids-0.02)),c(e060817mixN2PSTH$ciLow,rev(e060817mixN2PSTH$ciUp)),col=rgb(1,0,1,0.5),border=NA)
polygon(c(e060817citronN2PSTH$mids,rev(e060817citronN2PSTH$mids)),c(e060817citronN2PSTH$ciLow,rev(e060817citronN2PSTH$ciUp)),col=rgb(1,0,0,0.5),border=NA)
polygon(c(e060817terpiN2PSTH$mids-0.04,rev(e060817terpiN2PSTH$mids-0.04)),c(e060817terpiN2PSTH$ciLow,rev(e060817terpiN2PSTH$ciUp)),col=rgb(0,0,1,0.5),border=NA)
lines(e060817terpiN2PSTH$mids-0.04,e060817terpiN2PSTH$freq,col=rgb(0,0,1),lwd=2)
lines(e060817citronN2PSTH$mids,e060817citronN2PSTH$freq,col=rgb(1,0,0),lwd=2)
lines(e060817mixN2PSTH$mids-0.02,e060817mixN2PSTH$freq,col=rgb(0,0,0),lwd=2)
legend(8,0.9*max(e060817mixN2PSTH$ciUp),c("Terpineol","Citronellal","Mixture"),col=c(4,2,1),lwd=2)
## Neuron 3
plot(e060817mixN3PSTH$mids-0.02,e060817mixN3PSTH$ciUp,type="n",ylim=c(0,max(e060817mixN3PSTH$ciUp)),xlim=c(5,14),xlab="Time (s)",ylab="Firing rate (Hz)",main="Neuron 3 e060817")
rect(5.99,0,6.49,max(e060817mixN3PSTH$ciUp),col="grey80",border=NA)
abline(h=0)
polygon(c(e060817mixN3PSTH$mids-0.02,rev(e060817mixN3PSTH$mids-0.02)),c(e060817mixN3PSTH$ciLow,rev(e060817mixN3PSTH$ciUp)),col=rgb(1,0,1,0.5),border=NA)
polygon(c(e060817citronN3PSTH$mids,rev(e060817citronN3PSTH$mids)),c(e060817citronN3PSTH$ciLow,rev(e060817citronN3PSTH$ciUp)),col=rgb(1,0,0,0.5),border=NA)
polygon(c(e060817terpiN3PSTH$mids-0.04,rev(e060817terpiN3PSTH$mids-0.04)),c(e060817terpiN3PSTH$ciLow,rev(e060817terpiN3PSTH$ciUp)),col=rgb(0,0,1,0.5),border=NA)
lines(e060817terpiN3PSTH$mids-0.04,e060817terpiN3PSTH$freq,col=rgb(0,0,1),lwd=2)
lines(e060817citronN3PSTH$mids,e060817citronN3PSTH$freq,col=rgb(1,0,0),lwd=2)
lines(e060817mixN3PSTH$mids-0.02,e060817mixN3PSTH$freq,col=rgb(0,0,0),lwd=2)
legend(8,0.9*max(e060817mixN3PSTH$ciUp),c("Terpineol","Citronellal","Mixture"),col=c(4,2,1),lwd=2)
## End(Not run)
STAR documentation built on May 2, 2019, 11:44 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Description Usage Format Details Source References Examples
Description
Four (CAL1S and CAL1V), three (CAL2S and
CAL2C), three (e060517spont and e060517ionon),
three (e060817spont, e060817terpi, e060817citron
and e060817mix), two (e060824spont and
e060824citral) and four (e070528spont and
e070528citronellal) Cockroach (Periplaneta americana) antennal lobe neurons
(putative projection neurons) were recorded simultaneously
and extracellularly during spontaneous activity and odors (vanilin,
citral, citronellal, terpineol, beta-ionon)
responses from six different animals. The data sets contain the sorted
spike trains of the neurons.
Usage
1 2 3 4 5 6 7 8 9 10 11 12 13 14 |
Format
CAL1S is a named list with 4 components
("neuron 1", "neuron 2", "neuron 3",
"neuron 4"). Each
component contains the spike train (ie, action potentials occurrence
times) of one neuron recorded during 30 s of spontaneous
activity. Times are expressed in seconds.
CAL1V is a named list with 4 components
("neuron 1", "neuron 2", "neuron 3",
"neuron 4").
Each component is a named list
with 20 components: "stim. 1", ..., "stim. 20". Each
sub-list contains the spike train of one neuron during 1 stimulation
(odor puff) with vanillin
(http://en.wikipedia.org/wiki/Vanillin). Each acquisition was 10
s long. The command
to the odor delivery valve was on between sec 4.49 and sec 4.99.
CAL2S is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3"). Each
component contains the spike train (ie, action potentials occurrence
times) of one neuron recorded during 1 mn of spontaneous
activity. Times are expressed in seconds.
CAL2C is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3").
Each component is a named list
with 20 components: "stim. 1", ..., "stim. 20". Each
sub-list contains the spike train of one neuron during 1 stimulation
(odor puff) with citral
(http://en.wikipedia.org/wiki/Citral). Each acquisition was 14 s long. The command
to the odor delivery valve was on between sec 5.87 and sec 6.37.
e060517spont is a named list of with 3 components
("neuron 1", "neuron 2", "neuron 3"). Each
component is a spikeTrain object (ie, action potentials occurrence
times) of one neuron recorded during 61 s of spontaneous
activity. Times are expressed in seconds.
e060517ionon is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3").
Each component is a repeatedTrain object
with 19 spikeTrain objects: "stim. 1", ..., "stim. 19". Each
spikeTrain contains the spike train of one neuron during 1 stimulation
(odor puff) with beta-ionon
(http://commons.wikimedia.org/wiki/Image:Beta-Ionon.svg). Each acquisition was 15 s long. The command
to the odor delivery valve was on between sec 6.07 and sec 6.57.
e060817spont is a named list of with 3 components
("neuron 1", "neuron 2", "neuron 3"). Each
component is a spikeTrain object (ie, action potentials occurrence
times) of one neuron recorded during 60 s of spontaneous
activity. Times are expressed in seconds.
e060817terpi is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3").
Each component is a repeatedTrain object
with 20 spikeTrain objects: "stim. 1", ..., "stim. 20". Each
spikeTrain contains the spike train of one neuron during 1 stimulation
(odor puff) with terpineol
(http://en.wikipedia.org/wiki/Terpineol). Each acquisition was 15 s long. The command
to the odor delivery valve was on between sec 6.03 and sec 6.53.
e060817citron is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3").
Each component is a repeatedTrain object
with 20 spikeTrain objects: "stim. 1", ..., "stim. 20". Each
spikeTrain contains the spike train of one neuron during 1 stimulation
(odor puff) with citronellal
(http://en.wikipedia.org/wiki/Citronellal). Each acquisition was 15 s long. The command
to the odor delivery valve was on between sec 5.99 and sec 6.49.
e060817mix is a named list with 3 components
("neuron 1", "neuron 2", "neuron 3").
Each component is a repeatedTrain object
with 20 spikeTrain objects: "stim. 1", ..., "stim. 20". Each
spikeTrain contains the spike train of one neuron during 1 stimulation
(odor puff) with a mixture of terpinaol and citronellal
(the sum of the two previous stim.). Each acquisition was 15 s long. The command
to the odor delivery valve was on between sec 6.01 and sec 6.51.
e060824spont is a named list of with 2 components
("neuron 1", "neuron 2"). Each
component is a spikeTrain object (ie, action potentials occurrence
times) of one neuron recorded during 59 s of spontaneous
activity. Times are expressed in seconds.
e060824citral is a named list with 2 components
("neuron 1", "neuron 2").
Each component is a named list
with 20 components: "stim. 1", ..., "stim. 20". Each
sub-list contains the spike train of one neuron during 1 stimulation
(odor puff) with citral
(http://en.wikipedia.org/wiki/Citral). Each acquisition was 15 s long. The command
to the odor delivery valve was on between sec 6.01 and sec 6.51.
e070528spont is a named list of with 4 components
("neuron 1", "neuron 2", "neuron 3",
"neuron 4"). Each
component is a spikeTrain object (ie, action potentials occurrence
times) of one neuron recorded during 60 s of spontaneous
activity. Times are expressed in seconds.
e070528citronellal is a named list with 4 components
("neuron 1", "neuron 2", "neuron 3",
"neuron 4").
Each component is a repeatedTrain object
with 15 spikeTrain objects: "stim. 1", ..., "stim. 15". Each
spikeTrain contains the spike train of one neuron during 1 stimulation
(odor puff) with citronellal
(http://en.wikipedia.org/wiki/Citronellal). Each acquisition was 13 s long. The command
to the odor delivery valve was on between sec 6.14 and sec 6.64.
Details
Every repeatedTrain object of these data sets has an attribute named
stimTimeCourse containing the openng and closing times of the
odor delivery valve.
The data were recorded from neighboring sites on a NeuroNexus
(http://neuronexustech.com/) silicon probe. Sorting was done
with SpikeOMatic with superposition resolution which can AND
DOES lead to artifcats on cross-correlograms.
Source
Recording and spike sorting performed by Antoine Chaffiol antoine.chaffiol@univ-paris5.fr at the Cerebral Physiology Lab, CNRS UMR 8118: http://www.biomedicale.univ-paris5.fr/physcerv/physiologie_cerebrale.htm.
References
http://www.biomedicale.univ-paris5.fr/physcerv/C_Pouzat/Doc/ChaffiolEtAl_FENS2006.pdf
Examples
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 | ## load CAL1S data
data(CAL1S)
## convert the data into spikeTrain objects
CAL1S <- lapply(CAL1S,as.spikeTrain)
## look at the train of the 1sd neuron
CAL1S[["neuron 1"]]
## fit the 6 different renewal models to the 1st neuron spike train
compModels(CAL1S[["neuron 1"]])
## look at the ISI distribution with the fitted invgauss dist for
## this 1st neuron
isiHistFit(CAL1S[["neuron 1"]],model="invgauss")
## load CAL1V data
data(CAL1V)
## convert them to repeatedTrain objects
CAL1V <- lapply(CAL1V, as.repeatedTrain)
## look at the raster of the 1st neuron
CAL1V[["neuron 1"]]
## load e070528spont data
data(e070528spont)
## look at the spike train of the 1st neuron
e070528spont[["neuron 1"]]
## load e070528citronellal data
data(e070528citronellal)
## Get the stimulus time course
attr(e070528citronellal[["neuron 1"]],"stimTimeCourse")
## look at the raster of the 1st neuron
plot(e070528citronellal[["neuron 1"]],stim=c(6.14,6.64))
## Not run:
## A "detailed" analysis of e060817 were 2 odors as well as there mixtures
## were used.
## Load the terpineol, citronellal and mixture response data
data(e060817terpi)
data(e060817citron)
data(e060817mix)
## get smooth psths with gsspsth0
e060817terpiN1PSTH <- gsspsth0(e060817terpi[["neuron 1"]])
e060817terpiN2PSTH <- gsspsth0(e060817terpi[["neuron 2"]])
e060817terpiN3PSTH <- gsspsth0(e060817terpi[["neuron 3"]])
e060817citronN1PSTH <- gsspsth0(e060817citron[["neuron 1"]])
e060817citronN2PSTH <- gsspsth0(e060817citron[["neuron 2"]])
e060817citronN3PSTH <- gsspsth0(e060817citron[["neuron 3"]])
e060817mixN1PSTH <- gsspsth0(e060817mix[["neuron 1"]])
e060817mixN2PSTH <- gsspsth0(e060817mix[["neuron 2"]])
e060817mixN3PSTH <- gsspsth0(e060817mix[["neuron 3"]])
## look at them
## Neuron 1
plot(e060817terpiN1PSTH,stimTimeCourse=attr(e060817terpi[["neuron 1"]],"stimTimeCourse"),colCI=2)
plot(e060817citronN1PSTH,stimTimeCourse=attr(e060817citron[["neuron 1"]],"stimTimeCourse"),colCI=2)
plot(e060817mixN1PSTH,stimTimeCourse=attr(e060817mix[["neuron 1"]],"stimTimeCourse"),colCI=2)
## Neuron 2
plot(e060817terpiN2PSTH,stimTimeCourse=attr(e060817terpi[["neuron 2"]],"stimTimeCourse"),colCI=2)
plot(e060817citronN2PSTH,stimTimeCourse=attr(e060817citron[["neuron 2"]],"stimTimeCourse"),colCI=2)
plot(e060817mixN2PSTH,stimTimeCourse=attr(e060817mix[["neuron 2"]],"stimTimeCourse"),colCI=2)
## Neuron 3
plot(e060817terpiN3PSTH,stimTimeCourse=attr(e060817terpi[["neuron 3"]],"stimTimeCourse"),colCI=2)
plot(e060817citronN3PSTH,stimTimeCourse=attr(e060817citron[["neuron 3"]],"stimTimeCourse"),colCI=2)
plot(e060817mixN3PSTH,stimTimeCourse=attr(e060817mix[["neuron 3"]],"stimTimeCourse"),colCI=2)
## Make now fancier plots with superposed psths ####
## Take into account the fact that the stimuli onsets are not identical
## Neuron 1
plot(e060817mixN1PSTH$mids-0.02,e060817mixN1PSTH$ciUp,type="n",ylim=c(0,max(e060817mixN1PSTH$ciUp)),xlim=c(5,14),xlab="Time (s)",ylab="Firing rate (Hz)",main="Neuron 1 e060817")
rect(5.99,0,6.49,max(e060817mixN1PSTH$ciUp),col="grey80",border=NA)
abline(h=0)
polygon(c(e060817mixN1PSTH$mids-0.02,rev(e060817mixN1PSTH$mids-0.02)),c(e060817mixN1PSTH$ciLow,rev(e060817mixN1PSTH$ciUp)),col=rgb(1,0,1,0.5),border=NA)
polygon(c(e060817citronN1PSTH$mids,rev(e060817citronN1PSTH$mids)),c(e060817citronN1PSTH$ciLow,rev(e060817citronN1PSTH$ciUp)),col=rgb(1,0,0,0.5),border=NA)
polygon(c(e060817terpiN1PSTH$mids-0.04,rev(e060817terpiN1PSTH$mids-0.04)),c(e060817terpiN1PSTH$ciLow,rev(e060817terpiN1PSTH$ciUp)),col=rgb(0,0,1,0.5),border=NA)
lines(e060817terpiN1PSTH$mids-0.04,e060817terpiN1PSTH$freq,col=rgb(0,0,1),lwd=2)
lines(e060817citronN1PSTH$mids,e060817citronN1PSTH$freq,col=rgb(1,0,0),lwd=2)
lines(e060817mixN1PSTH$mids-0.02,e060817mixN1PSTH$freq,col=rgb(0,0,0),lwd=2)
legend(8,0.9*max(e060817mixN1PSTH$ciUp),c("Terpineol","Citronellal","Mixture"),col=c(4,2,1),lwd=2)
## Neuron 2
plot(e060817mixN2PSTH$mids-0.02,e060817mixN2PSTH$ciUp,type="n",ylim=c(0,max(e060817mixN2PSTH$ciUp)),xlim=c(5,14),xlab="Time (s)",ylab="Firing rate (Hz)",main="Neuron 2 e060817")
rect(5.99,0,6.49,max(e060817mixN2PSTH$ciUp),col="grey80",border=NA)
abline(h=0)
polygon(c(e060817mixN2PSTH$mids-0.02,rev(e060817mixN2PSTH$mids-0.02)),c(e060817mixN2PSTH$ciLow,rev(e060817mixN2PSTH$ciUp)),col=rgb(1,0,1,0.5),border=NA)
polygon(c(e060817citronN2PSTH$mids,rev(e060817citronN2PSTH$mids)),c(e060817citronN2PSTH$ciLow,rev(e060817citronN2PSTH$ciUp)),col=rgb(1,0,0,0.5),border=NA)
polygon(c(e060817terpiN2PSTH$mids-0.04,rev(e060817terpiN2PSTH$mids-0.04)),c(e060817terpiN2PSTH$ciLow,rev(e060817terpiN2PSTH$ciUp)),col=rgb(0,0,1,0.5),border=NA)
lines(e060817terpiN2PSTH$mids-0.04,e060817terpiN2PSTH$freq,col=rgb(0,0,1),lwd=2)
lines(e060817citronN2PSTH$mids,e060817citronN2PSTH$freq,col=rgb(1,0,0),lwd=2)
lines(e060817mixN2PSTH$mids-0.02,e060817mixN2PSTH$freq,col=rgb(0,0,0),lwd=2)
legend(8,0.9*max(e060817mixN2PSTH$ciUp),c("Terpineol","Citronellal","Mixture"),col=c(4,2,1),lwd=2)
## Neuron 3
plot(e060817mixN3PSTH$mids-0.02,e060817mixN3PSTH$ciUp,type="n",ylim=c(0,max(e060817mixN3PSTH$ciUp)),xlim=c(5,14),xlab="Time (s)",ylab="Firing rate (Hz)",main="Neuron 3 e060817")
rect(5.99,0,6.49,max(e060817mixN3PSTH$ciUp),col="grey80",border=NA)
abline(h=0)
polygon(c(e060817mixN3PSTH$mids-0.02,rev(e060817mixN3PSTH$mids-0.02)),c(e060817mixN3PSTH$ciLow,rev(e060817mixN3PSTH$ciUp)),col=rgb(1,0,1,0.5),border=NA)
polygon(c(e060817citronN3PSTH$mids,rev(e060817citronN3PSTH$mids)),c(e060817citronN3PSTH$ciLow,rev(e060817citronN3PSTH$ciUp)),col=rgb(1,0,0,0.5),border=NA)
polygon(c(e060817terpiN3PSTH$mids-0.04,rev(e060817terpiN3PSTH$mids-0.04)),c(e060817terpiN3PSTH$ciLow,rev(e060817terpiN3PSTH$ciUp)),col=rgb(0,0,1,0.5),border=NA)
lines(e060817terpiN3PSTH$mids-0.04,e060817terpiN3PSTH$freq,col=rgb(0,0,1),lwd=2)
lines(e060817citronN3PSTH$mids,e060817citronN3PSTH$freq,col=rgb(1,0,0),lwd=2)
lines(e060817mixN3PSTH$mids-0.02,e060817mixN3PSTH$freq,col=rgb(0,0,0),lwd=2)
legend(8,0.9*max(e060817mixN3PSTH$ciUp),c("Terpineol","Citronellal","Mixture"),col=c(4,2,1),lwd=2)
## End(Not run)
|
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.