How to force the last observation of a variable

Accepted Solution Solved
Reply
Super Contributor
Posts: 303
Accepted Solution

How to force the last observation of a variable

Hello, in the following data generation; I want to change the last observation for variable 'censored' to '0' when ever it is '1'. That is the 100th observation should have censored=0 , even if it is supposed to e 1

data test;

seed=-1;

  alpha1 = 2.50;

  beta1 = 9.50;

  do i = 1 to 100;

   lambdaT = 0.025; *baseline hazard;

   lambdaC= .03*light=0.15; *heavy=0.03;

   er=0+sqrt(0.0001)*rannor(1);

    t = rand("WEIBULL", 0.75, lambdaT);  * time of event;

    c = rand("WEIBULL", 1.25, lambdaC) ;* time of censoring;

    time = min(t, c);    * which came first?;

    censored = (c lt t);

      obs=(t lt c); * creating observation variable from censored when

      y= alpha1 + beta1*t + er;

    output;

  end;

run;


Accepted Solutions
Solution
‎07-15-2013 10:39 AM
Super Contributor
Posts: 1,636

Re: How to force the last observation of a variable

Posted in reply to desireatem

data test;

seed=-1;

  alpha1 = 2.50;

  beta1 = 9.50;

  do i = 1 to 100;

   lambdaT = 0.025; *baseline hazard;

   lambdaC= .03;  *light=0.15; *heavy=0.03;

   er=0+sqrt(0.0001)*rannor(1);

    t = rand("WEIBULL", 0.75, lambdaT);  * time of event;

    c = rand("WEIBULL", 1.25, lambdaC) ;* time of censoring;

    time = min(t, c);    * which came first?;

    censored =ifn(i ne 100,(c lt t),0);

      obs=(t lt c); * creating observation variable from censored when

      y= alpha1 + beta1*t + er;

    output;

  end;

run;

View solution in original post


All Replies
Solution
‎07-15-2013 10:39 AM
Super Contributor
Posts: 1,636

Re: How to force the last observation of a variable

Posted in reply to desireatem

data test;

seed=-1;

  alpha1 = 2.50;

  beta1 = 9.50;

  do i = 1 to 100;

   lambdaT = 0.025; *baseline hazard;

   lambdaC= .03;  *light=0.15; *heavy=0.03;

   er=0+sqrt(0.0001)*rannor(1);

    t = rand("WEIBULL", 0.75, lambdaT);  * time of event;

    c = rand("WEIBULL", 1.25, lambdaC) ;* time of censoring;

    time = min(t, c);    * which came first?;

    censored =ifn(i ne 100,(c lt t),0);

      obs=(t lt c); * creating observation variable from censored when

      y= alpha1 + beta1*t + er;

    output;

  end;

run;

🔒 This topic is solved and locked.

Need further help from the community? Please ask a new question.

Discussion stats
  • 1 reply
  • 150 views
  • 0 likes
  • 2 in conversation