BACKGROUND
rodents as well as many mammals are characterized by the presence of vibrissae or whiskers located on both sides of the muzzle  <cit> .
when rats acquire sensory information by actively moving their vibrissae, a neural code is manifested at different levels of the sensory system  <cit> .
research in the past  <dig> years has been extensive in the neural code at different stages of the whisker sensory pathway  <cit> .
arabzadeh et al, have demonstrated that a texture code would exist in peripheral afferent response of the vibrissal system and that spike rate would be the coding mechanism that underlies the textures discrimination in the primary somatosensory cortex  <cit> .
in a previous study we have shown that it is possible to characterize the vibrissal sensory information by analyzing the multifiber discharge of a selected vibrissal nerve.
we showed that the afferent activity amplitude  could be related to the physical properties of the surfaces  <cit> .
here, we used a similar protocol to that detailed by albarracín et al  <cit> , but in this case we used different roughness surfaces .
we have analyzed the afferent activity of a single vibrissa innervation  in order to relate the information obtained to the roughness of surfaces.
for this purpose, we have calculated the amplitude of afferent activity by using root mean square values  obtained during protraction and retraction phases.
here, we have also analyzed the effect of slip-resistance levels.
in all cases, the differences among experimental conditions were quantified through the maximum amount of information by using information theory  <cit> .
we found that the discrimination of lower roughness could be achieved at higher slip-resistance levels, whereas the low slip-resistance levels are more appropriate for higher roughness.
we also showed that the experimental conditions are better accurately discriminated in retraction phase.
thus, we have demostrated that discrimination of rough surfaces would require a strategy based on the slip-resistance levels variation.
methods
procedures
five wistar adult male rats , between  <dig> and  <dig> days, were used in our experiments.
they were deeply anesthetized with urethane  and their temperature was maintained at 37° by a servo-controlled heating pad.
surgery consisted of exposing the infraorbital nerve as well as the two branches of the facial nerve  on the right side.
the motor branches were dissected and transected proximally to avoid possible motor influences on the sensorial pathway.
the stimulation electrodes were placed on their distal stumps to produce the contraction of the mystacial muscles.
the deep vibrissal nerve innervating a vibrissal follicle  was identified with a high magnification of a dissecting microscope.
the dissected nerve was also transected proximally and this action allowed eliminating discharges arriving from higher level of the sensorial pathway.
to make sure that the nerve transection did not affect the functionality of the vibrissal nerve during our recording time, we test the falling of the nerve afferent activity throughout the time .
we conclude that the activity start decreasing  <dig> hour after the nerve section, so we never exceeded this space of time in our experiments.
we used a bipolar electrode  to record the multifiber afferent discharge of the vibrissal nerve selected.
the recording electrodes as well as the nerves were immersed in a mineral oil bath during all recording.
all these procedures were carry out in accordance with the recommendations of the guide for the care and use of laboratory animals .
depth of anesthesia was ascertained and controlled during the experiment by the lack of either withdrawal reflex to hindlimb pinching or blink reflex to a gentle stimulation of the cornea.
recording of the vibrissa electrical activity
in this study we have recorded the multifiber activity of the delta vibrissal nerve while the vibrissa was sweeping surfaces of different roughness.
the experimental protocol used in this paper has been previously described in detail by albarracín et al  <cit> .
the procedures are briefly described below.
vibrissa movements were induced by electrical stimulation of facial motor nerve .
square-wave pulses  simulated vibrissal whisking at its natural frequency.
a diagram of the experimental set up is shown in figure 1b.
nerve activity was recorded and digitized at  <dig> khz during a  <dig> ms window following onset of each cycle of whisker movement .
fifty whisker movement cycles were obtained for each surface, and an additional  <dig> cycles were recorded while whisker moved unobstructed in air .
three levels of slip-resistance were presented for each swept surface.
these levels were established by mounting the surface at different distances from the whisker base.
a minimal slip-resistance level was presented by placing the surface at a maximal distance from the whisker base so that the tip just barely contacted the surface throughout the entire movement cycle .
increased slip-resistance levels were presented by moving the surface  <dig> and  <dig> mm closer to the whisker base .
this procedure to increase the slip-resistance is similar to those described by albarracin et al  <cit> .
movements of the delta whisker were recorded simultaneously with nerve activity by using a custom-made photoresistive sensor .
the frequency response of the sensor was maximal in the range 0- <dig> hz, enabling direct identification of the protraction and retraction phases of the movement cycle  <cit> .
rough surfaces
the swept surfaces used in this paper were sandpapers of different grain size: p <dig>  p <dig>  p <dig>  p <dig>  and p <dig> .
we measured the surfaces roughness by using a hommel tester t <dig>  and we used the ra parameter  as a roughness estimation .
ra values obtained were:  <dig> ,  <dig> ,  <dig> ,  <dig>  and  <dig>  μm, respectively.
digital processing and statistics
the rms values and spectral estimation were previously used for the analysis of multifiber recordings  <cit> .
therefore, the theoretical details of both processing techniques can be found in the reference mentioned above.
here, the rms values were obtained and analyzed in:  the whole whisker movement ,  during the protraction and  during the retraction.
statistical analysis was carried out with kruskal wallis anova.
the kruskal wallis test is a nonparametric alternative to the one-way anova.
the test assesses the hypothesis that the different samples in the comparison were drawn from the same distribution or from distributions with the same median.
thus, the interpretation of the kruskal-wallis anova is basically identical to the parametric one-way anova, except that it is based on ranks rather than on mean values.
the dunn's method for multiple comparisons was used as post-hoc test, when anova repeat measurement revealed a significant difference.
these analyses were made with software sigmastat http://www.systat.com/.
information theory
to determine the amount of information in a biological system it is necessary to have at least a pair of stimuli/responses situations.
the stimulus may be a time series or simply belong to a class .
the response depends on the characteristic of the signal that is being examined.
thus, they can be real values  or integer values .
the information conveyed by the neuronal response about the stimulus can be quantified by shannon's mutual information formula  <cit> , abbreviated hereafter as information:
where p is the probability stimulus s presentation, p is the posterior probability of s given the observation of response r, and p is the probability of response r unconditional on the stimulus.
the information determines the maximum amount of knowledge  available to an observer who knows the posterior probabilities p and uses them to read off the signals available in a single observation of a spike train  <cit> .
for a better understanding of the procedure to determine the amount of information, we will consider a specific application in which the stimuli are given by the experimental situations .
the responses will be real values, for example, the rms values of the afferent activities.
these real values will be in the range .
where a and b are the maximum and minimum values of the response.
<dig>  first, the frequency diagrams  determination.
for each experimental situation, the number of occurrences of the response is calculated for each bin.
a graphical way to represent these diagrams is according to a histogram .
the number of bins of the histograms will be calculated directly by the root square of the elements number .
<dig>  determination of the joint probability distribution, p. the following table  <dig> is made.
the response and stimuli are r and s, and their marginal probabilities distributions are p and p, respectively
each element of table  <dig> is the joint probability value, p. it is conceptually defined as follows: p, it is the probability of obtaining a response r = a and that it corresponds to the sweep on p <dig>  this is:
where:
#occurrp <dig> a is the occurrences number obtained for a value of r = a .
#total_occurr is the total number of occurrences in all situations and all possible responses.
for frequency diagrams of figure  <dig>  the joint probability distribution is shown in figure  <dig>
<dig>  determination of the p and p probability distributions.
the probability of obtaining a response r, regardless of whether stimulus s did or not occur, is called the marginal probability, and it can be calculated by the sum of joint probabilities for a given response r. this is:
similarly, we obtain the probability function p.
<dig>  determination of conditional probability distribution p. by definition of conditional probability, is given by:
then, conditional probabilities distributions for each stimulus are obtained from joint and marginal probabilities distributions.
thus, the following probabilities distributions are obtained: p, p,..., p.
the conditional probability, p, conceptually answer the question: given a response r, what is the probability that this response has been produced by the p <dig> stimulus?.
figure  <dig> shows the conditional probability distributions for sweep situations.
<dig>  determination of the amount of information.
after obtaining all probability distributions, it is possible to obtain the mutual information using eq.
<dig>
in this paper, we quantified the information considering the rms values as biological responses of the system.
for the rms values, the information was calculated for each slip-resistance level, as well as for each whisking phases .
RESULTS
the recordings obtained are the average electric activity of about  <dig> myelinated axons which supply innervation to the delta vibrissa.
figure 5b shows the afferent discharge recorded during different sweep situations .
it is possible to distinguish amplitude changes in the signals in relation with the surface roughness.
this feature is more appreciable in figures 5c and 5d.
whisker displacements and evoked afferent activity were simultaneously recorded for all cases.
figure 5a shows the afferent activity evoked by vibrissal sweeps in the air  and its corresponding displacement.
vibrissal movement recordings were acquired with a bandwidth of 0- <dig> hz; thus, only slow movements were monitored.
therefore, displacement recordings were only used to identify the protraction and retraction phases.
we calculated one rms value for each recording obtained during a vibrissa sweeping.
each experimental condition consists of  <dig> rms values .
we first calculated the rms values for the whole-whisk recordings .
the distribution of rms values for each sweep situation is shown in figure  <dig>
the test computes statistic q, the number of rank sums, and shows whether p <  <dig>  or not, for the pair that are being compared.
p is the probability that the null hypothesis may be rejected and, thus, it helps conclude that there are differences between treatments.
diff of ranks is the difference in the rank sum orders that are being compared.
the rank sums are a measurement of the difference between two treatments.
the distribution of rms values for the situations p <dig>  p <dig> and p <dig>  in slip-resistance level  <dig>  are not significantly different among them .
however, p <dig> and p <dig> are significantly different when the slip-resistance level increased from  <dig> to  <dig>  p <dig> and p <dig> are not different statistically at slip-resistance level  <dig>  finally, p <dig>  p <dig> and p <dig> situations were statistically different at slip-resistance level  <dig>
the significant differences are more remarkable at slip-resistance level  <dig> for higher roughness .
however, these differences are not significant at slip-resistance levels  <dig> and  <dig>  probably because there are sudden shocks and stucks.
therefore, vibrissae movements are not repeatable, resulting in high dispersions of rms value distributions .
this effect is even more intensified when the slip-resistance level increases.
at slip-resistance level  <dig>  delta vibrissa sweeping on p <dig> and p <dig> results in low rms values due to improper movement of delta vibrissa.
since our experimental protocol allows us to identify the protraction and retraction phases, we have obtained the rms values in these two phases of the whisking cycle.
the distribution of rms values for each experimental condition and for the protraction phase is shown in figure 7a.
at slip-resistance level  <dig>  the statistically different situations were: p <dig> vs p600; p <dig> vs p180; p <dig> vs p180; p <dig> vs p <dig> and p <dig> vs p <dig>  in all these cases, p-values were lower than  <dig> .
the increase of slip-resistance level  produced statistical differences in p <dig> vs p <dig> and p <dig> vs p <dig>  which were not different at slip-resistance level  <dig>  at slip-resistance level  <dig> two significant differences were produced: p <dig> vs p <dig> and p <dig> vs p <dig>  table  <dig> shows the statistical results of these multiple comparisons.
these values are analyzed separately in protraction and retraction phases .
the parameters used in this analysis, such as diff of ranks, q and p, are described in the legend of table  <dig>
the information values obtained for the whole-whisking cycle  are shown in figure  <dig>  these values were calculated using three experimental situations .
the p <dig> and p <dig> situations  were not included in the calculation due to the differences in the kinematic characteristics of vibrissal movement described above.
the maximum amount of information was for slip-resistance level  <dig> .
this information value is given mainly by the difference between situations p <dig> vs p <dig> and p <dig> vs p <dig> .
even though situations p <dig>  p <dig> and p <dig> are statistically differents, the information value for slip-resistance level  <dig> is lower than for slip-resistance level  <dig>  this is due to that the quantitative differences, in slip-resistance level  <dig>  between p <dig> vs p <dig>  p <dig> vs p <dig> and p <dig> vs p <dig> do not exceed the differences between p <dig> vs p <dig> and p <dig> vs p <dig>  in slip-resistance level  <dig>  from this point of view, the information values provide a quantitative parameter of the differences between the sweep situations.
in the protraction phase, the information value was greater than in the whole-wisk in the slip-resistance level  <dig> .
in this phase, the information values were  <dig>  bits and  <dig>  bits for slip-resistance levels  <dig> and  <dig>  respectively .
in the retraction phase, the maximum amount of information was observed in slip-resistance level  <dig>  i =  <dig>  bits .
these results show that surfaces discrimination would require a strategy based on the slip-resistance levels variation.
furthermore, the discrimination of lower roughness has been achieved at higher slip-resistance levels, whereas the low slip-resistance levels are more appropriate for higher roughness.
discussion
texture encoding in the vibrissal system have being studied from the activity of trigeminal neurons and/or cortical neurons .
these responses usually consist of spikes or trains of spikes obtained from single or multi unit recordings  <cit> .
ideally, if the information from each constituent elements of the system is obtained, then the system would be fully described.
however, this would be practically impossible due to the complexity of necessary instruments.
the experimental procedures proposed here allow to observe  the system's global behaviour through a single recording .
in this paper the system's global behaviour refers to whole innervation of a vibrissal follicle.
contrary to the experimental protocols in which spikes are analyzed, this approach allows to makes inferences about activity of individual fibers.
here we analyzed compound action potentials  recorded at the vibrissal nerve.
the deep vibrissal nerve of the delta vibrissa consists of approximately  <dig> myelinated axons   <cit> .
these caps are the result of phase summation and cancellation of single fiber potentials with amplitudes that depend on fiber diameter, and the amplitude and shape of the cap is determined by the distribution of fiber diameters.
from the standpoint of digital processing these signals are considered to be stationary random process.
thus, for example, the temporal average would allow extracting the information from sensory nerves.
generally, the electrophysiological signals obtained with techniques of single-unit and multi-unit recordings are statistically weak.
it is easier to statistically analyze caps than single-unit recordings   <cit> .
we have proposed the rms method, because it presents temporal averaging.
another important aspect of this paper is the possibility to analyze the afferent activity in the preganglionic axons.
the preganglionic recordings show the afferents activity without pre-processing.
the hypothesis of this research is that the changes of afferent signal amplitude  would be related to different levels of mechanoreceptors activation  <cit> , and these last ones would be roughness surfaces dependent  <cit> .
as demonstrated by wolfe et al, each roughness pattern evokes kinetic signatures characterized by slip-stick events  <cit> .
these kinetic signatures would be encoded by spikes trains at both peripheral and central levels  <cit> .
thus, it is possible to speculate that changes in afferent signal amplitude may be related to temporal averages  and/or spatial averages   <cit> .
in this paper, the rms value was employed in the same way as they were used by albarracín et al  <cit> .
however, a more detailed analysis was performed.
the afferent activity was analyzed during the vibrissa protraction and retraction phases.
in addition, we calculated the maximum amount of information that was possible to extract from rms values.
we found that the lower roughness discrimination is achieved by using high slip-resistance levels .
on the other hand, higher roughness discrimination  is achieved by using the lower slip-resistance .
these results suggest that a behavioural strategy based on slip-resistance level changes could be used by rats in discrimination tasks.
we have also observed that as roughness increased, the rms values also increased in almost all cases.
these observations are consistent with those obtained by others authors  <cit> .
information values were calculated considering only the first three sweep situations , because the kinetic characteristics of vibrissal sweep over surfaces of greater roughness, changes dramatically with the increase of slip-resistance level.
in slip-resistance level  <dig>  irregularities in vibrissal sweep  were observed, thus avoiding complete the normal sweeping .
for this reason, high dispersions of rms values were observed .
under this slip-resistance level  <dig>  and during the sweeps on p <dig>  most of the sweeps were incomplete.
by increasing the slip-resistance level to  <dig>  the sweeps on p <dig> and p <dig>  were incomplete sweeps.
this particularity introduces unwanted factors in the evaluating of information when different slip-resistance levels are used.
to avoid this, the first three sweep situations were used to calculate the information.
under these considerations, information values variations with slip-resistance levels were observed, being the maximum value for the slip-resistance level  <dig> .
the information values obtained in protraction and retraction phases separately, shows that experimental situations are "better accurately" discriminated in retraction phase by using the slip-resistance level  <dig> .
authors such as arabzadeh et al.
<cit> , petersen et al.
<cit> , panzeri and diamond  <cit>  and others  <cit> , have used the information theory to investigate exactly which features of ensemble neuronal activity report information about whisker stimuli.
in other words, have investigated how sensory information is represented by neuronal populations in different stages of the sensory pathway.
unlike the researches mentioned above, here we used the information theory to demonstrate quantitatively the existence of texture information in the average afferent activity of a single vibrissa innervation and the dependence of this texture information with the slip-resistance.
methodological considerations
whisking in anesthetized rats
in this paper we induced artificial whisking in anesthetized rats and we recorded the afferent activity from vibrissal nerve during vibrissa sweeping on roughness surfaces.
although the principles of muscle-driven whisker movement and the basic patterns of movement trajectory are preserved in this kind of experimental protocols, the precise movements, mechanical interactions, and responses to touch may differ from those that occur during natural whisking  <cit> .
despite the inherent differences between natural and electrical whisking, the same kind of whisker motion was found in both experimental protocols  <cit> .
thus, experiments carried out in anesthetized animals, helps as the starting point for studies of multifiber activity in awake animals.
facial nerve stimulation - protocol
the vibrissa sweeping was produced by an electrical stimulation of the facial nerve.
we used square-wave pulse to simulate the vibrissal whisking at its natural frequency .
whisking movements during artificial whisking may differ from those produced during natural whisking in some details.
the vibrissa displacement angle observed in our experiments could differ from that of a behaving rat, as was previously described by carvell and simmons  <cit> .
however our experimental protocol allows to record and to analyze the afferent activity without stimulus artifacts.
the facial nerve stimulation with trains of pulses, as proposed by szwed et al  <cit> , it would contaminate the recording and could not be analyzed with current processing techniques.
CONCLUSIONS
the evidence of amplitude changes due to roughness surfaces and slip-resistance levels allows to speculate that texture information is slip-resistance dependent at peripheral level, although the amplitude of the afferent activity is not biologically a plausible neural coding schemes of texture.
finally, we have demonstrated that it is possible to characterize the texture information by analyzing the afferent activity of a single vibrissa innervation.
competing interests
the authors declare that they have no competing interests.
authors' contributions
fdf and ala designed and executed experimental protocols for acquisition of afferent signals, participated in the data interpretation and drafted the manuscript.
ala carried out the surgery procedures.
fdf conceived, designed and implemented the algorithms.
cjf assisted in digital processing and data interpretation.
all authors read and approved the manuscript.

 rodents as well as many mammals are characterized by the presence of vibrissae or whiskers located on both sides of the muzzle. when rats acquire sensory information by actively moving their vibrissae, a neural code is manifested at different levels of the sensory system. research in the past years has been extensive in the neural code at different stages of the whisker sensory pathway. arabzadeh et al, have demonstrated that a texture code would exist in peripheral afferent response of the vibrissal system and that spike rate would be the coding mechanism that underlies the textures discrimination in the primary somatosensory cortex. in a previous study we have shown that it is possible to characterize the vibrissal sensory information by analyzing the multifiber discharge of a selected vibrissal nerve. we showed that the afferent activity amplitude could be related to the physical properties of the surfaces. here, we used a similar protocol to that detailed by albarracín et al, but in this case we used different roughness surfaces. we have analyzed the afferent activity of a single vibrissa innervation in order to relate the information obtained to the roughness of surfaces. for this purpose, we have calculated the amplitude of afferent activity by using root mean square values obtained during protraction and retraction phases. here, we have also analyzed the effect of slip-resistance levels. in all cases, the differences among experimental conditions were quantified through the maximum amount of information by using information theory. we found that the discrimination of lower roughness could be achieved at higher slip-resistance levels, whereas the low slip-resistance levels are more appropriate for higher roughness. we also showed that the experimental conditions are better accurately discriminated in retraction phase. thus, we have demostrated that discrimination of rough surfaces would require a strategy based on the slip-resistance levels variation. methods procedures five wistar adult male rats, between and days, were used in our experiments. they were deeply anesthetized with urethane and their temperature was maintained at 37° by a servo-controlled heating pad. surgery consisted of exposing the infraorbital nerve as well as the two branches of the facial nerve on the right side. the motor branches were dissected and transected proximally to avoid possible motor influences on the sensorial pathway. the stimulation electrodes were placed on their distal stumps to produce the contraction of the mystacial muscles. the deep vibrissal nerve innervating a vibrissal follicle was identified with a high magnification of a dissecting microscope. the dissected nerve was also transected proximally and this action allowed eliminating discharges arriving from higher level of the sensorial pathway. to make sure that the nerve transection did not affect the functionality of the vibrissal nerve during our recording time, we test the falling of the nerve afferent activity throughout the time. we conclude that the activity start decreasing hour after the nerve section, so we never exceeded this space of time in our experiments. we used a bipolar electrode to record the multifiber afferent discharge of the vibrissal nerve selected. the recording electrodes as well as the nerves were immersed in a mineral oil bath during all recording. all these procedures were carry out in accordance with the recommendations of the guide for the care and use of laboratory animals. depth of anesthesia was ascertained and controlled during the experiment by the lack of either withdrawal reflex to hindlimb pinching or blink reflex to a gentle stimulation of the cornea. recording of the vibrissa electrical activity in this study we have recorded the multifiber activity of the delta vibrissal nerve while the vibrissa was sweeping surfaces of different roughness. the experimental protocol used in this paper has been previously described in detail by albarracín et al. the procedures are briefly described below. vibrissa movements were induced by electrical stimulation of facial motor nerve. square-wave pulses simulated vibrissal whisking at its natural frequency. a diagram of the experimental set up is shown in figure 1b. nerve activity was recorded and digitized at khz during a ms window following onset of each cycle of whisker movement. fifty whisker movement cycles were obtained for each surface, and an additional cycles were recorded while whisker moved unobstructed in air. three levels of slip-resistance were presented for each swept surface. these levels were established by mounting the surface at different distances from the whisker base. a minimal slip-resistance level was presented by placing the surface at a maximal distance from the whisker base so that the tip just barely contacted the surface throughout the entire movement cycle. increased slip-resistance levels were presented by moving the surface and mm closer to the whisker base. this procedure to increase the slip-resistance is similar to those described by albarracin et al. movements of the delta whisker were recorded simultaneously with nerve activity by using a custom-made photoresistive sensor. the frequency response of the sensor was maximal in the range 0- hz, enabling direct identification of the protraction and retraction phases of the movement cycle. rough surfaces the swept surfaces used in this paper were sandpapers of different grain size: p p p p and p. we measured the surfaces roughness by using a hommel tester t and we used the ra parameter as a roughness estimation. ra values obtained were:,,, and μm, respectively. digital processing and statistics the rms values and spectral estimation were previously used for the analysis of multifiber recordings. therefore, the theoretical details of both processing techniques can be found in the reference mentioned above. here, the rms values were obtained and analyzed in: the whole whisker movement, during the protraction and during the retraction. statistical analysis was carried out with kruskal wallis anova. the kruskal wallis test is a nonparametric alternative to the one-way anova. the test assesses the hypothesis that the different samples in the comparison were drawn from the same distribution or from distributions with the same median. thus, the interpretation of the kruskal-wallis anova is basically identical to the parametric one-way anova, except that it is based on ranks rather than on mean values. the dunn's method for multiple comparisons was used as post-hoc test, when anova repeat measurement revealed a significant difference. these analyses were made with software sigmastat http://www.systat.com/. information theory to determine the amount of information in a biological system it is necessary to have at least a pair of stimuli/responses situations. the stimulus may be a time series or simply belong to a class. the response depends on the characteristic of the signal that is being examined. thus, they can be real values or integer values. the information conveyed by the neuronal response about the stimulus can be quantified by shannon's mutual information formula, abbreviated hereafter as information: where p is the probability stimulus s presentation, p is the posterior probability of s given the observation of response r, and p is the probability of response r unconditional on the stimulus. the information determines the maximum amount of knowledge available to an observer who knows the posterior probabilities p and uses them to read off the signals available in a single observation of a spike train. for a better understanding of the procedure to determine the amount of information, we will consider a specific application in which the stimuli are given by the experimental situations. the responses will be real values, for example, the rms values of the afferent activities. these real values will be in the range. where a and b are the maximum and minimum values of the response. first, the frequency diagrams determination. for each experimental situation, the number of occurrences of the response is calculated for each bin. a graphical way to represent these diagrams is according to a histogram. the number of bins of the histograms will be calculated directly by the root square of the elements number. determination of the joint probability distribution, p. the following table is made. the response and stimuli are r and s, and their marginal probabilities distributions are p and p, respectively each element of table is the joint probability value, p. it is conceptually defined as follows: p, it is the probability of obtaining a response r = a and that it corresponds to the sweep on p this is: where: #occurrp a is the occurrences number obtained for a value of r = a. #total_occurr is the total number of occurrences in all situations and all possible responses. for frequency diagrams of figure the joint probability distribution is shown in figure determination of the p and p probability distributions. the probability of obtaining a response r, regardless of whether stimulus s did or not occur, is called the marginal probability, and it can be calculated by the sum of joint probabilities for a given response r. this is: similarly, we obtain the probability function p. determination of conditional probability distribution p. by definition of conditional probability, is given by: then, conditional probabilities distributions for each stimulus are obtained from joint and marginal probabilities distributions. thus, the following probabilities distributions are obtained: p, p,..., p. the conditional probability, p, conceptually answer the question: given a response r, what is the probability that this response has been produced by the p stimulus?. figure shows the conditional probability distributions for sweep situations. determination of the amount of information. after obtaining all probability distributions, it is possible to obtain the mutual information using eq. in this paper, we quantified the information considering the rms values as biological responses of the system. for the rms values, the information was calculated for each slip-resistance level, as well as for each whisking phases.  the recordings obtained are the average electric activity of about myelinated axons which supply innervation to the delta vibrissa. figure 5b shows the afferent discharge recorded during different sweep situations. it is possible to distinguish amplitude changes in the signals in relation with the surface roughness. this feature is more appreciable in figures 5c and 5d. whisker displacements and evoked afferent activity were simultaneously recorded for all cases. figure 5a shows the afferent activity evoked by vibrissal sweeps in the air and its corresponding displacement. vibrissal movement recordings were acquired with a bandwidth of 0- hz; thus, only slow movements were monitored. therefore, displacement recordings were only used to identify the protraction and retraction phases. we calculated one rms value for each recording obtained during a vibrissa sweeping. each experimental condition consists of rms values. we first calculated the rms values for the whole-whisk recordings. the distribution of rms values for each sweep situation is shown in figure the test computes statistic q, the number of rank sums, and shows whether p < or not, for the pair that are being compared. p is the probability that the null hypothesis may be rejected and, thus, it helps conclude that there are differences between treatments. diff of ranks is the difference in the rank sum orders that are being compared. the rank sums are a measurement of the difference between two treatments. the distribution of rms values for the situations p p and p in slip-resistance level are not significantly different among them. however, p and p are significantly different when the slip-resistance level increased from to p and p are not different statistically at slip-resistance level finally, p p and p situations were statistically different at slip-resistance level the significant differences are more remarkable at slip-resistance level for higher roughness. however, these differences are not significant at slip-resistance levels and probably because there are sudden shocks and stucks. therefore, vibrissae movements are not repeatable, resulting in high dispersions of rms value distributions. this effect is even more intensified when the slip-resistance level increases. at slip-resistance level delta vibrissa sweeping on p and p results in low rms values due to improper movement of delta vibrissa. since our experimental protocol allows us to identify the protraction and retraction phases, we have obtained the rms values in these two phases of the whisking cycle. the distribution of rms values for each experimental condition and for the protraction phase is shown in figure 7a. at slip-resistance level the statistically different situations were: p vs p600; p vs p180; p vs p180; p vs p and p vs p in all these cases, p-values were lower than. the increase of slip-resistance level produced statistical differences in p vs p and p vs p which were not different at slip-resistance level at slip-resistance level two significant differences were produced: p vs p and p vs p table shows the statistical results of these multiple comparisons. these values are analyzed separately in protraction and retraction phases. the parameters used in this analysis, such as diff of ranks, q and p, are described in the legend of table the information values obtained for the whole-whisking cycle are shown in figure these values were calculated using three experimental situations. the p and p situations were not included in the calculation due to the differences in the kinematic characteristics of vibrissal movement described above. the maximum amount of information was for slip-resistance level. this information value is given mainly by the difference between situations p vs p and p vs p. even though situations p p and p are statistically differents, the information value for slip-resistance level is lower than for slip-resistance level this is due to that the quantitative differences, in slip-resistance level between p vs p p vs p and p vs p do not exceed the differences between p vs p and p vs p in slip-resistance level from this point of view, the information values provide a quantitative parameter of the differences between the sweep situations. in the protraction phase, the information value was greater than in the whole-wisk in the slip-resistance level. in this phase, the information values were bits and bits for slip-resistance levels and respectively. in the retraction phase, the maximum amount of information was observed in slip-resistance level i = bits. these results show that surfaces discrimination would require a strategy based on the slip-resistance levels variation. furthermore, the discrimination of lower roughness has been achieved at higher slip-resistance levels, whereas the low slip-resistance levels are more appropriate for higher roughness. discussion texture encoding in the vibrissal system have being studied from the activity of trigeminal neurons and/or cortical neurons. these responses usually consist of spikes or trains of spikes obtained from single or multi unit recordings. ideally, if the information from each constituent elements of the system is obtained, then the system would be fully described. however, this would be practically impossible due to the complexity of necessary instruments. the experimental procedures proposed here allow to observe the system's global behaviour through a single recording. in this paper the system's global behaviour refers to whole innervation of a vibrissal follicle. contrary to the experimental protocols in which spikes are analyzed, this approach allows to makes inferences about activity of individual fibers. here we analyzed compound action potentials recorded at the vibrissal nerve. the deep vibrissal nerve of the delta vibrissa consists of approximately myelinated axons. these caps are the result of phase summation and cancellation of single fiber potentials with amplitudes that depend on fiber diameter, and the amplitude and shape of the cap is determined by the distribution of fiber diameters. from the standpoint of digital processing these signals are considered to be stationary random process. thus, for example, the temporal average would allow extracting the information from sensory nerves. generally, the electrophysiological signals obtained with techniques of single-unit and multi-unit recordings are statistically weak. it is easier to statistically analyze caps than single-unit recordings. we have proposed the rms method, because it presents temporal averaging. another important aspect of this paper is the possibility to analyze the afferent activity in the preganglionic axons. the preganglionic recordings show the afferents activity without pre-processing. the hypothesis of this research is that the changes of afferent signal amplitude would be related to different levels of mechanoreceptors activation, and these last ones would be roughness surfaces dependent. as demonstrated by wolfe et al, each roughness pattern evokes kinetic signatures characterized by slip-stick events. these kinetic signatures would be encoded by spikes trains at both peripheral and central levels. thus, it is possible to speculate that changes in afferent signal amplitude may be related to temporal averages and/or spatial averages. in this paper, the rms value was employed in the same way as they were used by albarracín et al. however, a more detailed analysis was performed. the afferent activity was analyzed during the vibrissa protraction and retraction phases. in addition, we calculated the maximum amount of information that was possible to extract from rms values. we found that the lower roughness discrimination is achieved by using high slip-resistance levels. on the other hand, higher roughness discrimination is achieved by using the lower slip-resistance. these results suggest that a behavioural strategy based on slip-resistance level changes could be used by rats in discrimination tasks. we have also observed that as roughness increased, the rms values also increased in almost all cases. these observations are consistent with those obtained by others authors. information values were calculated considering only the first three sweep situations, because the kinetic characteristics of vibrissal sweep over surfaces of greater roughness, changes dramatically with the increase of slip-resistance level. in slip-resistance level irregularities in vibrissal sweep were observed, thus avoiding complete the normal sweeping. for this reason, high dispersions of rms values were observed. under this slip-resistance level and during the sweeps on p most of the sweeps were incomplete. by increasing the slip-resistance level to the sweeps on p and p were incomplete sweeps. this particularity introduces unwanted factors in the evaluating of information when different slip-resistance levels are used. to avoid this, the first three sweep situations were used to calculate the information. under these considerations, information values variations with slip-resistance levels were observed, being the maximum value for the slip-resistance level. the information values obtained in protraction and retraction phases separately, shows that experimental situations are "better accurately" discriminated in retraction phase by using the slip-resistance level. authors such as arabzadeh et al., petersen et al., panzeri and diamond and others, have used the information theory to investigate exactly which features of ensemble neuronal activity report information about whisker stimuli. in other words, have investigated how sensory information is represented by neuronal populations in different stages of the sensory pathway. unlike the researches mentioned above, here we used the information theory to demonstrate quantitatively the existence of texture information in the average afferent activity of a single vibrissa innervation and the dependence of this texture information with the slip-resistance. methodological considerations whisking in anesthetized rats in this paper we induced artificial whisking in anesthetized rats and we recorded the afferent activity from vibrissal nerve during vibrissa sweeping on roughness surfaces. although the principles of muscle-driven whisker movement and the basic patterns of movement trajectory are preserved in this kind of experimental protocols, the precise movements, mechanical interactions, and responses to touch may differ from those that occur during natural whisking. despite the inherent differences between natural and electrical whisking, the same kind of whisker motion was found in both experimental protocols. thus, experiments carried out in anesthetized animals, helps as the starting point for studies of multifiber activity in awake animals. facial nerve stimulation - protocol the vibrissa sweeping was produced by an electrical stimulation of the facial nerve. we used square-wave pulse to simulate the vibrissal whisking at its natural frequency. whisking movements during artificial whisking may differ from those produced during natural whisking in some details. the vibrissa displacement angle observed in our experiments could differ from that of a behaving rat, as was previously described by carvell and simmons. however our experimental protocol allows to record and to analyze the afferent activity without stimulus artifacts. the facial nerve stimulation with trains of pulses, as proposed by szwed et al, it would contaminate the recording and could not be analyzed with current processing techniques.  the evidence of amplitude changes due to roughness surfaces and slip-resistance levels allows to speculate that texture information is slip-resistance dependent at peripheral level, although the amplitude of the afferent activity is not biologically a plausible neural coding schemes of texture. finally, we have demonstrated that it is possible to characterize the texture information by analyzing the afferent activity of a single vibrissa innervation. competing interests the authors declare that they have no competing interests. authors' contributions fdf and ala designed and executed experimental protocols for acquisition of afferent signals, participated in the data interpretation and drafted the manuscript. ala carried out the surgery procedures. fdf conceived, designed and implemented the algorithms. cjf assisted in digital processing and data interpretation. all authors read and approved the manuscript. 
