дифференцируя значения одной и той же переменной я

Question

Я выполняю логистическую регрессию в R, и я пытаюсь построить логит вероятности против вероятности получения 1. Я хотел бы построить все значения, предсказанные как 1 (положительные), используя один цвет, а значения, предсказанные как0 (негативы) другим цветом.Это означает, что значения> 0,5 отображаются одним цветом, а значения <0,5 - другим цветом.Есть идеи, как мне это сделать?Вот мой код: </p>

     pdgng<-data$pdgng
     ec<-data$ec
     logitp <- 0.497-1.699 * (log(pdgng)) +3.829 * (log(ec))
     logistic<-exp(logitp)/(1+exp(logitp))
     op5<-par(cex.lab=0.9,font.lab=2,cex.axis=1,bty="n")
     plot(logitp,logistic,ylab="Probability",xlab="logitp"       
     abline(h=0.5, col="blue",lwd=1,lty=2)

Я пытался использовать цикл, но не могу понять, как применить его к моему делу.С уважением, Антонела

Answer 1

Спасибо за ваш ответ ... Вот полный код и данные.

pregnancystate<-c("pregnancystate_nomales.csv")
data<-read.csv(file = "pregnancystate_nomales.csv", dec = ".",sep = ",", header = 

TRUE)

 cat(file = "preg.bug", "
  #Likelihood:
   model {
      for(i in 1: 69){
       pregnancy[i] ~ dbern(p[i]) 
       logit(p[i]) <- b0+ b1 * (log(pdgng[i])) + b2 * (log(ec[i])) 
           }
         #priors:
       b0 ~ dunif(0,1)
       b1 ~ dnorm(0, 0.0001)
       b2 ~ dnorm(0, 0.0001)

          }"
            )

          pdg<-data$pdgng
          ec<-data$ec
    obs<-69   #number of observations

   inits <- function() {
    list(b0=runif(0,0.05),b1 = rnorm(1,0,1), b2 = rnorm(1,0,1))
       }

      parameters <- c( "b0", "b1", "b2")
       ni <- 100000
        nt <- 1
         nb <- 50000
          nc <- 3

         library(jagsUI)
         pregn <- jags(data=data,  parameters.to.save = parameters, model.file 
          ="preg.bug" , n.chains = nc, n.thin = nt,
           n.iter = ni, n.burnin = nb)
           print(pregn)

     ###Logistic curve
     pdgng<-data$pdgng
       ec<-data$ec
       logitp <- 0.497-1.699 * (log(pdgng)) +3.829 * (log(ec))
       logistic<-exp(logitp)/(1+exp(logitp))

        op5<-par(cex.lab=0.9,font.lab=2,cex.axis=1,bty="n")
         plot(logitp,logistic,ylab="Probability",xlab="logitp")
          abline(h=0.5, col="blue",lwd=1,lty=2)

Мои данные: беременности, pdgng, ec, logit(р), вероятность 0,143997.3937,746.5102301,2.736380058,0.939139522 1,45109.3079,1418.995961,4.659342216,0.990616199 0,52683.58472,56.85769855, -0.802104453,0.309575536 1,138659.7743,852.4850646,2.984848323,0.951884912 0,52689.26541,47.78030436, -1,091276892, 0,251377908 0,35554.31237,54.64659066, -0.578210329,0.359344501 1,44957.50427,881.71069,3.870957685,0.979586972 0,17346.58536,34.3035689, -0.82320263,0.305084254 1,87056.00603,959.375944,3.524232772,0.971369456 0,65611.00906,75.41128029, -0.494460715,0.3788433020,40403.69619,50.41093373, -0.806534415,0.308629482 1,70574.46512,201.8534225,1.088218252,0.748046058 0,42819.85124,63.11200154, -0.47587338,0.383227039 0,24432.8854,55.57544032, -0.273726259,0.431992534 0,40603.1634,66.64685283, -0.346121395,0.414323288 0, 42019,87914,52,39963524, -0,771133528,0,316233952 0,31035,09,74,1115375,0,028378308,0.507094101 1,72415.54675,245.6925039,1.395895237,0.801531717 1,84035.46785,125.7413611,0.172889837,0.543115117 0,54189.70386,72.94486225, -0.408776022,0.399205644 1,27718.06513,59.68819086, -0.248043693,0.438305071 0,34963.48171,61.44647007, -0,370936316,0,408314794 0,107577.8631,100.9496181, -0.374125172,0.407544612 0,45300.54732,74.97611356, -0.23107752,0.442486313 1,87096.67182,967.8444156,3.538495187,0.97176345 0,51185.37663,139.2228995,0.707560153,0.669861819 0,55756.69828,54.10187088, -0.926462682,0.28364291 0, 62001.79489,72.06355191, -0.528232186,0.370929297 1,90068.53164,1229.299153,3.911209047,0.980376504 1,49585.43419,167.4432442,1.037725477,0.738410897 1,49404.23826,89.63236934,0.001765223,0.500441306 1,61502.94127,665.0569165,3.171324365,0.959740787 1,67251.66979,233.6405586, 1.366821986,0.796866212 1,92243.7537,327.5783641,1.695613274,0.844960934 0,54199.09589,49.26495746, -1.06124114,0.257072343 1,123323.012,1444.722047,3.94799264,0.98107228, 12247644582,0.206255721 0,29770.49904,54.5633406, -0.449895094,0.389385709 1,109766.4443,780.8332001,3.011143597,0.953075026 1,98604.30654,172.0454198,0.576164889,0.640184471 1,312081.4201,215.5575672,0.101773998,0.52542156 0,57012.73092,55.32276656, -0.905791033,0.2878618941,60997,32874,384.736727,2.267748818,0.906170555 1,97002.76256,214.6832766,0.956219147,0.722364177 1,153642.8724,1119.772165,3.362507665,0966512036 0,38540.3286,225 525 532 560660660560560560560560560560560560560560560560560560560560560560560560560560560560560% Папка городского телевидения73.12316073, -0.182908088,0.454400038 1,231120.767,555.4809947,1.896414313,0.869485158 1,69794.21866,192.9932648,1.021808213,0.735324218 1,38155,6815,155.986,239,225, 0,837255778 1,81997,44506,602.4558705,2,795054752,0,942408006 0,51168.54791,54,91864886, -0,838270405,0,301899182 1,45249,61128,277.0660431,1.942187893,0.87459231, 91,8198666, 195, 195, 19, 19, 19, 19, 19, 197548945,0.275780463 1,102338.8059,986.4792596,3.451336152,0.969270963 0,36239.49787,24.65023113, -1.915433194,0.128371691, -1.93347035,18.41201813, -1.93, 231, 231, 231, 231, 231, 231, 231, 231, 231, 231, 231, 231, 1932, 17, 30, 231, 231, 231, 231, 1932, 1932, 1900, 153, 171, 173, 171, 171, 171, 1900, 153, 171, 1900, 153, 171, 171, 171, 171, 1900, 1700, 153, 153, 161, 163, 163, 163, 163, 163, 163, ППК, 0,233192011 0,30118,85279,31.41423431, -1,376085775,0.201638377 0,24570,3708,34.2774031, -1.081053394,0.253306723 0,24154.21332,36,994861, -0,94166181,0,280564787 1,925034690068868560570560560560560560560560560560560560560560560560560560560560605605605485605485485485485485485485485485485485485485485485485485295295295295295295291302, 0,233192011,30118,85279,31,414234311,16917,84884,260,5582078,2,565165807,0,928585784 1,36961,88734,338.2647117,2.422986594,0.918563435

1 - беременные женщины, а 0 - не беременные женщины.Я намереваюсь войти против вероятности и использовать другой цвет для точек, которые принадлежат беременным женщинам, и другой цвет для тех, которые принадлежат небеременным.Надеюсь, это прояснит мой вопрос.С уважением, Антонела