About CathyVI

CathyVI · ‎03-24-2023

@pink_poodle Am assuming you didn't get the AIC as part of your output but if you did let me know. I did not get AIC and this is one of the output am looking for.

CathyVI · ‎03-22-2023

Hello, Am getting a warning from my 'Heterogenous Compound Symmetry' code and ods output will not generate the output (AIC, BIC etc.) I wanted. I have used the same code for other covariance structures and it worked perfectly. I used the same code I posted below for another dataset with CSH and it worked. Here is a sample of my dataset. I would appreciate any help. Thanks data have; input Obs ID Group pre post1 post2 post3 time wk @@; cards; 1 14789 2 10 10 15 15 1 10 2 14789 2 10 10 15 15 2 10 3 14789 2 10 10 15 15 3 15 4 14789 2 10 10 15 15 4 15 5 15170 2 16 20 20 20 1 16 6 15170 2 16 20 20 20 2 20 7 15170 2 16 20 20 20 3 20 8 15170 2 16 20 20 20 4 20 9 15611 1 10 20 18 22 1 10 10 15611 1 10 20 18 22 2 20 11 15611 1 10 20 18 22 3 18 12 15611 1 10 20 18 22 4 22 13 15697 1 10 10 10 11 1 10 14 15697 1 10 10 10 11 2 10 15 15697 1 10 10 10 11 3 10 16 15697 1 10 10 10 11 4 11 17 15844 1 21 26 26 27 1 21 18 15844 1 21 26 26 27 2 26 19 15844 1 21 26 26 27 3 26 20 15844 1 21 26 26 27 4 27 21 15846 2 26 29 28 24 1 26 22 15846 2 26 29 28 24 2 29 23 15846 2 26 29 28 24 3 28 24 15846 2 26 29 28 24 4 24 25 15998 2 10 10 10 12 1 10 26 15998 2 10 10 10 12 2 10 27 15998 2 10 10 10 12 3 10 28 15998 2 10 10 10 12 4 12 29 16021 1 11 13 12 15 1 11 30 16021 1 11 13 12 15 2 13 31 16021 1 11 13 12 15 3 12 32 16021 1 11 13 12 15 4 15 33 16961 2 14 20 18 18 1 14 34 16961 2 14 20 18 18 2 20 35 16961 2 14 20 18 18 3 18 36 16961 2 14 20 18 18 4 18 37 17129 2 14 16 18 15 1 14 38 17129 2 14 16 18 15 2 16 39 17129 2 14 16 18 15 3 18 40 17129 2 14 16 18 15 4 15 41 17571 2 17 21 23 . 1 17 42 17571 2 17 21 23 . 2 21 43 17571 2 17 21 23 . 3 23 44 17571 2 17 21 23 . 4 . 45 17578 1 21 40 37 43 1 21 46 17578 1 21 40 37 43 2 40 47 17578 1 21 40 37 43 3 37 48 17578 1 21 40 37 43 4 43 49 17939 1 16 13 14 14 1 16 ; run; SAS code: /* Heterogenous Compound Symmetry CSH */ proc mixed data=have; class id group time; model wk=id time time*group; repeated time / type = csh subject = id r rcorr; ods output FitStatistics=FitCSH(rename=(value=CSH)) Dimensions=ParmCSH(rename=(value=NumCSH)); run; WARNING: Output 'FitStatistics' was not created. Make sure that the output object name, label, or path is spelled correctly. Also, verify that the appropriate procedure options are used to produce the requested output object. For example, verify that the NOPRINT option is not used.

CathyVI · ‎03-13-2023

@Kurt_Bremser There was an error when i ran your code 61 data long; 62 set have; 63 array all {*} a--d; 64 do i 0 1 to dim(all); - 73 200 ERROR 73-322: Expecting an =. ERROR 200-322: The symbol is not recognized and will be ignored. 65 x = all{i}; 66 output; 67 end; 68 keep x; 69 run; NOTE: The SAS System stopped processing this step because of errors. WARNING: The data set WORK.LONG may be incomplete. When this step was stopped there were 0 observations and 1 variables. NOTE: DATA statement used (Total process time): real time 0.02 seconds cpu time 0.00 seconds

CathyVI · ‎03-13-2023

@PaigeMiller I want the total number of 1s and the total number of zeros. This is what I showed in my output table. The table shows total number of 1's and 0's. This is how I want my output to look like. Thanks

CathyVI · ‎03-13-2023