the evidence of amplitude changes due to roughness surfaces and slip-resistance levels allows to speculate that texture information is slip-resistance dependent at peripheral level. here, we analyzed the multifiber afferent discharge of one vibrissal nerve during active whisking. vibrissae movements were induced by electrical stimulation of motor branches of the facial nerve. as roughness increased, the rms values also increased in almost all cases. in addition, we observed a better discrimination performance in the retraction phase. the amplitude of afferent activity was analyzed according to its root mean square value. we used sandpapers of different grain size as roughness discrimination surfaces and we also consider the change of vibrissal slip-resistance as a way to improve tactile information acquisition. we found that the change of the vibrissal slip-resistance is a way to improve the roughness discrimination of surfaces. in particular, we focus on recent evidence of how neurons at different levels of the sensory pathway carry information about tactile stimuli. the comparisons among experimental situation were quantified by using the information theory.
 studies in tactile discrimination agree that rats are able to learn a rough-smooth discrimination task by actively touching objects with their vibrissae. in particular, we focus on recent evidence of how neurons at different levels of the sensory pathway carry information about tactile stimuli. here, we analyzed the multifiber afferent discharge of one vibrissal nerve during active whisking. vibrissae movements were induced by electrical stimulation of motor branches of the facial nerve. we used sandpapers of different grain size as roughness discrimination surfaces and we also consider the change of vibrissal slip-resistance as a way to improve tactile information acquisition. the amplitude of afferent activity was analyzed according to its root mean square value. the comparisons among experimental situation were quantified by using the information theory.  we found that the change of the vibrissal slip-resistance is a way to improve the roughness discrimination of surfaces. as roughness increased, the rms values also increased in almost all cases. in addition, we observed a better discrimination performance in the retraction phase.  the evidence of amplitude changes due to roughness surfaces and slip-resistance levels allows to speculate that texture information is slip-resistance dependent at peripheral level. 