Skip to content
Snippets Groups Projects
Commit 9ecbaca3 authored by Mildred Alexandra Arias Yanez's avatar Mildred Alexandra Arias Yanez
Browse files

Tarea_de_Física_Médica

parent 1dfd84bc
No related branches found
No related tags found
No related merge requests found
%% Cell type:markdown id: tags:
#### Con los datos de la tabla 3 del artículo, verifique que la frecuencia de dicéntricos más anillos por célula aumenta con el cuadrado de la dosis de radiación inmediatamente después del accidente y 15 meses después del accidente.
%% Cell type:code id: tags:
``` python
import numpy
import matplotlib.pyplot as plt
import scipy.stats as stats
```
%% Cell type:code id: tags:
``` python
x = [7,6.2,4.4,4.3,3,2.9,2.7,2.1,1.3,1.1,1,1,0.3,0.2]
y= [1.253,1.030,0.824,0.540,0.270,0.262,0.228,0.160,0.070,0.060,0.053,0.050,0.02,0.005]
z= [0.26,0.18,0.06,0.14,0.04,0.08,0.04,0.04,0.04,0,0.02,0.02,0.02,0]
poly = numpy.polyfit(x, y, 2)
polynomial = numpy.poly1d(poly)
plt.scatter(x, y)
xs = numpy.linspace(x[0], x[-1], len(x))
y_int = polynomial(xs)
plt.xlabel("dosis")
plt.ylabel("Inmediatamente despues del accidente")
plt.plot(xs, y_int)
plt.show()
print(polynomial)
poly = numpy.polyfit(x, z, 2)
polynomial = numpy.poly1d(poly)
plt.scatter(x, z)
xs = numpy.linspace(x[0], x[-1], len(x))
z_int = polynomial(xs)
plt.xlabel("dosis")
plt.ylabel("15 meses después del accidente")
plt.plot(xs, z_int)
plt.show()
print(polynomial)
```
%% Output
2
0.01723 x + 0.06768 x - 0.03289
2
0.004728 x + 0.0002015 x + 0.01276
%% Cell type:markdown id: tags:
Entonces el modelo indica que la frecuencia del dicéntrico aumenta efectivamente con el cuadrado de la dosis inmediatamente.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment