as can be appreciated , osteopenia and osteoporosis , i . e ., the loss of bone mass , arises to a large degree from the natural aging process , but also to a lesser degree from a decrease in muscle activity , such as due to bed rest . this loss can be detected using vibrational stimulation arising from a source external to the body , which creates mechanical , frequency specific , low level oscillations in the subjacent bones . compared to imaging techniques such as dxa that neglect the volume of the measured bones , or ultrasonometer that measure localized speed of sound , acoustic vibrations around the natural frequency of an object excites the whole object , independent of size and shape . changes in shape and bone mass content will alter its responsive frequency , and can thus be detected by the disclosed device and / or using the disclosed methods . in more technical terms if one considers trabecular bone as a random network of struts , it is known from percolation theory that when a fraction vo of struts are randomly removed from a large network , it will spontaneously fragment into multiple segments . the strength of the network , as a whole , will vanish at that point . the value ( at which fragmentation occurs ) vo , known as the bond - percolation threshold , depends on the class of networks being evaluated . for three dimensional voronoi structures ( structures assuming a random network of struts ) its value is v 0 = 0 . 5 . thus , when half the struts on a voronoi network are randomly removed , its fracture load ( strength ) vanishes . because trabecular perforation has been shown to be the dominant cause of bone damage , it can be expected that the fracture load of trabecular bone will vanish when it loses about half its mass . consequently , the commonly used power - law relationship ( regression fit that correlates fracture load with bone mineral density ) between fracture load and density of bone ( bell , g . h ., o . dunbar , et al . ( 1967 ). “ variations in strength of vertebrae with age and their relation to osteoporosis .” calcified tissue research 1 ( 1 ): 75 - 86 ); mcelhaney , j . h ., j . l . fogle , et al . ( 1970 ). “ mechanical properties on cranial bone .” j biomech 3 ( 5 ): 495 - 511 ; carter , d . r . and w . c . hayes ( 1977 ). “ the compressive behavior of bone as a two - phase porous structure .” journal of bone and joint surgery 59 - a : 954 - 962 ; mcbroom , r . j ., w . c . hayes , et al . ( 1985 ). “ prediction of vertebral body compressive fracture using quantitative computed tomography .” j bone joint surg am 67 ( 8 ): 1206 - 14 ; rice , j . c ., s . c . cowin , et al . ( 1988 ). “ on the dependence of the elasticity and strength of cancellous bone on apparent density .” journal of biomechanics 21 ( 2 ): 155 - 168 )— all references listed above are hereby incorporated by reference .) can only be an approximation that is invalid for weak bones ; because under a power - law both the strength and density will vanish simultaneously . the power - law relationship does not take into account the spatial organization of the trabecular tissue . accordingly , an estimate of the percolation threshold v 0 for trabecular bone at different anatomic sites is a need that still exists in the art , because trabecular bone resembles disordered cubic networks more closely than voronoi networks . mosekilde , l . ( 1988 ). “ age - related changes in vertebral trabecular bone architecture — assessed by a new method .” bone 9 ( 4 ): 247 - 50 . in accordance with the present invention , a method for detecting osteopenia and osteoporosis comprises the step of inducing a relatively low level deformation in the bone tissue within a relatively large frequency range , preferably by mechanically loading the bone . the present invention includes both apparatus and methods for testing or characterizing tissue , particularly bone . certain embodiments also include using the apparatus described herein as a therapeutic tool , particularly in the area of chiropractics , pain relieve , bone growth stimulation , bone healing ( particularly after surgery ) and osteoporosis treatment at selected anatomic sites . in reconstructive surgery , requiring implants , the present invention can be utilized as an intra - operative and post - operative diagnostic tool to evaluate implant fixation and implant loosening over time . such embodiments generally comprise finding the natural frequency for a specified tissue in order to induce the maximum excitation of the tissue . other embodiments include business methods of using the devices and methods described herein as a means for generating a diagnostic database and generating a revenue stream based on subscriptions or access to diagnostic data . in general , the apparatus embodiments tissue quality assessment comprises a signal generator capable of generating a signal , a means for converting the signal into a vibration over a frequency range into one or more tissues at a designated point on a patient &# 39 ; s body , and a means for measuring a tissue response at that same point on the patient &# 39 ; s body . the device is capable of housing all of the necessary components within a single original device . the device should be capable of inducing a vibration over the frequency range of about 20 hz to about 1 . 5 mhz . the vibration generation may normally comprise a shaker having one or more devices selected from the group consisting of a piezoelectric device , an electromagnetic device , an electromechanical device , a pneumatic device and a hydraulic device . the means for measuring a tissue response comprises an impedance transducer , which may comprise one or more selected from the group consisting of an accelerometer , a displacement transducer and a load cell . further , the impedance transducer in certain embodiments are capable of measuring an input and output signal . other components are contemplated for various embodiments as would be understood by one of ordinary skill in the art such as power amplifiers , signal converters , data storage , power sources , data transmitters , input and output components , and the like . other embodiments comprise one or more additional portable and / or remote sensors for measuring localized tissue damage . each remote sensor should be capable of measuring the tissue response at a different region from the primary vibration point of origin using one or more accelerometer , load cell and / or a displacement transducer . one advantage to the present invention is that tissue can be evaluated without the need for stabilizers and / or the use of the 3 point testing procedure , which is typical in determining bone stiffness . referring now to fig1 , a tissue characterization device 100 is shown . tissue characterization device 100 includes a housing 110 , an actuator 120 , an accelerometer 130 , and a data acquisition system 160 . actuator 120 and accelerometer 130 are preferably located within housing 110 . data acquisition system 160 may be located within housing 110 or may be external to housing 110 . in some embodiments , tissue characterization device 110 further includes a tissue contact pressure sensor 150 located on the proximal end 112 of housing 110 . housing 110 may be constructed of any suitable material or materials . additionally , housing 110 may be a continuous piece ( not shown ) or shaped into a plurality of modules 114 . housing 110 preferably encloses the components of tissue characterization device 100 such that tissue characterization device 100 is a self - contained assembly . for example , in some embodiments , tissue characterization device 100 may additionally include a power means ( not shown ), such as a battery , so that tissue characterization device 100 is freely mobile , ie . cordless . actuator 120 preferably comprises a piezoelectric element 122 . when mechanical load is applied to a piezoelectric element , the piezoelectric element generates electric charge separation , resulting in electric field or voltage . conversely , when an electric field or voltage is applied to a piezoelectric element , the piezoelectric element undergoes a mechanical deformation . in some embodiments , piezoelectric element 122 consists of a plurality of thin layers ( not shown ) of electroactive ceramic material electrically connected in parallel . the total displacement d t piezoelectric element 122 moves is the sum of the displacements d i of the individual layers . the thickness t i of each individual layer determines the maximum operating voltage for actuator 120 . suitable piezoelectric materials include , for example , quartz , tourmaline , rochelle salt , barium titanate , and lead zirconate titanate ( pzt ). in some embodiments , piezoelectric element 122 may be replaced with an acoustic actuator ( now shown ). acoustic actuators may be desirable in some instances , such as those where a load - controlled device is required . data acquisition system 160 preferably includes a data logger 162 . in some embodiments , data logger 162 records data or measurements made by tissue characterization device 100 . in other embodiments , data logger 162 records and manipulates data . in still other embodiments , data logger 162 records and manipulates data and presents results . data logger 162 may transfer data , manipulated data , and / or results to a computer ( not shown ), e . g . via a wireless connection or cable . in a preferred embodiment , an internal control and / or measurement unit 180 applies a positive voltage (+) across piezoelectric element 122 . in response , piezoelectric element 122 mechanically deforms and forces a piston 132 and accelerometer 130 to move along x - axis 190 , as indicated by arrow 192 . in addition to piston 132 , accelerometer 130 additionally includes a load cell 134 connected to piston 132 . load cell 134 preferably comprises a weight element 135 . when the voltage across piezoelectric element 122 is kept at its current value , piezoelectric element 122 remains in its most recent position . in order to deform piezoelectric element 122 in the opposite direction , measurement unit 180 applies an inverse or negative voltage (−) across piezoelectric element 122 . this allows piston 132 to move back toward its original position along x - axis 190 , as indicated by arrow 194 . preferably , the voltage is an alternating voltage , such that piezoelectric element 122 is continuously changing between its original and deformed shapes . this in turn , causes piston 132 to move back and forth along x - axis 190 at a variable piston frequency and magnitude . applying an offset voltage can change the motion profile from oscillation to pulsation . as stated above , a pressure sensor 150 may be located on the proximal end 112 of housing 110 . pressure sensor 150 preferably comprises a compressible pad 152 such that upon contact with a patient &# 39 ; s body ( not shown ), pressure sensor 150 compresses accordingly . when a pressure sensor 150 is present , pressure sensor 150 is preferably in communication with internal control unit 180 . internal control unit 180 is capable of receiving pressure sensor data from pressure sensor 150 and indicating adequate pressure levels for taking repeatable measurements . when pressure sensor 150 is not present , tissue characterization device 100 may have a power control means 140 that is in communication with internal control unit 180 . in some embodiments , power control means 140 is a knob or wheel 142 such that at a first position ( not shown ) tissue characterization device 100 is “ off ” and when power wheel 142 is turned clockwise to a second position ( not shown ), tissue characterization device is “ on .” power control means 140 is preferably located at least partially external to housing 110 such that a user may physically access it . in a preferred embodiment , tissue characterization device 100 is configured such that the amount a user turns power wheel 142 counterclockwise corresponds to the voltage magnitude being sent to piezoelectric element 122 . it has been contemplated that internal control unit 180 has a number of predetermined settings ( i . e . frequencies and magnitudes ), which may correspond to visual marking on power wheel 142 . for example , power wheel 142 may have a “ 0 ” or “ off ” marking ( not shown ) to let the user know tissue characterization device 100 is off . power wheel 142 may additionally have markings , for example from 1 - 10 ( not shown ), to let the user know that tissue characterization device 100 is on and that there are different settings that tissue characterization device 100 can operate at ( i . e . different settings according to chosen anatomic sites — femur , tibia , spine , and others ). regardless of the mechanism for turning tissue characterization device 100 on , tissue characterization device 100 also preferably comprises an emergency off means 185 . in a preferred embodiment , emergency off means 185 comprises a feedback loop ( not shown ) within or in communication with internal control unit 180 . the feedback loop preferably monitors the amount of load placed on actuator 120 , and if the load / pressure exceeds a predetermined threshold , the feedback loop forces or instructs internal control unit 180 to stop sending voltage to piezoelectric element 122 . in other embodiments , emergency off means 185 comprises a button 188 located at least partially external to housing 110 so that a user may easily access it . similar to feedback loop 186 , button 188 is in communication with internal control unit 180 and forces or instructs internal control unit 180 to stop sending voltage to piezoelectric element 122 when button 188 is pushed . as described above , piston 132 moves back and forth along x - axis 190 at a variable piston frequency and magnitude . the piston frequency and magnitude directly correspond to the voltage frequency and magnitude applied to piezoelectric element 122 . as can be appreciated , piston 132 is the portion of tissue characterization device 100 that contacts a patient &# 39 ; s body and applies a mechanical deformation to an anatomical region of the patient &# 39 ; s body . therefore , the piston frequency and magnitude also directly correspond to the patient contact frequency and magnitude . the frequency range of the mechanical stimulus applied to the patient &# 39 ; s body or tissue is generally between about 1 hertz and 1 megahertz , preferably between about 20 hertz and about 600 kilohertz and is more preferably between about 100 and about 600 hertz . the magnitude of the deformation induced in the tissue is preferably between about 0 . 1 and about 5 millimeters and is more preferably between about 0 . 5 and about 2 millimeters at the above - mentioned frequency ranges . the optimal frequency of the mechanical deformation is between about 100 hertz and about 600 hertz , and the optimal peak - to - peak level of the deformation induced in the tissue is about 0 . 5 millimeter at the optimal frequency range . the optimal contact load applied to the patient through tissue characterization device 100 and monitored with pressure sensor 150 is in the range of about 10 newtons to about 20 newtons . however , the contact load need only be maintained at a range between 10 newtons and 100 newtons to result in a valid test . in a preferred embodiment , a stepwise frequency “ sweep ” is applied to the patient from about 20 hz to about 1 . 5 mhz , preferably about 20 hz to about 600 khz . in a preferred embodiment , piston 132 or pressure sensor 150 of tissue characterization device 100 is placed in contact with a patient &# 39 ; s body at a particular location , e . g . spine ( fig2 ). upon contact with the patient &# 39 ; s body , pressure sensor 150 sends a signal to internal control unit 180 , which turns tissue characterization device 100 on . in the on mode , internal control unit 180 applies a voltage to piezoelectric element 122 at a predetermined frequency and magnitude . this causes piston 132 to oscillate back and forth , generating a mechanical deformation on the patient &# 39 ; s local anatomical tissue at a specific frequency . data logger 162 then records the body &# 39 ; s response to the mechanical deformation by measuring the output signal from the load cell and accelerometer as a function of frequency . the data logger is an off - the - self device driven by a microprocessor . a internal algorithm transforms the data preferably using fast - fourier - transformation and determines the peak response of the tissue versus frequency . these preliminary data will be used to determine the validity of the test . indicated on the display 311 . in some embodiments , tissue characterization device 100 is used to characterize hard tissues including bone and joint tissues . bone and joint tissues generate different responses from each other when a mechanical stimulus is applied to them . the present inventor has discovered that by applying a frequency sweep to the hard and soft tissues , the different responses can be separated from the complex tissue response . in a preferred embodiment , the frequency spectrum indicates the damping factor of the tissues . this enables a user to characterize multiple tissues , such as articular cartilage and the surrounding bone simultaneously . it has been contemplated that in order to determine changes within tissues such as bone and joint tissues , data sets could be used . these data sets may be established in a laboratory , i . e . in vitro testing and simulation of different injury mechanism , or be a previously recorded data set of a patient . when the data set comprises a previously recorded data set of a patient , the progress of a patient may be monitored . this is especially useful in osteoporotic patients , for which imaging techniques such as computerized tomography ( ct ) do not have the necessary resolution and are too expensive to repeat at every physical exam . the present invention overcomes the inherent disadvantages of conventional bone quality assessment tools . it is less dangerous and more easily applied to the patient , especially the elderly . the low level of the mechanical deformation minimizes the chance of injury to the patient and the bone tissue being characterized , and the relatively high frequency range of loading significantly reduces the period of time required for such mechanical exposure . finally , it would appear that bone and joint tissue is acutely responsive to stimuli induced in the above - mentioned frequency ranges , thereby having a therapeutic effect . in addition to characterizing hard tissues and joints , it has been contemplated that the present invention may be used to characterize other tissues including skin and lymphoid tissue . for example , the present invention may be used to detect abnormalities in both of these tissues such as melanomas . the same technique described above to characterize hard or soft tissue non - invasively can alternately be used to evaluate the stability of foreign implants in the body , such as fracture fixation devices , as well as monitor wound healing . in characterizing the dynamic response of patients implants in load - bearing anatomic sites the present invention can be utilized to monitor implant longevity . furthermore , during surgery , the pre - load generally applied on bone - implant constructs may be monitored and used as an indicator for adequate adjustment of the implant . additionally , during minimally invasive procedures such as a vertebroplasty or kyphoplasty , tissue characterization device 100 may be used to monitor the effectiveness of the treatment during operation before the patient applies his own weight onto the bone . corrections of the surgical approach while the patient is still in the operating room are feasible using the present device / technique . the devices described herein are preferably portable , inexpensive and easy to use . it is envisioned that the devices would be available to individuals as well as professionals . professionals could use the instant invention as quick diagnostic tools for detecting bone integrity , e . g ., fractures , bone density , etc ., in an ordinary examine room without the need for sending the patient to a separate room or facility ( as with certain existing technology ). the objective for individuals could be monitoring bone quality at home , thereby avoiding expensive doctor visits . fig3 shows another embodiment comprising of a portable device 200 comprising a piston 210 , an acceleration means 220 , an electromagnet 230 , an optional accelerometer 240 , a load cell 250 , a contact tip 260 , circuitry 270 , an energy source 280 and an external housing 290 . the operating piston 210 could be made of any suitable material such as plastic or metal , preferably metal . the basic components are known in the art and readily understood by those of ordinary skill . the contact tip 260 may be made of any suitable material such as plastic or metal , preferably plastic . the contact tip 260 and piston 210 may be integrally connected such that the movement of the piston 210 necessarily moves the contact tip 260 . alternatively , the piston 210 may not be integrally connected to the contact tip 260 . in other words , piston 210 would move independent of contact tip 260 through a through bore in the contact tip 260 such that piston 210 could directly contact the designated area being tested . in either embodiment , the contact tip 260 is in communication with load cell 250 so that a predetermined amount of contact force may be applied to a designated area for testing . load cell 250 may be used to determine if an appropriate amount of force is being applied to the target area prior to initiated the test procedure . if piston 210 and contact tip 260 are connected , then piston 210 must overcome the force being applied to the target area via the contact tip 260 for operation and data collection . if piston 210 moves independent of contact tip 260 then piston 210 would not have to overcome the pre - applied force at the target area . one skilled in the art will appreciate that such a design alternative may be a consideration in selecting the appropriate mechanism for moving piston 210 during operation . the acceleration means 220 should be able to accelerate piston 210 over a frequency range of about 20 hz to about 600 khz . circuitry 270 may comprise a processor and data storage means ( not shown ). the processor , preferably a microprocessor , may comprise any known in the art and should be capable of performing the functions described herein . the data storage means may comprise any means capable of storing data collected during operation including the parameters of the testing device for a given data sample collected . also the circuitry must be capable of allowing the stored data to be transferred to a remote device , preferably a stand alone computer or similar device . accordingly , portable device 200 may comprise a means ( wireless or otherwise ) for transferring data ( either live or stored ) to a remote device such as a computer . suitable non - wireless alternatives may comprise , but are not limited to , direct wiring to the remote device , ports that allow for connective wiring to the remote device and / or ports that allow communication through a cradle device ( wherein the cradle device is in communication with the remote device . circuitry 270 may also comprise external components such as a display means and one or more input means . for example , one embodiment may comprise one or more of the following : ( a ) an external toggle switch or button for powering on and off the device , ( b ) one or more entry buttons for selecting operating parameters , e . g ., anatomic area ( arm , hip , spine , etc ) to be tested , male vs female , age , etc ., ( c ) one or more led &# 39 ; s , lights , speakers or other audio device to indicate readiness , that the preload has been achieved , confirmation / rejection of data collection , battery status , etc ., and / or ( d ) a display panel , audio means or other means for displaying or indicating a message , such as a number or text message . energy source 280 may comprise a portable energy device , such as one or more batteries , or a means for transferring electricity to the circuitry 270 , such as electrical connections and necessary devices for making a standard connection with an electrical outlet . for embodiments using portable energy , the energy source 280 may include a receptacle for receiving the one or more batteries and optionally a means for recharging said batteries if applicable . all such simple electrical circuitry is readily known in the art . fig4 shows another embodiment comprising a device 300 including contact tip 301 , piston 302 , impedance head 303 , vibration generator 304 , vibration insulator 305 , vibration generator housing 306 , outer housing 307 , preload cell 308 , hinge 309 ( optional ), manual inputs 310 , display 311 , indicator outputs 312 , power amplifier 313 , charge amplifier for the receiving signal 314 , circuitry 315 , power pack ( energy source ) 316 , data transmitter 317 , data storage 318 and various electrical connectors 319 as needed . all housing and / or non - electrical type components , e . g ., contact tip 301 , piston 302 , vibration generator housing 306 , outer housing 307 , hinge 309 , etc ., may comprise any suitable material such as plastic or metal , preferably plastic . contact tip 301 may comprise any shape , e . g ., flat or concave , and may be easily exchangeable with various shaped tips designed for particular anatomical regions . piston 302 is preferably a threaded rod . impedance head 303 preferably comprises a means for measuring both the vibration ( frequency sweep ) produced by device 300 to excite the tissue being tested and the tissue response to the excitation . vibration generator 304 preferably comprises a shaker , wherein the shaker induces vibration in the tissue or excites the tissue being tested using one or more devices selected from the group consisting of a piezoelectrical device , an electromagnetic device , an electromechanical device , a pneumatic device and a hydraulic device . vibration insulator 305 and vibration generator housing 306 may alternatively comprise a single entity that is able to perform both functions . vibration insulator 305 dampens the generated vibration or frequency sweep to improve the integrity of the signal and prevent transfer to the outer housing 307 and user hand . preload cell 308 measures the deformation of the outer housing 307 , thereby indicating contact force between device and testing object . manual inputs 310 may comprise one or more buttons or keypads for powering device 300 and / or selecting operational parameters . typical operational parameters may include anatomical site to be tested , age of patient , sex of patient , body weight , shoe size , file name for test data , and the like . display 311 may comprise any display means known in the art and is not critical to the invention . for example , display 311 may comprise an lcd screen . likewise , indicator outputs 312 may comprise any technology known in the art such as led technology . indicator outputs 312 may flash or otherwise indicate when the appropriate load is being applied for a test run , when a valid set of data has been collected , and the like . inputs 310 , display 311 and outputs 312 may alternatively be combined as technology allows and comprise one or more of the following : ( a ) an external toggle switch or button for powering on and off the device , ( b ) one or more entry buttons for selecting operating parameters , e . g ., anatomic area ( arm , hip , spine , etc ) to be tested , sex of patient , age , etc ., ( c ) one or more led &# 39 ; s , lights , speakers or other audio device to indicate readiness , that the preload has been achieved , confirmation / rejection of data collection , battery status , etc ., and / or ( d ) a display panel , audio means or other means for displaying or indicating a message , such as a number , text or graphic type message . power amplifier 313 increases the signal from the signal generator ( not shown ) contained in circuitry 315 and drives the vibration generator 304 . charge amplifier 314 converts the charge output from a piezoelectric , capacitive or other charge - producing sensor to a signal such as analog voltage or current . as with all embodiments , circuitry 315 may comprise a processor and data storage means . as shown , data storage means may comprise a separate storage device 318 . the processor , preferably a microprocessor , may comprise any known in the art and should be capable of performing the functions described herein . the data storage means ( whether part of circuitry 315 or as a separate component 318 may comprise any means capable of storing data collected during operation including the parameters of the testing device for a given data sample collected . also the circuitry must be capable of allowing the stored data to be transferred to a remote device , preferably a stand alone computer or similar device . accordingly , circuitry 315 or component 317 may comprise a transmitter ( wireless or otherwise ) for transferring data ( either live or stored ) to a remote device such as a computer . suitable non - wireless alternatives may comprise , but are not limited to , direct wiring to the remote device , ports that allow for connective wiring to the remote device and / or ports that allow communication through a cradle device ( wherein a cradle device is in communication with the remote device . as stated previously , power pack 316 may comprise a portable energy component , such as one or more batteries , or a means for transferring electricity to the circuitry 315 , such as electrical connections and / or necessary devices for making a standard connection with an electrical outlet e . g ., cables or cords . for embodiments using portable energy , the energy source 316 may include a receptacle for receiving the one or more batteries and optionally a means for recharging said batteries if applicable . all such simple electrical circuitry is readily known in the art . fig5 shows a portable extension 400 comprising a sensor 401 and circuitry 402 connected by wire 407 . for all embodiments of the present invention , wireless connections ( where technologically available ) are interchangeable for hardwired connections . circuitry 402 comprises a housing or platform 403 , processor 404 and transmitter 405 . transmitter 405 may be integrated as part of processor 404 ( not shown ). alternatively , sensor 401 and circuitry 402 may be integrated into a single structure . other embodiments of the present invention comprise multiple sensors each with its own corresponding circuitry or , alternatively , multiple sensors connected to a single portable circuitry . regardless , the processor 404 , transmitter 405 and energy source 406 comprise similar components as described in connection with other embodiments described herein . the data collected by the portable extension should be capable of transmitting via transmitter 405 or processor 404 said data from the portable extension 400 to either the primary device ( such as device 300 described in fig4 ) or to a remote device such as a computer . fig6 shows a schematic of one embodiment for a portable device 500 in accordance with of the general principles of the present invention . it is worth noting that the invention is not limited to the exact representation in fig6 , but instead fig6 is provided to give greater clarity by way of a specific example . one of ordinary skill will understand that many modifications to such an exact schematic are possible without leaving the spirit and scope of the invention . for example , charge amplifiers and / or buffers are easily moved within the device from one unit to another or combined to a single circuitry . accordingly , portable device 500 comprises a processor unit 501 and mechanical unit 502 . processor unit 501 comprises signal generator 503 , digital to analog converter 504 , a plurality of buffers 505 , 521 , 522 , 523 , 524 , a plurality of analog to digital converters 506 , 525 , 526 , voltage to signal converter 507 , a processor 508 , and display 509 . mechanical unit 502 comprises power amplifier 510 , shaker 511 , preload load cell 512 , trigger 513 , indicator 514 , and impedance transducer 515 . impedance transducer 515 comprises accelerometer 516 and load cell 517 . mechanical unit 502 may also comprise charge amplifiers 518 , 519 to correspond to accelerometer 516 and load cell 517 , respectively . optionally , device 500 comprises one or more stand alone output device 520 such as a monitor , lcd screen , plotter , printer , storage unit , and the like . processor unit 501 communicates with mechanical unit 502 and output device ( s ) 520 using either hardwired or wireless connections . another embodiment of the present invention comprises a method for collecting data regarding tissue characteristics using one or more of the device embodiments or obvious variations thereof described herein . collecting the data sample may be accomplished by applying a load to the target area , i . e ., pressing the contact tip of the device to the skin , bone , or implant until an indicator on the device displays that sufficient force has been applied . once the load is applied , the shaker would be activated , either manually or automatically , to cause the piston to accelerate at a frequency range of about 20 hz to about 600 khz . an accelerometer and / or other device may then measure the resulting force and / or acceleration of the piston and the data would be stored in the circuitry storage means . once data has been collected , an output signal can confirm the sample collected appears useable and a new sample can be collected . one or more samples may be collected and 3 or more are preferable . the collected data samples can then be evaluated . evaluation may take place within the device itself depending upon the sophistication of the software / hardware . alternatively , the device may transfer the data to a computer or other device for evaluation . the data may be evaluated against previous measurements taken by the patient or with a database of known measurements . in general , the data should be manipulated using an algorithm to analyze the data and produce a value . it is preferred that evaluation be carried out by using the calculated value and comparing it with other values calculated in a similar manner . once evaluated an output may be generated such as a simple yes / no , a percentage value of bone integrity , an extrapolated value for bone integrity and / or a normalized value having significance with the patient . in one embodiment , the database values would be stored in the remote computer . periodic updates may be downloaded to enlarge the number of values available for comparison . alternatively , the database values may be stored on a second remote device accessible via the internet , such that evaluation may comprise sending the collected data and / or calculated values via the internet to the second remote device for evaluation and / or comparison . in these embodiments , the second remote device would reply with the output signal described above . fig6 schematic may also be used to illustrate , in general , one embodiment for a testing method of the present invention . a load is placed on an anatomic region of a patient , preferably a location that has a thin layer of soft tissue over a bone . the preload load cell 512 measures the load and activates a trigger 513 that in turn can activate an indicator 514 so that the user may identify when the appropriate load is achieved for testing . the load range should be in the range of about 10 newtons ( n ) to about 100 n , preferably from about 10 n to about 20 n . alternatively , trigger 513 may automatically activate the testing sequence by activating signal generator 503 and optionally initialize software . signal generator 503 produces a signal capable of activating shaker 511 . as shown , the signal is a digital voltage signal that is converted via the digital to analog converter 504 . after passing through buffer 521 , the signal is amplified ( if necessary ) via power amplifier 510 , which in turn powers shaker 511 . shaker 511 generates a vibration or frequency sweep that can be transmitted into the tissue at the point under load via a contact tip 301 in connection with shaker 511 , preferably by means of a rod and casing ( not shown ). the contact tip is also operative connected to an impedance transducer 515 having components 517 , 516 for measuring the signal going out ( i . e ., signal being transferred into the tissue ) and the responsive signal , respectively . as shown , component 515 comprises an accelerometer 516 and comprises a load cell 517 . the signal generator 503 , power amplifier 510 and shaker 511 produce a vibration or frequency sweep over a specified range based on a given set of operational parameters such as anatomic site , age and sex of the patient . one or more of these operational parameters may be inputted by the user prior to , during or after the preload load is executed , preferably prior to placing the region under load . low and high frequencies are utilized . a preferred frequency sweep range is from about 20 hz to about 1 . 5 mhz . the preferred low frequency is from about 20 hz to about 600 khz and high frequency is from about 800 hz to about 1 . 5 mhz . the input and output signals are measured via load cell 517 and accelerometer 516 , respectively , and transmitted to the voltage to signal converter 507 . the signals may be amplified using charge amplifiers 518 ( responsive signal ), 519 ( outgoing signal ), optionally buffered via buffers 524 , 523 and converted in analog to digital converters 526 ( responsive signal ), 525 ( outgoing signal ) prior to reaching the voltage to signal converter 507 . the converted signal may then be processed by processor 508 to produce a user response . a user response may be communicated via an output device , such as a display panel , audio means or other means for displaying or indicating a message , such as a number , text or graphic type message . any output device known in the art is sufficient . in one embodiment , generating a user response may comprise processor 508 calculating a relative value based on comparative data from a database . alternatively , generating a user response may comprise processor 508 calculating a relative value based on an analytical method . in yet another embodiment , the user response may comprise a diagnostic value for osteoporosis or bone fractures . other embodiments may include one or more of the basic steps of ( a ) selecting an anatomical region for testing ; ( b ) selecting operational parameters ; ( c ) apply a force to the anatomical region to achieve constant contact pressure with the region ; ( d ) inducing a frequency sweep at the anatomical region while maintaining contact pressure ; ( e ) collecting tissue response data ; and ( f ) providing a response based on the data from step ( e ). step ( f ) may comprise one or more of the following steps ( i ) confirming that the data collected in step ( e ) was a valid sample ; ( ii ) calculating one or more raw values based on the data collected in step ( e ); ( iii ) comparing the raw values from ( ii ) to either : ( 1 ) one or more previous calculated raw values for a given patient at the same anatomic site , ( 2 ) one or more previous calculated raw values for a given patient of his / her contralateral side ( if available ), or ( 3 ) a database of similar values ; ( iv ) determining a relative diagnosis based on the comparison in step ( iii ); and ( v ) communicating the relative diagnosis to a user . further , the response of step ( f ) may comprise one or more forms selected from the group consisting of alpha - numeric characters , electronic data , graphs , pictures , video , charts , icons or sounds . fig7 shows a flow chart 600 illustrating a method for data collection in accordance with one embodiment of the present invention . the user or some other trigger initiates 601 the start measurement sequence . the user is allowed to choose a new test sequence or open an existing set of data 602 . if the use does not choose a new test , she can select from a one or more old data sets and optionally other choices 603 . opening an old set of data prompts additional options to manipulate , transfer , save or otherwise use the data 611 . if the new test option is selected at 602 then the user is prompted to input a file name 604 . the user can then select operational parameters 605 as have been described herein . the testing or data generation is then conducted by initiating an auto sequence program 606 . once collected the data is saved to temporary or permanent storage 607 and the user is prompted to select whether or not more tests are to be run 608 . if more tests are selected the user is then taken back 609 to action 604 and the sequence repeats . if no additional tests are selected the user is prompted for additional options to manipulate , transfer , save or otherwise use the data collected . fig8 shows a flow chart 650 illustrating a method for the testing sequence in accordance with one embodiment the present invention . once the auto sequence program is initiated 612 a processor reads the operating parameters 613 ( comprising at least one or more user inputs ). the system is checked and indicates active sensors 615 . at this point the user should be applying a force against the tissue by pressing the tip against the anatomic region to be tested . a load cell override 616 may be selected 614 . if selected the testing sequence will be triggered 621 . alternatively , the user may depress a manual trigger 623 , which would trigger the testing sequence 621 . if no load cell override command is activated , then the static force measured within the housing will be evaluated 617 to determine if the force is too low 618 or too high 620 and result in a too low 619 or too high 622 contact force alert . static force that maintains within a specified range will allow the release trigger to start the testing sequence 621 . once initiated , a frequency sweep will start 624 over a specified range . contact force is monitored 625 throughout the frequency sweep duration . if the contact force falls below the minimum 626 the sweep is terminated 629 and a low - force indicator is activated 619 . if the contact force becomes too great 627 the sweep is terminated 630 and a low - force indicator is activated 622 . a frequency sweep that is achieved without a contact force violation continues over the entire specified range 628 . the responsive signal from the tissue is collected as sensor data 631 . fig9 shows a flow chart 700 illustrating a method for data analysis in accordance with one embodiment of the present invention . the user may select from options 611 . option may comprise , among other things , obtaining a raw response value or a comparative value . obtaining a raw response value comprises the following steps : ( 1 ) data is loaded / read into a processor 701 ; ( 2 ) the user selects an analysis algorithm 702 ( alternatively , this step may be omitted if an algorithm has been pre - selected ); ( 3 ) constants are initialized 703 ; ( 4 ) an optional data reduction or filter process is carried out 704 ; ( 5 ) a value is calculated according to an analytical method 705 ; ( 6 ) a parameter convergence test is carried out 706 to determine if the value ( s ) are in an expected range 707 ( if out of range the sequence is restarted at step ( 3 )); converged data is used to either ( a ) calculate a user response value 708 , or ( b ) produce a graphic response 709 ; ( 7 ) optionally create a prognosis based on comparing ( 6 )( a ) or ( 6 )( b ) with data from a database or look - up table ; ( 8 ) prompt the user to select another analysis 611 or to end the analysis sequence 711 . selecting to end the analysis sequence at 711 will end the analysis sequence 712 . obtaining a comparative value comprises the following steps : ( 1 ) data is loaded / read into a processor 713 ; ( 2 ) the user selects an analysis algorithm 714 ( alternatively , this step may be omitted if an algorithm has been pre - selected ); ( 3 ) constants are initialized 715 ; ( 4 ) an optional data reduction or filter process is carried out 716 ; ( 5 ) a value is calculated according to an analytical method 717 ; ( 6 ) a parameter convergence test is carried out 718 to determine if the value ( s ) are in an expected range 719 ( if out of range the sequence is restarted at step ( 3 )); converged data is compared with to a previous data set from the same patient to calculate ( a ) a comparative response value at 708 , or ( b ) produce a graphic response 709 ; ( 7 ) optionally create a prognosis based on comparing ( 6 )( a ) or ( 6 )( b ) with data from a database or look - up table ; ( 8 ) prompt the user to select another analysis 611 or to end the analysis sequence 711 . selecting to end the analysis sequence at 711 will end the analysis sequence 712 . differences measured between various data sets obtained ( i . e ., transfer function results obtained at different time points ) for the same subject are indicative of changes within the soft and / or hard tissue . in a similar manner , comparison to a standard database can be made , or on the contralateral side if applicable . certain embodiments comprise the use of an analytical method . the following are basic analytical methods that are within the spirit and scope of one or more embodiments of the present invention . the methods are not intended to be limiting , but are provided by way of example for clarity . suitable analytical methods comprise calculating a ratio between a high frequency response of the tissue being tested and a low frequency response . based on preliminary numerical analysis , the high frequency response can be used to determine the amount of trabeculae ( bone tissue ) present , while the low frequency response indicates the amount of those trabeculae being utilized for load transfer . a fraction close to one indicates healthy bone , whereas a fraction close to 0 . 5 is indicative of a very fragile bone . the low frequency response is measured at the first natural frequency mode , while the high frequency is in the upper kilohertz range . it has been known that the inefficiency of load transmission is the principle cause for reduction of fracture load . this inefficiency is dependent on tissue level damage , which can be estimated by the fraction of struts belonging to a term called stress backbone . the required fraction can be obtained if the number of trabeculae ( struts ) that belong to the stress backbone and the total number of trabeculae on a sample of bone are known . the following observations provide a method to estimate these last two quantities . first , the elastic modulus e ( 0 ) of a sample depends on the number of stress pathways . each pathway allows additional load to be transmitted , and hence increases e ( 0 ). second , the response e ( ω ) of the sample to ( sonic or ultrasonic ) vibrations of a sufficient high frequency ω can be used to estimate the number of struts on the network . under these conditions , signals are attenuated very quickly inside the sample , and stresses are limited to the immediate neighborhood of the surface that is driven . as a result , the presence of large fractures in the interior of the network , which reduce the extend of the stress backbone , play no role in the response to high - frequency driving ; i . e ., all available trabecular elements near the surface contribute to e ( ω ). hence , assuming that the trabecular perforation is ( fairly ) uniform within each layer of trabecular elements , e ( ω ) can be used to estimate the total number of trabecular elements in the sample . thus , the ratio γ ≡ e ( 0 )/ e ( ω ), which can be obtained from vibrational analysis , provides an estimate of the fraction of struts belonging to the stress backbone . if , as we assume , the reduction in fracture load is related to the extend of the stress backbone , then the former is related to γ . stated differently , this analytical method comprises calculating a ratio of low frequency response over high frequency response of a tissue sample . the low frequencies excite the trabeculae or struts throughout the entire bone length , i . e ., low frequencies are able to penetrate the whole tissue being tested . this identifies the integrity of the sample or more correctly the amount of tissue that is still functioning properly , i . e ., the mass of tissue that can still bear a load . the high frequencies do not excite or penetrate the entire tissue being tested . instead , higher frequencies penetrate only a thin layer of the tissue , but identify a total mass across in that plane . the total mass can then be projected out for the tissue being tested . it is assumed that a nearly uniform mass distribution exists in the tissue . accordingly , the ratio of low frequency over high frequency responses gives a percentage of the amount of active “ load bearing ” tissue . another embodiment comprises using an analytical method comprising a step - wise frequency sweep being applied to the object being tested . the frequency sweep starts at a low frequency , preferably 20 hz , and stops at a high frequency , generally about 2000 hz ( preferably 100 - 1000 hz , with a predetermined step function . the data recorded from an accelerometer , load cell , and driving signal in the time domain are converted to a frequency domain using a discrete fourier transformation at each driving frequency . various filters ( known and understood by one of ordinary skill in the art ) are used for data smoothing and reduction of noise . the data are normalized to the maximum peak measured at any frequency . a two - dimensional plot is generated with one axis being the driving frequency ( step - wise frequency ), and the other axis the frequency spectrum as analyzed using the fourier transformation . the amplitude may be color coded to visualize variations in signal intensity . a log or linear color scheme can be used , preferably log scale . a pattern analysis approach is used to identify regions that change according to the pathology , for example osteoporosis . different regions of the plot represent responses of various tissues . the low frequency response represents the interaction between soft and hard tissue , an intermediate region includes the resonance of bone and information of its mass , while the high frequency region contains information about the soft tissue properties . therefore , variations in patterns obtained from the aforementioned 2d plot can provide information on bone quality , apparent damage , and joint properties . in addition , a database of plots from known tissue samples ( patients with known problems ) can provide comparative results . further , it is believed that the plots may also be used to separate joint problems ( i . e ., herniated disc or facet joint pathologies ) from bone tissue problems . in an alternative approach , the fourier transformed signals described in the methods above may be used to calculate other mechanical transfer functions such as accelerance , effective mass , mobility , impedance , compliance , and stiffness in addition to cumulative energy . table 1 below lists the various functions that may be calculated from the data and fourier transformed signals generated in accordance with the various embodiments of the present invention . a similar pattern recognition method as identified above is utilized to identify changes between measurements and to provide an indication for the progression of the pathology or healing / treatment . mechanical vibration transfer functions definition accelerance acceleration / force effective mass force / acceleration mobility velocity / force impedance force / velocity compliance deformation / force stiffness force / deformation an alternative approach to method 2 is to utilize fast fourier analysis of the recorded signal for a particular excitation frequency range , preferably at a width of 200 hz . multiple analysis at various frequency intervals may be conducted . for each frequency range , one or more of accelerance , effective mass , mobility , impedance , compliance , stiffness , and cumulative energy may be determined , as they represent mechanical vibration transfer functions derived from various excitation ( input ) and response ( output ) signals . variations in these signals / spectra for the different frequency range s are indicative of damage to the different tissue types , such as hard or soft tissue . a variation of method 4 comprises analyzing only single frequency steps . a discrete fourier transformation can be utilized to convert the time dependent signal into a frequency domain . mechanical transfer functions and cumulative energy may be calculated . peaks of the different modes of the mechanical transfer functions can be plotted cumulatively in one or more graphs . shape changes , ratios between first and consecutive peaks , slopes , and more may be calculated from the graph . variations in these signals / spectra for the different input frequencies selected are indicative of damage to the different tissue types , such as hard or soft tissue . a log or linear frequency sweep is used to test a patient &# 39 ; s bone . after fast fourier analysis of the output signal ( accelerometer ) a spectral analysis is performed on the low frequency range to determine bone mass . based on known anthropometric data , bone size can be estimated from the calculated bone mass . the natural frequency is determined from the peak analysis of the frequency spectra . a well established correlation between bone size times natural frequency and age is used to calculate a parameter that can separate normal aging patients with patients suffering from excessive bone loss or bone micro - damage . alternatively , the value for bone mass is multiplied with the natural frequency directly and used as a predictor of pathology . this technique differs from previous techniques in that it uses much more information from the vibrational behavior of tissue . the models include mass elements , springs , and dashpots , representing the various tissues and joints . a simplified viscoelastic model will be used to fit the experimental data after curve fitting and regression analysis . this model may contain only three elements ( mass , spring , dashpot ) to describe a differential equation for acceleration , velocity and displacement behavior over various frequencies . utilizing a least - square method , parameter fitting will be conducted to estimate values for the three elements ( mass , spring , dashpot ). the viscoelastic parameters obtained , will be used as initial guess values for a more complex viscoelastic model , which will include the flesh , bone , and the joints . a similar least - square method will be utilized to obtain values for the differential equation representing the more complex viscoelastic model . the viscoelastic model will assist in separating the dynamic response measured with the accelerometer and load cell for the different tissue components . the parameter values are directly related to biomechanical properties of the tissue . pathological disease and sources for pain can easily be identified using this approach . in predicting fracture risk , one has to know the strength of a material and the load it is subjected to during regular activities . the following technique calculates the ratio of predicted bone strength and normal postural load experienced by human vertebral bodies within an anatomic region . subjects are tested in an unloaded state and subsequently in a loaded state ; for the spine , this would correspond to a laid down position and standing up , respectively . based on the vibration response of the tissue and the data obtained through this methodology an estimate of the physical load applied to the bone tissue in vivo is feasible . in applying principles mentioned in method 7 , material properties such as stiffness of the bone tissue can be determined for the unloaded and loaded case . the difference in vibrational reading between the two cases is to the largest extend due to the loading of the tissue in moving from they laid down position to an upright position . this difference can be expressed as deformation of the tissue due to loading , which in tern is a measure of stiffness . using the stiffness − strength correlation , a fracture risk prediction can be performed . again , a viscoelastic model may be used to fit the recorded data and deliver the physical quantities to describe the loading and tissue behavior . currently , there is no known method to assess loading conditions in vivo . published data are back calculations based on body mass distribution or inverse kinematics analysis . a similar procedure may be developed for the femoral neck . another embodiment of the present invention comprises using the apparatus described herein as a therapeutic tool , particularly in the area of chiropractics , pain relieve , bone growth stimulation , bone healing ( particularly after surgery ) and osteoporosis treatment at selected anatomic sites . such embodiments generally comprise finding the natural frequency for a specified tissue in order to induce the maximum excitation of the tissue . based on current knowledge known to the field , the human bodies capability to regenerate tissue is based on dynamic mechanical loading , this is particularly the case of hard tissue like bone and cartilage . cartilage is an avascular tissue , which through compression and relaxation transports nutrition and removes waste products . in bone , sportive activities such as bicycling or swimming have shown to effect bone tissue to a minor extend , while high impact loading such as jogging can stabilize or even increase bone mass . the exact mechanism behind the dynamic loading response is not known yet . the objective of our device is to regionally stimulate tissue through vibration , thereby promote and increase the regenerative and healing response . in applying a vibrational stimulus to tissue at the resonance frequency , the maximal amount of energy is transmitted to the tissue . one embodiment of the invention comprises a business method wherein the techniques described above are carried out and a fee is charged either for the updates to additional data and / or for access to the second remote computer having a dynamic database via the internet . the fee may include ( 1 ) a one - time fee , ( 2 ) an upfront or reoccurring multiple access fee , ( 3 ) a lifetime access fee , or ( 4 ) some combination thereof . another embodiment includes a method for building a database via the internet . for example , a user may utilize a portable tissue characterization device in accordance with one or more of the embodiments described herein to obtain individual testing data for one or more specific anatomic regions . the user may then transmit via the internet that data to a remote device ( such as a computer ) via a hardwired or wireless connection . alternatively , the device may dock in a cradle that can facilitate the transmission of the data to the remote device . the remote device may perform a data analysis ( 1 ) based on one or more analytical methods as described herein , ( 2 ) based on comparative data ( from the earlier testing by the same patient ) stored on the remote device , and / or ( 3 ) based on comparative data from a look - up table downloaded from a central server or generated in house from multiple patients . a response value can then be generated and communicated to the user or other party based on the data analysis . alternatively , the user transmits testing data to a central server for data analysis . the central server receives the data and perform the data analysis ( 1 ) based on one or more analytical methods as described herein , ( 2 ) based on comparative data ( from the earlier testing by the same patient ) previously or contemporaneously transmitted to the central server , and / or ( 3 ) based on comparative data from a look - up table stored at the central server . a response value can then be generated and communicated to the user or other party based on the data analysis . this embodiment also describes a method for generating a database of patient testing data . as more subscribers / users send in patient information the database grows and provides more accurate response values . a patient &# 39 ; s information may be stored indefinitely , such that the patient may at a subsequent date submit new testing data for comparison with his / her previous data and / or with the entire database . access is preferably based on a subscription fee . in order to receive a response value , new data ( including biographical information about the corresponding patient ) would need to be submitted in along with the access fee . it should be understood that prolific users , such as a doctor &# 39 ; s office , may pay a one time fee or set fee for multiple analyses . another embodiment allows the patient to access his / her data for monitoring , further analysis or review at anytime via the internet . the patient may use the portable device in the privacy of his / her own home without the need for nurses , doctors or technicians . thus , in an alternative embodiment , the patent &# 39 ; s data may be transmitted either as described above directly to the central server from the patient &# 39 ; s remote device to receive a response value or to a doctor who can monitor and forward the data to the central server for the response value . alternatively , the patient may submit the data directly to the central server and the doctor may monitor by accessing the central server with the permission of the patient . this allows remote diagnosis , which may be more desirable for both the patient and doctor . while preferred specific embodiments of this invention have been shown and described , modification thereof can be made by one skilled in the art without departing from the spirit or teaching of this invention . the embodiments described herein are exemplary only and are not limiting . many variations and modifications of the compositions and methods are possible and are within the scope of this invention . for example , while a piezoelectric element has been described as the actuator , a spring - loaded device , such as sold by activator methods inc . for chiropractors , may alternatively be used as the actuator . depending on in vivo measurements , another alternative , an acoustic actuator , may be used . additionally , the data acquisition system may comprise a system built on a cart , which includes a computer , software , storage and networking capabilities . alternatively , the data acquisition system may comprise a portable system connectable to a laptop computer for field usage and possible on - site data analysis . thus , the scope of protection is not limited to the embodiments described herein , but is only limited by the claims , the scope of which shall include all equivalents of the subject matter of the claims .
methods and apparatus for non - invasively assessing physiological hard of soft tissue of human and other species are described . in a preferred embodiment , tissue is vibrationally stimulated in vivo through a frequency spectrum . the tissue reacts against the stimulus and the reaction is preferably measured and recorded . based on analytical algorithms or comparisons with previously taken measurements , changes within the tissue can be detected and used for diagnostic purposes . further embodiments describe the usage of the device and methods for in vivo intra - operative and post - operative implant evaluations and as a therapeutic tool .