Skip to content
Snippets Groups Projects
Commit 311cd154 authored by José Antonio López Rodríguez's avatar José Antonio López Rodríguez
Browse files

update

parent 2062bff6
No related branches found
No related tags found
No related merge requests found
File added
File added
%% Cell type:markdown id: tags: %% Cell type:markdown id: tags:
<img src="img/isotipo-RGB-pequeno.png" alt="" width="15%" />
# Mezcla de poblaciones celulares y resistencia a la radiación # Mezcla de poblaciones celulares y resistencia a la radiación
%% Cell type:code id: tags: %% Cell type:code id: tags:
``` python3.6 ``` python3.6
from pandas import read_csv from pandas import read_csv
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
import numpy as np import numpy as np
import random import random
``` ```
%% Cell type:markdown id: tags: %% Cell type:markdown id: tags:
## Ejercicio ## Ejercicio
Se dispone de un archivo de texto (poblacion.csv) con valores en un arreglo de dimensión 101,101 Se dispone de un archivo de texto (poblacion.csv) con valores en un arreglo de dimensión 101,101.
Los valores posibles del arreglo son los enteros del 0 al 3. Este arreglo podría estar relacionado con una simulación usando un modelo de Potts celular. Los valores posibles del arreglo son los enteros del 0 al 3. Cada valor corresponderá con un estado de una célula presente en el sitio.
En nuestro caso de interés los valores tienen la siguiente interpretación:
Este arreglo podría estar relacionado con una simulación usando un modelo de Potts celular.
* 0: no hay células
* 1: células aireadas En el caso de interés los valores tienen la siguiente interpretación:
* 2: células hipóxicas
* 3: células muertas * 0: no hay célula
* 1: célula aireada
* 2: célula hipóxica
* 3: célula muerta
El objetivo de este ejercicio es simular el efecto de la radiación sobre este conjunto de células, usando un modelo de blanco simple. El objetivo de este ejercicio es simular el efecto de la radiación sobre este conjunto de células, usando un modelo de blanco simple.
### Modelo: ### Modelo:
La probabilidad de daño a una determinada célula es $p=1-e^{-D/Do}$, donde $Do=1$ Gy para las células tipo 1 y $Do=5$ Gy para las células tipo 2. La probabilidad de daño a una determinada célula es $p=1-e^{-D/Do}$, donde $Do=1$ Gy para las células tipo 1 y $Do=5$ Gy para las células tipo 2.
### Instrucciones ### Instrucciones
* Represente gráficamente los datos del archivo, de forma que se pueda visualizar en que zona se encuentra cada tipo de célula. * Represente gráficamente los datos del archivo, de forma que se pueda visualizar en que zona se encuentra cada tipo de célula.
* Simule la aplicación de diferentes dosis a réplicas de su población de prueba en un rango de 1 a 25 Gy. Utilice el modelo proporcionado para decidir si cada célula sobrevive o muere. * Simule la aplicación de diferentes dosis a réplicas de su población de prueba en un rango de 1 a 25 Gy. Utilice el modelo proporcionado para decidir si cada célula sobrevive o muere.
* Calcule la fracción de supervivencia en cada caso, tabule y grafique $\ln(S)$ en función de $D$ * Calcule la fracción de supervivencia en cada caso, tabule y grafique $\ln(S)$ en función de $D$
* Tabule y grafique la proporción de células sobrevivientes tipo 2 con respecto al total de células sobrevivientes tipo 1 y 2 en función de la dosis. * Tabule y grafique la proporción de células sobrevivientes tipo 2 con respecto al total de células sobrevivientes tipo 1 y 2 en función de la dosis.
* Indique qué propiedad de la curva resultante evidencia la resistencia diferencial de las sub poblaciones 1 y 2. * Indique qué propiedad de la curva resultante evidencia la resistencia diferencial de las sub poblaciones 1 y 2.
* Usted puede usar esta hoja para desarrollar la respuesta o cualquier otra herramienta a su alcance, siempre que el resultado sea razonablemente reproducible. * Usted puede usar esta hoja para desarrollar la respuesta o cualquier otra herramienta a su alcance, siempre que el resultado sea razonablemente reproducible.
%% Cell type:code id: tags: %% Cell type:code id: tags:
``` python3.6 ``` python3.6
df = read_csv('poblacion.csv',header=None) df = read_csv('poblacion.csv',header=None)
``` ```
%% Cell type:code id: tags: %% Cell type:code id: tags:
``` python3.6 ``` python3.6
df df
``` ```
%% Output %% Output
0 1 2 3 4 5 6 7 8 9 ... 91 92 93 \ 0 1 2 3 4 5 6 7 8 9 ... 91 92 93 \
0 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ... 0 0 0
1 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 1 0 0 0 0 0 0 0 0 0 0 ... 0 0 0
2 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 2 0 0 0 0 0 0 0 0 0 0 ... 0 0 0
3 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 3 0 0 0 0 0 0 0 0 0 0 ... 0 0 0
4 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 4 0 0 0 0 0 0 0 0 0 0 ... 0 0 0
.. ... ... ... ... ... ... ... ... ... ... ... ... ... ... .. ... ... ... ... ... ... ... ... ... ... ... ... ... ...
96 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 96 0 0 0 0 0 0 0 0 0 0 ... 0 0 0
97 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 97 0 0 0 0 0 0 0 0 0 0 ... 0 0 0
98 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 98 0 0 0 0 0 0 0 0 0 0 ... 0 0 0
99 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 99 0 0 0 0 0 0 0 0 0 0 ... 0 0 0
100 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 100 0 0 0 0 0 0 0 0 0 0 ... 0 0 0
94 95 96 97 98 99 100 94 95 96 97 98 99 100
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0
2 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0
4 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0
.. ... ... ... ... ... ... ... .. ... ... ... ... ... ... ...
96 0 0 0 0 0 0 0 96 0 0 0 0 0 0 0
97 0 0 0 0 0 0 0 97 0 0 0 0 0 0 0
98 0 0 0 0 0 0 0 98 0 0 0 0 0 0 0
99 0 0 0 0 0 0 0 99 0 0 0 0 0 0 0
100 0 0 0 0 0 0 0 100 0 0 0 0 0 0 0
[101 rows x 101 columns] [101 rows x 101 columns]
%% Cell type:code id: tags: %% Cell type:code id: tags:
``` python3.6 ``` python3.6
plt.imshow(df, cmap='viridis') plt.imshow(df, cmap='viridis')
plt.colorbar() plt.colorbar()
plt.show() plt.show()
``` ```
%% Output %% Output
%% Cell type:code id: tags: %% Cell type:code id: tags:
``` python3.6 ``` python3.6
``` ```
......
2/img/isotipo-RGB-pequeno.png

18.1 KiB

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