it is to be understood that the figures and descriptions of the present invention have been simplified to illustrate elements that are relevant for a clear understanding of the invention , while eliminating , for purpose of clarity , other elements that may be well known . the detailed description will be provided herein below with reference to the attached drawings . the term “ subject ” in this application refers to both animals and humans . the term “ stable rem ” refers visually to the bottom portion of the pattern as in the bimodal distribution of rem . the term “ unstable rem ” refers visually to the top portion of the pattern in the bimodal distribution of rem . the methods described herein are disclosed in detail in pct / us2006 / 018120 , the disclosure is fully incorporated herein by reference . the present invention provides a system and method to obtain and classify eeg data in both animals and humans . obtained eeg signals are low - power frequency signals and follow a 1 / f distribution , whereby the power in the signal is inversely related , e . g ., inversely proportional , to the frequency . eeg signals have typically been examined in time in series increments called epochs . for example , when the eeg signal is used for analyzing sleep , sleep may be segmented into one or more epochs to use for analysis . the epochs can be segmented into different sections using a scanning window , where the scanning window defines different sections of the time series increment . the scanning window can move via a sliding window , where sections of the sliding window have overlapping time series sequences . an epoch can alternatively span an entire time series , for example . according to the present application , different forms of sleep state of a subject may be monitored . a sleep state is described as any distinguishable sleep or wakefulness that is representative of behavioral , physical or signal characteristics . sleep states which are referred to in this application include slow wave sleep or sws , rapid eye movement sleep or rem , intermediate sleep states also called inter or is states , and awake states . awake states may actually be part of the sleep state , and the awake states can be characterized by vigilance into attentiveness or levels of alertness . the intermediate sleep can also be characterized as intermediate - 1 sleep and intermediate - 2 sleep . an artifact may also be obtained during acquisition of an eeg . an artifact is data that misrepresents the eeg . for example , movement within a user that registers on the eeg may be an artifact . example artifacts include muscle twitches and the like . referring now to fig1 which is a flow diagram of an exemplary system 100 for determining sleep state information of a subject . the eeg data 102 is received from the subject . in any of the embodiments described herein , a variety of source data can be analyzed including electroencephalography ( eeg ) data , electrocardiography data ( ekg ), electrooculography data ( eog ), electrocorticographic ( ecog ) data , intracranial data , electromyography data ( emg ), local field potential ( lfp ) data , magnetoencephalograhic data ( meg ), spike train data , wave data including sound and pressure waves , and any data exhibiting where there are differences in dynamic range of power for various frequencies across a frequency spectrum of the data e . g ., a 1 / f distribution . source data can include encoded data stored at low power frequency within source data . in one embodiment of the invention , the data 102 once received from the subject is transmitted to a software program 104 for analysis . exemplary system for determining low - power frequency information from source data with at least one low power frequency range source data 102 with at least one low power frequency range is obtained and input into software 104 to determine low power frequency information . source date with at least one low power frequency range 102 is received . for example , electroencephalography source data for a subject can be received . source data can be received via a single channel or multiple channels . in a preferred embodiment of this invention a single channel of eeg was sufficient to decouple sleep and waking states . source data is adjusted to increase the dynamic range for power within at least one low power frequency range of the frequency spectrum of the source data as compared to a second higher power frequency range . a number of adjustment techniques described herein , including normalization and frequency weighting can be used . in an embodiment , electroencephalography source data is normalized to increase the low power , higher frequency range data relative to the higher power , lower frequency range data or , more generally , to normalize the powers of the different signal parts . after the source data is adjusted , various other processing can be done . for example , a visualization of the adjusted source data can be presented . further , low power frequency information can be extracted from the adjusted source data . for example , low power frequency information can be extracted from adjusted electroencephalography source data . higher power frequency information can also be extracted from the adjusted source data . the method described in this or any of the other examples can be a computer - implemented method performed via computer - executable instructions in one or more computer - readable media . any of the actions shown can be performed by software incorporated within a signal processing system or any other signal data analyzer system . referring again to fig1 , electroencephalography data 102 for a subject is obtained and input into software 104 to determine sleep state information for the subject 106 . the software can employ any combination of technologies , such as those described herein , to determine sleep state information for the subject . referring now to fig2 , a block diagram of an exemplary system 200 for determining sleep states of a subject wherein the data can be normalized to compute a spectrogram 202 . another embodiment uses multiple normalizations for even further dynamic range increase . normalizations can be performed by normalizing frequency across time or time across frequency . exemplary method for adjusting source data to account for differences in power over a spectrum of frequencies over time for example , electroencephalography data with at least one low power frequency range can be received . artifacts in the data can be removed from the source data . for example , artifact data can be manually removed from the source data or automatically filtered out of source data via a filtering ( e . g ., dc filtering ) or data smoothing technique . the source data can also be pretreated with component analysis 204 . the source data is segmented into one or more epochs ; where each epoch is a portion of data from the series . for example , the source data can be segmented into a plurality of time segments via a variety of separating techniques . scanning windows and sliding windows can be used to separate the source data into time series increments . the one or more epochs are normalized for differences in power of the one or more epochs across time . for example , the power of each epoch at one or more frequencies can be normalized across time to determine appropriate frequency windows for extracting information . such normalization can reveal low power , statistically significant shifts in power at one or more frequencies ( e . g ., delta , gamma , and the like ). any frequency range can be revealed and utilized for analysis . information can be calculated for each of the one or more epochs after appropriate frequency windows have been established . such information can include low frequency power ( e . g ., delta power ), high frequency power ( e . g ., gamma power ), standard deviation , maximum amplitude ( e . g ., maximum of the absolute value of peaks ) and the sort . further calculations can be done on the information calculated for each of the one or more epochs creating information such as gamma power / delta power , time derivative of delta , time derivative of gamma power / delta power and the like . time derivatives can be computed over preceding and successive epochs . after calculating the information , that information can then be normalized across the one or more epochs . a variety of data normalization 202 techniques can be conducted including z - scoring and other similar techniques . results of the adjustment of source data to account for differences in power over a spectrum of frequencies over time can be presented as one or more epochs of data . for example , frequency weighted epochs can be presented as adjusted source data . electroencephalography data for a subject is obtained and input into segmenter to segment the data into one or more epochs . in practice , epochs are of similar ( e . g ., the same ) length . epoch length can be adjusted via a configurable parameter . the one or more epochs , in turn , are input into normalizer 202 to normalize frequency data in the one or more epochs across time , thereby frequency weighting the one or more epochs of electroencephalography data . the one or more frequency weighted epochs are then input into classifier to classify the data into sleep states , thereby generating sleep state information for the subject 208 . methods for determining sleep state information for a subject are described in detail below . electroencephalography ( eeg ) data for a subject is received . for example , electroencephalography data , which exhibits lower dynamic range for power in at least one low power first frequency range in a frequency spectrum as compared to a second frequency range in the frequency spectrum , can be received . the electroencephalography data for the subject is segmented into one or more epochs . for example , the eeg data can be segmented into one or more epochs via a variety of separating techniques . scanning windows and sliding windows can be used to separate the eeg data into one or more epochs . the source data can also be filtered via direct current tiitermg during , prior to , or after segmenting . the source data can also be pretreated with component analysis 204 ( e . g ., principle or independent component analysis ). in entire night eeg data the higher frequencies ( e . g ., gamma ) exhibit lower power than the lower frequencies ( e . g ., delta , theta and the like ) in the whole night eeg data . frequency power of the one or more epochs is weighted across time . for example , the power of each epoch at one or more frequencies can be normalized 202 across time to determine appropriate frequency windows for extracting information . such normalization can reveal low power , statistically significant shifts in power at one or more frequencies ( e . g ., delta , gamma , and the like ). additionally , each epoch can be represented by the frequency with the highest relative power over time to determine appropriate frequency windows for extracting information . alternatively , component analysis ( e . g ., principle component analysis ( pca ) or independent component analysis ( ica )) 204 can be utilized after normalization 202 to further determine appropriate frequency windows for extracting information . any frequency range can be revealed and utilized for analysis . information can be calculated for each of the one or more epochs after appropriate frequency windows have been established ( e . g ., after weighting frequency ). such information can include low frequency power ( e . g ., delta power ), high frequency power ( e . g ., gamma power ), standard deviation , maximum amplitude ( e . g ., maximum of the absolute value of peaks ) and the sort . further calculations can be done on the information calculated for each of the one or more epochs creating information such as gamma power / delta power , time derivative of delta , time derivative of gamma power / delta power and the like . time derivatives can be computed over preceding and successive epochs . after calculating the information , it can then be normalized across the one or more epochs . a variety of data normalization techniques can be conducted including z - scoring and the like . the higher frequency data is now more clearly visible . sleep states 208 in the subject are classified based on the one or more frequency weighted epochs . for example , the one or more frequency weighted epochs can be clustered 206 by any variety of clustering techniques including k - means clustering . the clustering can be done on information calculated from the epochs ( e . g ., delta power , gamma power , standard deviation , maximum amplitude ( gamma / delta ), time derivative of delta , time derivative - of ( gamma / delta , and the sort ). component analysis ( e . g ., pca or ica ) can be used to determine the parameter space ( e . g ., types of information used ) in the clustering . subsequent to clustering 206 , sleep state designations can be assigned to the epochs . sleep state designated epochs can then be presented as representations of sleep states in the subject for the period of time represented by the epoch . classification can also incorporate manually determined sleep states ( e . g ., manually determined “ awake ” versus “ sleeping ” sleep states ). additionally , artifact information ( e . g . movement data , poor signal data , or the like ) can be utilized in the classification . epochs can be classified according to the sleep states they represent . an epoch can be classified according to normalized variables ( e . g ., information calculated for an epoch ) based on high frequency information , low frequency information , or both high and low frequency information . for example , rem sleep state epochs can have higher relative power than sws at higher frequencies and lower relative power than sws at lower frequencies . similarly , sws sleep state epochs can have lower relative power than rem at higher frequencies and higher relative power than rem at lower frequencies . additionally , epochs initially classified as both nrem and nsws sleep ( e . g ., epochs having low relative power at both higher and lower frequencies ) can be classified as intermediate sleep and epochs classified as both rem and sws sleep ( e . g ., epochs having high relative power at both higher and lower frequencies ) can be classified as outliers . further , epochs initially classified as both nrem and nsws sleep can be classified as intermediate stage i sleep and epochs initially classified as both rem and sws sleep can be classified as intermediate stage ii sleep . additionally , sleep states can be split in the classifying to look for spindles , k - complexes , and other parts . any group of epochs initially classified as one sleep state can be split into multiple sub - classified sleep states according to increasing levels of classification detail . for example , a group of epochs classified as sws can be reclassified as two distinct types of sws . artifact data ( e . g . movement data , poor signal data , and the like ) can also be used in sleep state classification . for example , artifacts can be used to analyze whether epochs initially assigned a sleep state designation should be reassigned a new sleep state designation due to neighboring artifact data . for example , an epoch assigned a sleep state designation of rem that has a preceding movement artifact or awake epoch can be reassigned a sleep state designation of awake . further , for example , an artifact epoch that has a succeeding sws epoch can be reassigned a sleep state designation of sws because there is a high likelihood that the epoch represents a large sws sleep epoch rather than a large movement artifact which is more common during wakefulness . in such ways , for example , artifact data can be utilized in a data smoothing technique . any variety of data smoothing techniques can be used during the assigning of sleep states . for example , numbers ( e . g ., 0 and 1 ) can be used to represent designated sleep states . neighboring epochs &# 39 ; sleep state designation numbers can then be averaged to determine if one of the epochs is inaccurately assigned a sleep state designation . for example , abrupt jumps from sws - nsws - sws ( and rem - nrem - rem ) are rare in sleep data . therefore , should a group of epochs be assigned sleep state designations representing abrupt jumps in sleep states , smoothing techniques can be applied to improve the accuracy of the assigning . referring now to fig3 , a block diagram of an exemplary system 300 for determining sleep states of a subject . the data is received from the subject 302 either manually or automatically . the preferred frequency analysis , temporal fragmentation or spectral fragmentation 304 can be performed on the data in order to determine at least one parameter of sleep . this information can be further classified to determine a sleep state 306 . previous embodiments have shown how normalization , for example using z scoring , allowed analysis of more information from the brainwave signal . the analysis which was previously carried out normalized power information across frequencies . the normalization preferably used z scoring , but any other kind of data normalization can be used . the normalization which is used is preferably unitless , like z scoring . as well - known in the art , z scoring can be used to normalize a distribution without changing a shape of the envelope of the distribution . the z scores are essentially changed to units of standard deviation . each z score normalized unit reflects the amount of power in the signal , relative to the average of the signal . the scores are converted into mean deviation form , by subtracting the mean from each score . the scores are then normalized relative to standard deviation . all of the z scored normalized units have standard deviations that are equal to unity . while the above describes normalization using z scores , it should be understood that other normalizations can also be carried out , including t scoring , and others . multiple normalizations may also be employed . normalizations can be performed by normalizing frequency across time or time across frequency . the above embodiments describe normalizing the power at every frequency within a specified range . the range may be from 0 , to 100 hz , or to 128 hz , or to 500 hz . the range of frequencies is only restricted by the sampling rate . with an exemplary sampling rate of 30 khz , an analysis up to 15 khz can be done . according to the present embodiment , additional normalizations are carried out which normalizes the power across time for each frequency . this results in information which has been normalized across frequencies and across time being used to create a normalized spectrogram . this embodiment can obtain additional information from brainwave data , and the embodiment describes automatically detecting different periods of sleep from the analyzed data . the periods of sleep that can be detected can include , but are not limited to , short wave sleep ( sws ), rapid eye movement sleep ( rem ), intermediate sleep ( iis ) and wakefulness . according to an important feature , a single channel of brainwave activity ( that is obtained from a single location on the human skull ) is used for the analysis . as described above , the obtained data can be one channel of eeg information from a human or other subject . the eeg data as obtained can be collected , for example , using a 256 hz sampling rate , or can be sampled at a higher rate . the data is divided into epochs , for example 30 second epochs , and characterized according to frequency . a first frequency normalization is carried out . the power information is normalized using a z scoring technique on each frequency bin . in the embodiment , the bins may extend from one to 100 hz and 30 bins per hertz . the normalization occurs across time . this creates a normalized spectrogram or ns , in which each frequency band from the signal has substantially the same weight . in the embodiment , each 30 second epoch is represented by a “ preferred frequency ” which is the frequency with the largest z score within that epoch . this creates a special frequency space called the preferred frequency space . analysis of how those patterns are formed and analysis of the characteristics of the patterns can be done . different sleep states , therefore , can be defined according to a discrimination function , where the discrimination function looks for certain activity in certain areas , and non - activity in other areas . the function may evaluate sleep states according to which of the frequency at areas have activity and which do not have activity . more generally , however , any form of dynamic spectral scoring can be carried out on the compensated data . the discrimination function may require specific values , or may simply require a certain amount of activity to be present or not present , in each of a plurality of frequency ranges . the discrimination function may simply match envelopes of frequency response . the discrimination function may also look at spectral fragmentation and temporal fragmentation . a second normalization which is carried out across frequencies . the second normalization produces a doubly normalized spectrogram . this produces a new frequency space , in which the bands become even more apparent . the doubly normalized spectrogram values can be used to form filters that maximally separate the values within the space . a clustering technique which is carried out on the doubly normalized frequency . for example , the clustering technique may be a k means technique as described in the previous embodiments . each cluster can represent a sleep state . the clusters are actually multi dimensional clusters , which can themselves be graphed to find additional information . the number of dimensions can depend on the number of clustering variables . this illustrates how the doubly normalized spectrogram also allows many more measurement characteristics . measurement of the average spread in normalized power across frequency which illustrates the spectral fragmentation is also possible . fragmentation values can alternatively be based on temporal fragmentation for the different states may also be used as part of the discrimination function . these two functions are evaluated on the doubly normalized spectrum , relying on homogeneous increases in gain at all frequencies as caused movement artifacts in nrem sleep and w would lead to abnormally elevated fragmentation values in the singly normalized spectrum . these fragmentation values may be used as part of the discrimination function . importantly , and as described above , this discrimination function is typically not apparent from any previous analysis technique , including manual techniques . the computation may be characterized by segmenting , or may use overlapping windows or a sliding window , to increase the temporal registration . this enables many techniques that have never been possible before . by characterizing on - the - fly , this enables distinguishing using the dynamic spectral scoring , between sleep states and awake states using the brainwave signature alone . the exemplary methods for data analysis described above were combined with a standard non - invasive eeg method for humans . the result is the ability to non - invasively extract attenuated rhythms in animals , automatically analyze the brain activity from a single channel of eeg , and sufficiently classify the sleep parameters for the animals . rats were anesthetized with isoflurane . the scalp was gently shaved . conductive electrogel was applied and a standard 6 mm gold plated electrode was secured with collodion . the resulting data were analyzed using advanced computational techniques , which are described above , by using software and techniques described in p . c . t . application wo2006 / 1222201 . voltage signal from the rat brain is collected by the electrodes and sent to the computer for analysis . the signal is broken down into roughly three second epochs of signal . the frequency spectra for each epoch are calculated to produce a whole recording spectrum . the resulting spectrum is then normalized across frequencies which allows for the detection of previously unidentified frequencies . at each time epoch , only the frequency with the highest shift with respect to the baseline is mapped . the resulting map shows different signatures in this space relative to the baseline . referring again to fig2 , these signatures can be used to create variables used on a multiply normalized ( normalizations across time and frequency ) spectrogram 202 to create a parameter space to separate stages . component analysis 204 can also be used on the multiply normalized spectrogram to create clusters 206 . there are a wide variety of clustering and classification methods used in computational signal processing to differentiate data into distinct classes . as described herein , the clustering method used is k - means clustering but any computational signal processing method for differentiating groups of data could be used . similarly , classification methods such as component analysis ( e . g ., principle and independent component analysis ) are used as described herein . clustering ( or cluster analysis ) is unsupervised learning where the classes are unknown a priori and the goal is to discover these classes from data . for example , the identification of new tumor classes using gene expression profiles is a form of unsupervised learning . classification ( or class prediction ) is a supervised learning method where the classes are predefined and the goal is to understand the basis for the classification from a set of labeled objects and build a predictor for future unlabeled observations . for example , the classification of malignancies into known classes is a form of supervised learning . clustering procedures commonly fall into two categories : hierarchical methods and partitioning methods . hierarchical methods can be either divisive ( top - down ) or agglomerative ( bottom - up ). hierarchical clustering methods produce a tree or dendrogram . hierarchical methods provide a hierarchy of clusters , from the smallest , where all objects are in one cluster , through to the largest set , where each observation is in its own cluster partitioning methods usually require the specification of the number of clusters . then , a mechanism for apportioning objects to clusters must be determined . these methods partition the data into a prespecified number k of mutually exclusive and exhaustive groups . the method iteratively reallocates the observations to clusters until some criterion is met ( e . g . minimize within - cluster sums of - squares ). examples of partitioning methods include k - means clustering , partitioning around medoids ( pam ), self organizing maps ( som ), and model - based clustering . most methods used in practice are agglomerative hierarchical methods , in a large part due to the availability of efficient exact algorithms . however both clustering methods have their advantages and disadvantages . hierarchical advantages include fast computation , at least for agglomerative clustering , and disadvantages include that they are rigid and cannot be corrected later for erroneous decisions made earlier in the method . partitioning advantages include that such methods can provide clusters that ( approximately ) satisfy an optimality criterion , and disadvantages include that one needs an initial k and the methods can take long computation time . in summary , clustering is a more difficult problem than classifying for a variety of reasons including the following : there is no learning set of labeled observations the number of groups is usually unknown implicitly , one must have already selected both the relevant features and distance measures used in clustering methods . techniques involving statistics , machine learning , and psychometrics can be used . examples of classifiers include logistic regression , discriminant analysis ( linear and quadratic ), principle component analysis ( pca ), nearest neighbor classifiers ( k - nearest neighbor ), classification and regression trees ( cart ), prediction analysis for microarrays , neural networks and multinomial log - linear models , support vector machines , aggregated classifiers ( bagging , boosting , forests ), and evolutionary algorithms . logistic regression is a variation of linear regression which is used when the dependent ( response ) variable is a dichotomous variable ( i . e ., it takes only two values , which usually represent the occurrence or non - occurrence of some outcome event , usually coded as 0 or 1 ) and the independent ( input ) variables are continuous , categorical , or both . for example , in a medical study , the patient survives or dies , or a clinical sample is positive or negative for a certain viral antibody . unlike ordinary regression , logistic regression does not directly model a dependent variable as a linear combination of dependent variables , nor does it assume that the dependent variable is normally distributed . logistic regression instead models a function of the probability of event occurrence as a linear combination of the explanatory variables . for logistic regression , the function relating the probabilities to the explanatory variables in this way is the logistic function , which has a sigmoid or s shape when plotted against the values of the linear combination of the explanatory variables . logistic regression is used in classification by fitting the logistic regression model to data and classifying the various explanatory variable patterns based on their fitted probabilities . classifications of subsequent data are then based on their covariate patterns and estimated probability ] discriminant analysis : in summary discriminant analysis represents samples as points in space and then classifies the points . linear discriminant analysis ( lda ) finds an optimal plane surface that best separates points that belong to two classes . quadratic discriminant analysis ( qda ) finds an optimal curved ( quadratic ) surface instead . both methods seek to minimize some form of classification error . lda finds linear combinations ( discriminant variables ) of data with large ratios of between - groups to within - groups sums of squares and predicts the class of an observation x by the class whose mean vector is closest to x in terms of the discriminant variables . advantages of lda include that it is simple and intuitive where the predicted class of a test case is the class with the closest mean and it is easy to implement with a good performance in practice . nearest neighbor methods are based on a measure of distance between observations , such as the euclidean distance or one minus the correlation between two data sets . k - nearest neighbor classifiers work by classifying an observation x as follows : find the k observations in the learning set that are closest to x . predict the class of x by majority vote , i . e ., choose the class that is most common among these k neighbors . simple classifiers with k = 1 can generally be quite successful . a large number of irrelevant or noise variables with little or no relevance can substantially degrade the performance of a nearest neighbor classifier . referring now to fig4 , an exemplary system for determining sleep states for a subject utilizing either automated data or manual data 400 . automated date 402 as well as manually scored data 404 can be used to compute the spectrogram 406 . the methods described above can be applied to analyze the data 408 and subsequently determine sleep state information for the subject . example 2 illustrates how the exemplary methods can be applied to determine sleep patterns from a single channel of eeg using either automated or manual data . one channel of eeg ( c3 - a2 derivation ) from twenty - six nights ( 8 hours each ) of sleep was obtained from twenty - six different polysomnographic recordings conducted in twenty - six healthy human subjects . the eeg data and manual scoring was provided by the experimental procedures were approved by the institutional review boards at each institution . eeg data were collected at 256 hz and bandpassed at 0 . 3 - 100 hz with a 60 hz notch filter ( ucsd ) or collected at 250 hz and bandpassed at 0 . 53 - 70 hz ( mpi ). these recordings were amplified at 10 k and manually scored in 30 sec epochs in accordance with r - k . for each recording , the whole night spectrogram was computed over 2 orthogonal tapers on 30 sec epochs using a standard multitaper technique . the power information was then normalized by z - scoring for each frequency bin ( from 1 to 100 hz , 30 bins per hz ) across time . this normalized spectrogram ( ns ) weighed each frequency band equally . each 30 second segment was represented by the frequency with the largest z - score . in this preferred frequency space ( pfs ), sleep and waking states broadly separated into different patterns ( fig2 , 22 .) w was always characterized by a band in alpha ( 7 - 12 hz ) and sometimes by a band in beta ( 15 - 25 hz ). is exhibited prominent activity in the spindle frequencies ( 12 - 15 hz ). surprisingly , rem was defined by compact bands in theta ( 4 - 8 hz ) and sometimes beta ( 15 - 25 hz ) frequencies whereas sws was dominated by delta activity . when computed over overlapping 3 sec windows and a 1 sec sliding window , similar trends were visible in the pfs except that beta activity emerged in rem . at that resolution , rem appears more “ awake - like ” than at a 30 sec resolution . however , at that resolution , all the sleep states whether they were identified manually or automatically had distinct signatures in the preferred frequency space . at each time point , z - scoring the normalized spectrogram across frequencies creates a doubly normalized spectrogram . in this space , bands apparent in the pfs still had positive values whereas dark regions tended to have negative values . by adding the double normalized spectrogram values of frequencies that show up as bands in the pfs and subtracting those that do not , filters can be constructed that maximally separate states . one maximizes w (‘ w filter ”), another separates nrem from w and rem (‘ nrem filter ’) and a third distinguishes is from sws (‘ sws filter ’). the output of these three filters spans a space in which the three broad sleep stages and w tend to separate . interestingly , stage i did not cluster in either space and sws formed only one cluster ( rather than two , one for stage iii and one for stage iv ). the latter is in accordance with the recent revision of r - k which abandoned the stage iii / iv distinction . manual scoring of stages i and iii was done in 30 sec increments . at that resolution , epochs manually labeled as stage iii could not be disambiguated from epochs manually labeled as stage ii or stage iv in the majority of recordings and epochs manually labeled as stage i could not be distinguished from epochs manually labeled as stage ii , rem or w in most recordings in the pfs . thus it is conceivable that stages i and iii are not stationary sleep states per se but rather are transitional . however rem was easily distinguishable from waking . thus , human rem sleep should no longer be thought of as “ awake - like ” or “ paradoxical ”. a k - means clustering algorithm ( fig2 ) was applied to the normalized data in the spaces above to classify sleep states . even though the va and mpi data were filtered differently , the general position of the sleep and waking clusters was similar across sets . moreover , although the algorithm was optimized on the mpi data set , it performed at 80 . 6 % on the va data , which is unprecedented using a single channel of data and is similar to the performance of other algorithms using many more channels . ( flexer , a ., et al ., artif intell med . 33 , 199 ( 2005 ). the standard error of the mean was also lower for the va set than the mpi set even though the former had 6 subjects and the latter had 20 subjects ( 1 . 73 % vs . 1 . 78 %, respectively ). the average agreement rate with human scoring on the full data set was 77 . 58 % on 4 stages . this striking concordance can be visualized by overlapping automated and manually derived hypnograms , which plot sleep stages for a given subject over a given night . in two out of twenty - six recordings , it appeared that the algorithm was mislabeling the data and in these cases . while that data appeared different when compared to the rest of the data set , visualization of the manual scoring on the preferred frequency map did however show separate signatures for sleep and waking stages . on the va data , when the algorithm &# 39 ; s performance was compared against data rescored by the same person or scored by a more experienced scorer , the average agreement rate with the algorithm increased and was in the 82 . 4 - 83 . 3 % range . further normalizations in time and frequency can be applied to the whole night spectrogram , at both a 30 sec ( fig7 a , c , fig . s 8 ) and a 1 sec resolution ( fig7 b , d , fig9 ). here sleep and waking stages tile the entire 1 - 100 hz spectrum with rem , w and is exhibiting broadband patterns ( fig8 , fig9 c - d ). in this space , one can measure the fragmentation in normalized power across time ( temporal fragmentation ) ( fig1 , 21 - 22 ). this analysis revealed a bimodal distribution for rem sleep . this pattern persisted when the frequency range was narrowed to 4 - 40 hz ( data not shown ). the more fragmented part of rem accounted for ( mean ± s . e . m ) 26 . 18 ± 1 . 7 % of rem at a rate of 37 . 42 ± 2 . 70 epochs per night lasting an average of 36 . 18 ± 1 . 27 seconds and separated by an average of 129 . 08 ± 11 . 04 seconds of stable rem ( fig2 ). these components of rem do not correspond to tonic and phasic rem ( fig2 ) and exhibit different spectral signatures ( fig2 ). this unstable part of rem sleep was more likely to be confused with stage ii than the stable part ( fig3 - 32 ). in these cases some spindles and k - complexes in the presence of rem caused these epochs to be scored as stage ii ( fig2 ) even though they would have been scored as rem at a finer temporal resolution . according to r - k rules , no spindles or k - complexes can be separated by less than 3 minutes in rem . while k - complexes and spindles can be found in rem , according to the analysis presented here , these signals are not responsible for the bimodal temporal fragmentation pattern observed in rem since manually scored rem , presumably devoid of spindles and k - complexes , still exhibits this pattern ( fig1 a - b , 21 - 22 , 31 right columns ). moreover , rem still exhibited a bimodal distribution on a spectrum without spindle frequency power . the temporal fragmentation is sensitive to sudden changes in normalized power . such changes can also be brought about by artifacts and the changes they produce will be enhanced in the background of a low power eeg . therefore , artifacts of some sort could be responsible for most if not all of the bimodal temporal fragmentation of rem . when epochs adjacent to epochs known to contain movement artifacts were discarded from the analysis as well as any epoch having a preferred frequency greater than 25 hz , the percentage of unstable rem epochs was diminished even if the bimodal pattern could still be seen . the bimodal pattern was even less apparent when more artifacts were isolated . however when these artifacts were included in the fragmentation analysis , in 4 out of 6 cases ( 5 out of 6 cases when rem was visually identified by a second scorer ), they accounted for a higher percentage of the non - fragmented portion of rem ( 6 out of 6 for automated scoring ) and in all but two cases for manual scoring ( non - fragmented portion of rem — 71 . 91 % in subject 9 and 50 . 73 % and 52 . 24 % in subject 20 , depending on the scorer ) and in all but one case for automated scoring ( non - fragmented portion of rem — 75 . 9 % in subject 9 ), they accounted for less than 50 % of either portion of rem . a nearest - neighbor analysis was performed on epochs which did not themselves include artifacts ( fig3 ). the fragmented portion of rem had almost in all cases more neighbors which contained an artifact than the non - fragmented portion , according to manual scoring ( 5 / 6 subjects for one scorer 6 / 6 subjects for the other ). when rem was detected automatically , in most subjects , the majority of both the fragmented and non - fragmented epochs were devoid of neighboring artifacts . further analysis of these data will be necessary to identify eeg features that might be responsible for the observed patterns and possibly a new state of sleep . nevertheless , temporal fragmentation provides yet another variable that easily distinguishes rem from both w and stage i ( fig3 ). in any of the technologies described herein , any variety of statistics can be generated from adjusted source data . for example , sleep statistics can be generated from adjusted source eeg data that has been classified into sleep states . exemplary sleep statistics can include information including sleep stage densities , number of sleep stage episodes , sleep stage average duration , cycle time , interval time between sleep stages , sleep stage separation statistics , onset of sleep , rapid eye movement sleep latency , regression coefficients of trends , measures of statistical significance of trends , and the like . in any of the examples herein , an electronic or paper - based report based on sleep state data can be presented . such reports can include customized sleep state information , sleep state statistics , pathological conditions , medication and / or chemical effects on sleep , and the like for a subject . recommendations for screening tests , behavioral changes , and the like can also be presented . although particular sleep data and low frequency information results are shown in some examples , other sleep data presenters and visualizations of data can be used . any of the computer - implemented methods described herein can be performed by software executed by software in an automated system ( for example , a computer system ). fully - automatic ( for example , without human intervention ) or semiautomatic operation ( for example , computer processing assisted by human intervention ) can be supported . user intervention may be desired in some cases , such as to adjust parameters or consider results . such software can be stored on one or more computer - readable media comprising computer - executable instructions for performing the described actions . such media can be tangible ( e . g ., physical ) media . the above has described how information can be used to determine sleep states . these techniques may also be used for other applications including characterizing sleep states , and other techniques . applications may include determination of whether a patient has taken certain kinds of drugs based on their sleep state , and based on variables that were previously determined as changing in brain function based on those sleep states . referring now to fig5 , which is a block diagram showing an exemplary system for determining a pathological condition of a subject from sleep states 500 . electroencephalography data for an animal is obtained and input into sleep state analyzer to determine a pathological condition of the subject . a pathological condition can be detected in an animal based on the sleep states 506 . for example , sleep states can be acquired for an animal 502 and analyzed 504 to determine whether the sleep states represent normal sleep or abnormal sleep . abnormal sleep could indicate a pathological condition 508 . for example , sleep states can be acquired from animals with pathological conditions and analyzed for common attributes to generate an exemplary distinctive “ pathological condition ” sleep state profile and / or sleep state statistics representative of having the pathological condition . such a profile or statistics can be compared to sleep states determined for an animal in order to detect whether the subject has the pathological condition or any early indicators of the pathological condition . any variety of pathological conditions can be detected and / or analyzed . for example , sleep related pathological conditions can include epilepsy , alzheimer &# 39 ; s disease , depression , brain trauma , insomnia , restless leg syndrome , and sleep apnea . for example , polysomnographically , subjects with alzheimer &# 39 ; s can show decreased rapid eye movement sleep in proportion to the extent of their dementia . narcolepsy is associated with sudden transitions into rem . it has recently been reported that there are instability patterns in the eeg of narcoleptic animals . if these apply to rem and humans as well , narcoleptics may have a marked difference in their rem fragmentation patterns as well . many other diseases have been linked to sleep disorders . for example , depression is associated with short rem latency and increased rem sleep . parkinson &# 39 ; s disease is also associated with rem behavior disorder . alzheimer &# 39 ; s patients already have unstable sleep patterns . these conditions and their treatment ( maois , used against depression block rem ; cholinesterase inhibitors , used against alzheimer &# 39 ; s disease , affect rem as well ) may be associated with new expressions of stable and unstable rem , which could be used to assess both pathology and treatment . the preferred frequency and iterated preferred frequency plots could also help to extract biomarkers of pathology and treatment . in any of the technologies described herein , the effect of medications and chemicals on sleep states of an animal can be determined via analyzing source data obtained for an animal . for example , sleep states can be modified by alcohol , nicotine , and cocaine use . exemplary medications that affect sleep include steroids , theophylline , decongestants , benzodiazepines , antidepressants , monoamine oxidase inhibitors ( e . g ., phenelzine and moclobemide ), selective serotonin reuptake inhibitors ( e . g ., fluoxetine ( distributed under the prozac ® name ) and sertralie ( distributed under the zoloft ® name ), thyroxine , oral contraceptive pills , antihypertensives , antihistamines , neuroleptics , amphetamines , barbiturates , anesthetics , and the like . sleep patterns may be used as a diagnostic as described above for pathological conditions and medication effects . the example below illustrates how sleep patterns may be used as a biomarker to identify individuals . sleep data for four pairs of twins were analyzed utilizing the exemplary sleep staging techniques described above . each column in 1 - 4 corresponds to 4 pairs of twins ( pair 1 is fraternal , pairs 2 - 4 is identical ). only rem is shown ( temporal fragmentation across time ). twins exhibit a similar temporal fragmentation pattern ( fig3 ). the general structure and techniques , and more specific embodiments which can be used to effect different ways of carrying out the more general goals are described herein . throughout this application , various publications , patents , and / or patent applications are referenced in order to more fully describe the state of the art to which this invention pertains . the disclosures of these publications , patents , and / or patent applications are herein incorporated by reference in their entireties , and for the subject matter for which they are specifically referenced in the same or a prior sentence , to the same extent as if each independent publication , patent , and / or patent application was specifically and individually indicated to be incorporated by reference . although only a few embodiments have been disclosed in detail above , other embodiments are possible and the inventors intend these to be encompassed within this specification . the specification describes specific examples to accomplish a more general goal that may be accomplished in another way . this disclosure is intended to be exemplary , and the claims are intended to cover any modification or alternative which might be predictable to a person having ordinary skill in the art . for example , other applications are possible , and other forms of discrimination functions and characterization is possible . while the above extensively described characterizing the frequency in terms of its “ preferred frequency ”, it should be understood that more rigorous characterization of the information may be possible . also , while the above only refers to determining sleep states from the eeg data , and refers to only a few different kinds of determination of sleep states , it should be understood that other applications are contemplated . having illustrated and described the principles of the invention in exemplary embodiments , it should be apparent to those skilled in the art that the described examples are illustrative embodiments and can be modified in arrangement and detail without departing from such principles . techniques from any of the examples can be incorporated into one or more of any of the other examples . also , the inventors intend that only those claims which use the words “ means for ” are intended to be interpreted under 35 usc 112 , sixth paragraph . moreover , no limitations from the specification are intended to be read into any claims , unless those limitations are expressly included in the claims .
traditional analysis of sleep patterns requires several channels of data . this analysis can be useful for customized analysis including assessing sleep quality , detecting pathological conditions , determining the effect of medication on sleep states and identifying biomarkers , and drug dosages or reactions .