Hello, I would like to sum my data set and use a proc tabulate to get sum for 1's and 0's. For example: data Have; input A B C D; datalines; 1 0 1 1 0 1 1 1 1 1 0 0; run; I used the code below but its only gives me value when total =1 or 0, not the value of all totals. data want; set have; if a=1 | b=1 |c=1|d=1 then total=1; else total=0; run; I would like my output to have total of a,b,c,d when its 1, and when its 0. For example Output Total N % 1 8 66.6 0 4 33.4 All 12 100

CathyVI · ‎03-03-2023

@PaigeMiller Thank you for your observation. This is the sample data. I have more than 100 rows. How do I find the average post-intervention umaqs score for each patient? It is the "each patient" that's getting me confused. data wide1; input ID$ Group umaqs0 umaqs6 umaqs12 umaqs26; datalines; 14789 2 10 10 15 15 15170 2 16 20 20 20 15611 1 10 20 18 22 15697 1 10 10 10 11 15844 1 21 26 26 27 ;run;

CathyVI · ‎03-03-2023

Hi, I have a dataset called wide1 with the following id variables which represent each individual. Here is a sample of my data below. How would I find the average post-intervention umaqs score for each patient? umaqs0 is measure before intervention, umaqs6- umaqs26 is measure post intervention. I want using the code below but am not sure if this is right. proc means data=wide1; var umaqs6 umaqs12 umaqs26; run; Obs ID Group _NAME_ umaqs0 umaqs6 umaqs12 umaqs26123456 14789 2 umaqs 10 10 15 15 15170 2 umaqs 16 20 20 20 15611 1 umaqs 10 20 18 22 15697 1 umaqs 10 10 10 11 15844 1 umaqs 21 26 26 27 15846 2 umaqs 26 29 28 24

CathyVI · ‎03-01-2023

@Reeza The dates does not have to be in all but I think since I need first visit date + 90 days to identify the second visit date, i will need first date to be in dates. I was trying my hands on this code but I don't know how to make first visit date in to date. What am i missing here? data want; set have; if (visit_dt)>index_dt then first_visit=1; else first_visit=0; run;

CathyVI · ‎03-01-2023

@Reeza @PeterClemmensen Am sorry, it a typo. You are right the first visit should be 4/5/2020, which is the first visit after the index date. How do I get this first visit and 90 days after first visit?

CathyVI · ‎03-01-2023

@Reeza @PeterClemmensen Thank you for your responses. This is what I expected in my output. For example, Jones first visit is 05/08/2019 and if the time between the first and second visit is more than 90days from the first visit then second visit =1. This will be my retention in care outcome because it met our criteria, if not then second visit=0. Expected output Obs pt_id Name index_dt visit_dt followup_365 first_visit second_visit 1 3 Jones 03/08/2019 05/08/2019 03/07/2020 05/08/2019 0 2 3 Jones 03/08/2019 02/19/2020 03/07/2020 05/08/2019 1 3 3 Jones 03/08/2019 09/08/2020 03/07/2020 05/08/2019 1 4 1 Smith 04/30/2019 04/30/2019 04/29/2020 04/30/2019 0 5 1 Smith 04/30/2019 05/03/2019 04/29/2020 04/30/2019 0 6 1 Smith 04/30/2019 07/23/2020 04/29/2020 04/30/2019 1 7 1 Smith 04/30/2019 12/30/2020 04/29/2020 04/30/2019 1 8 2 Williams 03/05/2020 02/05/2020 03/05/2021 02/05/2020 0 9 2 Williams 03/05/2020 04/05/2020 03/05/2021 02/05/2020 0 10 2 Williams 03/05/2020 10/15/2020 03/05/2021 02/05/2020 1 Thank you again, looking forward to your responses.

CathyVI · ‎02-28-2023

Hello, I am working on a data that I combined data and would like to identify the first visit date and second visit (90 days apart) from the index date. I have created a dummy data for assistance. I have created a index date, visit date and follow_up date (index date+365). I want to find the first visit date after the index date for each patient then find the second visit which will be 90 days from the first visit. this is my attempt but am not getting it. I will appreciate any help with the coding. data have; input pt_id $ Name $ index_dt :mmddyy10. followup_365 :mmddyy10. visit_dt :mmddyy10.; format followup_365 index_dt visit_dt mmddyy10.; datalines; Obs pt_id Name index_dt visit_dt followup_365 1 3 Jones 03/08/2019 05/08/2019 03/07/2020 2 3 Jones 03/08/2019 02/19/2020 03/07/2020 3 3 Jones 03/08/2019 09/08/2020 03/07/2020 4 1 Smith 04/30/2019 04/30/2019 04/29/2020 5 1 Smith 04/30/2019 05/03/2020 04/29/2020 6 1 Smith 04/30/2019 07/23/2020 04/29/2020 7 1 Smith 04/30/2019 12/30/2020 04/29/2020 8 2 Williams 03/05/2020 02/05/2020 03/05/2021 9 2 Williams 03/05/2020 04/05/2020 03/05/2021 10 2 Williams 03/05/2020 10/15/2020 03/05/2021 ; run; This is my two attempts but not working proc sql; create table chc_master as select pt_id, Name, index_dt,followup_365, visit_dt, min(visit_dt) as first_visit format mmddyy10. from have; quit; data chc_master; set have; if first.pt_id then do; if first.index_dt ne . & first.visit_dt ne . then first_visit=1; else first.visit=0; end; run;

