doi
/
NHP-pRF
forked from ChrisKlink/NHP-pRF


			
							12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091
							<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"
                "http://www.w3.org/TR/REC-html40/loose.dtd">
<html>
<head>
  <title>Description of lfex1</title>
  <meta name="keywords" content="lfex1">
  <meta name="description" content="Local Estimation a spike firing rate (in spikes per unit time).">
  <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
  <meta name="generator" content="m2html &copy; 2005 Guillaume Flandin">
  <meta name="robots" content="index, follow">
  <link type="text/css" rel="stylesheet" href="../../../m2html.css">
  <script type="text/javascript">
    if (top.frames.length == 0) { top.location = "../../../index.html"; };
  </script>
</head>
<body>
<a name="_top"></a>
<!-- ../../menu.html chronux_2_10 --><!-- # locfit --><!-- menu.html Neuro -->
<h1>lfex1
</h1>

<h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
<div class="box"><strong>Local Estimation a spike firing rate (in spikes per unit time).</strong></div>

<h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
<div class="box"><strong>This is a script file. </strong></div>

<h2><a name="_description"></a>DESCRIPTION <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
<div class="fragment"><pre class="comment"> Local Estimation a spike firing rate (in spikes per unit time).

 The estimation procedure approximates the log of the spike firing
 rate by a quadratic polynomial, within a sliding window.

 The first argument to locfit() is a column vector of the spike times.
 'family','rate' specifies that the output of the density estimate should
 be in terms of events per unit time. (instead of the default, statistical
 density estimation).

 'alpha',0.6 specifies the width of the sliding windows, as a fraction of
 the total spikes. 

 xlim gives the limits of the observation interval. Correct specification
 of this is critical to avoid introducing bias at end-points.

 The lfband(fit) line adds confidence bands (based on pointwise 95%
 coverage) to the plot.

 The spike time data is from Hemant Bokil.</pre></div>

<!-- crossreference -->
<h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
This function calls:
<ul style="list-style-image:url(../../../matlabicon.gif)">
<li><a href="../../../chronux_2_10/locfit/m/lfband.html" class="code" title="function lfband(fit,varargin)">lfband</a>	adds confidence bands around the plot of a locfit() fit.</li><li><a href="../../../chronux_2_10/locfit/m/lfplot.html" class="code" title="function lfplot(varargin)">lfplot</a>	Plot (for one or two dimensions) a locfit() fit.</li><li><a href="../../../chronux_2_10/locfit/m/locfit.html" class="code" title="function fit=locfit(varargin)">locfit</a>	Smoothing noisy data using Local Regression and Likelihood.</li></ul>
This function is called by:
<ul style="list-style-image:url(../../../matlabicon.gif)">
</ul>
<!-- crossreference -->


<h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
<div class="fragment"><pre>0001 <span class="comment">% Local Estimation a spike firing rate (in spikes per unit time).</span>
0002 <span class="comment">%</span>
0003 <span class="comment">% The estimation procedure approximates the log of the spike firing</span>
0004 <span class="comment">% rate by a quadratic polynomial, within a sliding window.</span>
0005 <span class="comment">%</span>
0006 <span class="comment">% The first argument to locfit() is a column vector of the spike times.</span>
0007 <span class="comment">% 'family','rate' specifies that the output of the density estimate should</span>
0008 <span class="comment">% be in terms of events per unit time. (instead of the default, statistical</span>
0009 <span class="comment">% density estimation).</span>
0010 <span class="comment">%</span>
0011 <span class="comment">% 'alpha',0.6 specifies the width of the sliding windows, as a fraction of</span>
0012 <span class="comment">% the total spikes.</span>
0013 <span class="comment">%</span>
0014 <span class="comment">% xlim gives the limits of the observation interval. Correct specification</span>
0015 <span class="comment">% of this is critical to avoid introducing bias at end-points.</span>
0016 <span class="comment">%</span>
0017 <span class="comment">% The lfband(fit) line adds confidence bands (based on pointwise 95%</span>
0018 <span class="comment">% coverage) to the plot.</span>
0019 <span class="comment">%</span>
0020 <span class="comment">% The spike time data is from Hemant Bokil.</span>
0021 
0022 load lmem5312.mat;
0023 fit = <a href="../../../chronux_2_10/locfit/m/locfit.html" class="code" title="function fit=locfit(varargin)">locfit</a>(data(6).times{1},<span class="string">'xlim'</span>,[78.9 80.30],<span class="string">'family'</span>,<span class="string">'rate'</span>,<span class="string">'nn'</span>,0.6);
0024 <a href="../../../chronux_2_10/locfit/m/lfplot.html" class="code" title="function lfplot(varargin)">lfplot</a>(fit);
0025 title(<span class="string">'Spike Firing Rate Estimation'</span>);
0026 <a href="../../../chronux_2_10/locfit/m/lfband.html" class="code" title="function lfband(fit,varargin)">lfband</a>(fit);</pre></div>
<hr><address>Generated on Fri 12-Aug-2011 11:36:15 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/" target="_parent">m2html</a></strong> &copy; 2005</address>
</body>
</html>