SAS/STAT 15.3 User's Guide The NLMIXED Procedure Example 89.1 One-Compartment Model with Pharmacokinetic Data https://go.documentation.sas.com/doc/en/pgmsascdc/9.4_3.5/statug/statug_nlmixed_examples01.htm Details --> Compartment Models https://go.documentation.sas.com/doc/en/pgmsascdc/9.4_3.5/statug/statug_nlmixed_details07.htm Syntax --> CMPTMODEL Statement https://go.documentation.sas.com/doc/en/pgmsascdc/9.4_3.5/statug/statug_nlmixed_syntax05.htm The MCMC Procedure (Bayesian nonmem) Example 79.22 One-Compartment Model with Pharmacokinetic Data https://go.documentation.sas.com/doc/en/pgmsascdc/9.4_3.5/statug/statug_mcmc_examples31.htm Moreover SAS can do : Non-Compartmental (Model Independent) methods (but no examples found in doc)!   Also, ... found an old paper (2007) ... probably still useful. Paper 148-2007 Population Pharmacokinetics Studies with Nonlinear Mixed Effects Modeling David Shen, Alps Biostatistics Inc. ; Zaizai Lu, AstraZeneca Pharmaceuticals https://support.sas.com/resources/papers/proceedings/proceedings/forum2007/148-2007.pdf Good luck, Koen ... View more

Is your question in the realm of population pharmacokinetics/pharmacodynamics (pop PK/PD) analysis used in Clinical Pharmacology research? Or are you talking about Nonlinear Mixed Effects Modeling (NONMEM) in general? Are you comparing the NONMEM software (ICON Development Solutions) to SAS/Stat PROC NLMIXED? Do you need Ordinary Differential Equations (ODE)? Are you fitting compartment Models? If yes, how many compartments? Which specific dosing (aka forcing functions) do you need?   If you give a bit more insight in your use case, then the answers provided will be more insightful for you as well.   Koen ... View more

@bhr-q wrote: I agree with you, I actually tried to do that with right click on your answer but couldn't find that. if there is anyway let me know to do that. Thank you so much all of you for your help and sorry for the mistake. You should not look too hard (to change solutions). For me, it's okay. It's fine as it is. Most important thing is that your original question is solved! Moreover, the addition of @Quentin (on the counter-intuitive aspect) is also relevant. Many people will recognize themselves in it. 😉   Koen ... View more

You can also use this formula to get predictions on new data (use a data step where t_value is your input-time-value to be scored): predicted_target = 1/(1+exp(-((t_value - intercept) / scale_parameter))); Please verify by scoring the training data and comparing with proc lifereg predictions. I won't put my hand up to it, but I think that's the right formula.   Koen ... View more

What you CANNOT do: use PROC SCORE use the STORE statement (PROC LIFEREG) to store the fitted model in an item store. Then use PROC PLM to create predictions on new data. I think the PLM scoring possibility is disabled for models coming from LIFEREG. If not, don't trust the predictions (they are wrong). What you CAN do: specify a missing value for the response variable for each observation that you want to score but not use for fitting the model.  Since such observations are not used to fit the model, it does not matter whether those observations contain a non-zero weight. specify a weight of 0 for observations that you want score but not use for fitting the model. (PROC LIFEREG has a WEIGHT statement) I agree it's not optimal since every time you want to score another set of data ... you need to fit your model again. See also this blog: The missing value trick for scoring a regression model By Rick Wicklin on The DO Loop February 17, 2014 https://blogs.sas.com/content/iml/2014/02/17/the-missing-value-trick-for-scoring-a-regression-model.html   The below blog does not contain a solution to your question, but I want you to know about it anyway: Interpret estimates for a Weibull regression model in SAS By Rick Wicklin on The DO Loop October 27, 2021 https://blogs.sas.com/content/iml/2021/10/27/weibull-regression-model-sas.html   Koen ... View more

https://communities.sas.com/t5/Statistical-Procedures/Difference-in-difference-for-gamma-distribution/td-p/623102   For difference in difference analysis with non-normal responses, see this note.    Koen ... View more

When I submit proc setinit; run; I get ---SAS/GRAPH 14JUL2024 (CPU A)   When I submit proc product_status; run; I get For SAS/GRAPH ... Custom version information: 9.4_M8 M8 = Maintenance Release 8 (another maintenance release will do as well of course). M6 or above is preferred though.   I think you need to see something similar in order to be able to turn ods graphics ... off .   Koen ... View more

This is the same question as posted in the programming board.   Home > Programming > Programming > Difference-in-difference https://communities.sas.com/t5/SAS-Programming/Difference-in-difference/m-p/927128 Please do not duplicate your question on multiple boards. Also take the right board. I will move this one to "Statistical Procedures"-board in a minute.   Koen ... View more

stacked bar chart in SAS site:blogs.sas.com stacked bar chart in SAS site:support.sas.com Google Search.   Koen ... View more

Hello, https://support.sas.com/sassamples/graphgallery/PROC_SGPLOT.html Scroll down to: Sample 55866 - Use the SEGLABEL option to label bar segments in a bar chart https://support.sas.com/sassamples/graphgallery/PROC_GCHART_Graph_Types_Charts_Bar.html Scroll down to:  Sample 48451 - Center the subgroup labels in a horizontal bar chart created with PROC GCHART https://support.sas.com/kb/24/880.html Sample 24880: Reorder the subgroups in a stacked bar chart produced with PROC GCHART Koen ... View more

Hello, What SAS version are you using? Is it SAS 9.4 or SAS Viya?   Submit  %put &=SYSVLONG4; to find out (see LOG-screen for result).   To see the expiration dates for your SAS modules, you can submit: proc setinit; run;  /*SAS 9.4*/  ... or cas; proc cas; getLicensedProductInfo; quit; /*SAS Viya*/ When you cannot turn ODS GRAPHICS ... OFF , you can still use the ODSTITLE= option behind the forward slash on a pchart (p-chart) statement to get what you want.   Koen ... View more

I cannot see the LOG of your code after submission, but 2 things might have gone wrong here: 1. are your quarters (the character version) strings like "2022Q2" (yyyy for the year then a letter Q and then the quarter number)? 2. you need vref=0.9 instead of href=0.9.  I do not see how 0.9 can be a quarter-indication, but 0.9 can be a proportion. On the VREF= option, you specify vertical axis values, on the HREF= option you specify horizontal axis values. Like @Quentin , you might be confused about this 😉.   Koen ... View more