je ne suis qu'un jeune doctorant qui débute avec R et vous remercie par ailleurs pour votre indulgence sur la maîtrise du langage R.
Mon problème:
J'ai 96 observations d'AUC suite à l'administration de doses croissantes d'un médicament sur 6 chiens. Chaque dose étant donnée chaque semaine, la concentration plasmatique étant < LOQ en moins de 24h, on considère qu'il n'y a pas d'influence de la dose précédente sur l'AUC suivant.
Le processus d'absorption étant via transporteurs actifs, on s'attend à une saturation des transporteurs au-delà d'un seuil à définir. Il n'a pas été trop compliqué de démontrer la non linéarité AUC/Dose. Maintenant, je cherche à savoir si à partir de ces données, je pourrais déterminer le seuil en question. En fouillant différents forums, je ne suis pas parvenu à le programmer.
J'ai essayé ce modèle, qui fonctionne très bien mais il est faux, ne prenant pas en considération la variable aléatoire "chien". Et si j'applique cette formule à chaque chien pris séparément,le modèle ne fonctionne pas.
Code : Tout sélectionner
LP<- function (b, Theta0, Theta1, Theta2)
{
Y<-Theta0+Theta1*(b-Theta2)
Y[b>=Theta2]<-Theta0
return(Y)
}
LP
Fit<-nls(log(AUC)~LP(POSO, Theta0, Theta1, Theta2), start=list(Theta0=5, Theta1=3, Theta2=0.5), data=da5)
Merci pour votre aide svp