CathyVI · ‎02-20-2023

Hello, I have a dataset with over 10,000 drugs list. I will like to categorize each drug into specific group. However, some drugs are likely to be in 2 or 3 groups e.g., Prochlorperazine could be antipsychotics or antiemetics. I want Prochlorperazine to be group as antipsychotics if found and also antiemetics if found so that Prochlorperazine will be in both groups just like the table below. I have tried to use the FIND function or the UPCASE =: to group the drugs but this methods could not assign a drug into different group if that's the case for the drug. SAS only output the drug into the last category. Here is what I want to achieve with table of drugs and the categories. Drug Category Prochlorperazine Antipsychotics Prochlorperazine Antiemetics Celecoxib Analgesics Celecoxib Anti-inflammatory Agents Diclofenac Anti-inflammatory Agents Diclofenac Dermatological Agents Diclofenac Ophthalmic Agents Diclofenac Analgesics This is my code but did not solve this problem. Or is this logic impossible to do? data check2; length category $80; set check1; if upcase (drug) =: "Celecoxib" then do; category = "Analgesics"; end; else if upcase (drug) =: "Celecoxib" then do; category = "Anti-inflammatory Agents"; end; run;

CathyVI · ‎02-19-2023

Hello, I created 3 variables x1,x2,x3 with variance =9 from a simulate of 1000 copies. Now I want to: 1. Confirm that the empirical variances of (x1+x2) is consistent with Var x1 +x2 =Var (x1) +Var (x2) + 2Cov (x1, x2) 2. Confirm that the empirical covariance of (x1+x2) and (x1 + x3) is consistent with Cov (x1 +x2, x1 +x3 ) = Cov (x1, x1) +Cov (x1, x3) +Cov (x2, x1) +Cov (x2,x3). It seems I can use PROC GLIMMIX but am not sure how to use the procedure. How do I achieve this in sas. Thanks

CathyVI · ‎02-19-2023

Hello, Am trying to simulate 1000 copies of three variables x1, x2, and x3 with variances equal to 9. I can choose any mean. How do I do this in sas? I created this simulation code but I don't know how to include variance =9 or chose any mean. Thanks for your help. data a; do i=1 to 1000; y1=rannor(0); y2=rannor(0); y3=rannor(0); output; end; run;

CathyVI · ‎02-09-2023

Thank you both. I wish I could accept both solutions because they both provided me a better understand of how to address the problem. @ballardw Thank you, you provide the solution and refined my initial thought process. @SASJedi Thank you so much. I accepted your solution because when I looked into the entire data, I found some drug name with weird name. For example 0.9%ESOMEPRAZOLE or HCL ESOMEPRAZOLE. I was able to use the find function to identify these ones.

Online Status	Offline
Date Last Visited	‎10-16-2025 09:08 AM

Re: check if the content of a variable is in another variable

Re: check if the content of a variable is in another variable

check if the content of a variable is in another variable

Find the count using proc tabulate

Re: Creating subtotal in Proc report

Creating subtotal in Proc report

Re: Bold subheading in proc report

Re: Bold subheading in proc report

Bold subheading in proc report

Re: merging data sets using set statement

Re: check if the content of a variable is in another variable

Re: check if the content of a variable is in another variable

Re: check if the content of a variable is in another variable

Re: check if the content of a variable is in another variable

Re: Bold subheading in proc report

Using proc report

How do I analyze the effect of repeated dose

Proc freq with multiple variables

variable categorization

Identifing date of encounters

Re: Heterogenous Compound Symmetry

Heterogenous Compound Symmetry

Re: sum multiple variables

Re: sum multiple variables

sum multiple variables

Re: Find the average post-intervention umaqs score for each patient

Find the average post-intervention umaqs score for each patient

Re: Identify first and second visit after index date

Re: Identify first and second visit after index date

Re: Identify first and second visit after index date

Identify first and second visit after index date

categorize a response into multiple group

Empirical variances

Simulate 1000 copies of three variables with variances

Re: Identify variables in the dataset

SAS Analytics Explorers