About rogersaj

rogersaj · ‎12-14-2022

Hi, I have a dataset of women in labor. We have a repeated measure, cervical dilation for each woman, as well as the time between cervical dilation measurements. My goal is repeat what a similar paper did to find how the rate of cervical dilation differs with different patient characteristics. One prior study used Proc Mixed to model increasing cervical dilation parameter over time using an eighth-degree polynomial model. I'd like to: 1) Fit the data to the correct polynomial model 2) Impute the value of 6cm (clinically significant) for all women for whom that value was not recorded 3) Find out the average time taken to traverse from the first recorded cervical dilation to 6cm (eg. took 120 minutes to go from 3cm to 6cm). Here is a sample of my data; Patient ID Cervical dilation Difference in time between two consecutive cervical dilation checks (Minutes) 1 1 1 2 60 1 3 30 1 5 60 1 6 30 1 7 45 1 10 30 Thanks in advance!

rogersaj · ‎05-11-2018

Well, I think I found the answer to my own question: LMS values (measures of skew, the median, and the standard deviation) were computed from the interpolated cubic splines at weekly intervals. Cole’s procedures and an iterative least squares method were used to derive the LMS parameters (L = Box-Cox power, M = median, S = coefficient of variation) from the multicentre meta-analyses for weight, head circumference and length. The LMS splines were smoothed slightly while maintaining data integrity as noted above. Cole TJ, Green PJ: Smoothing reference centile curves: the LMS method and penalized likelihood. Stat Med. 1992, 11: 1305-1319. 10.1002/sim.4780111005.

rogersaj · ‎05-11-2018

This is slightly more in the realm of math questions than SAS, but maybe someone can help me figure out how to use SAS to compute it: The CDC publishes US infant growth charts for weight and length -- complete with percentile, z-scores, and "LMS parameters: the median (M), the generalized coefficient of variation (S), and the power in the Box-Cox transformation (L)" to allow us to calculate our own percentiles from weight and length. The chart starts at birth. However, I want these percentiles from BEFORE birth. The WHO just published these values for expected fetal weight at different gestational ages. These percentiles (5th, 10th, 25th, etc.) aren't good enough for me. I need a much more granular percentile -- i.e. I need the L, M, and S values so that I can calculate my own percentile from birthweight data. Any idea how to get those LMS values from the above table? The authors don't publish them. Do I need to contact them for raw data? That might take forever... Thanks SO much!

rogersaj · ‎02-07-2018

Thank you so very much. SAS had been installed on my new machine and by default was using SAS 9.4 (English with DCBS). Now that I've changed it to SAS 9.4 (Unicode support), my issue has been solved.

rogersaj · ‎02-07-2018

I really just want to set my dataset in my WORK library and analyze it without getting errors and warnings. My main concern is this line and the fact that it will say an inordinately small number of observations are in the dataset (When I know there are many, many observations). Data file CSL.CSL1.DATA is in a format that is native to another host, or the file encoding does not match the session encoding. Thanks for your comments on the UTF8 encoding. I took this part (which was a "solution" that I found on another person's post regarding the same error) out of my code. proc datasets library=csl nodetails nolist; modify csl1/ correctencoding='utf8'; quit; I'm still getting the same errors. libname csl 'C:\Users\Anna Joy Rogers\Documents\CSL'; proc contents data =csl.csl1; run; proc freq data =csl.csl1; tables ab_any; run; 4 proc contents data =csl.csl1; run; NOTE: Writing HTML Body file: sashtml.htm NOTE: Data file CSL.CSL1.DATA is in a format that is native to another host, or the file encoding does not match the session encoding. Cross Environment Data Access will be used, which might require additional CPU resources and might reduce performance. WARNING: Some character data was lost during transcoding in the dataset CSL.CSL1. Either the data contains characters that are not representable in the new encoding or truncation occurred during transcoding. NOTE: PROCEDURE CONTENTS used (Total process time): real time 1.36 seconds cpu time 0.54 seconds 5 proc freq data =csl.csl1; NOTE: Data file CSL.CSL1.DATA is in a format that is native to another host, or the file encoding does not match the session encoding. Cross Environment Data Access will be used, which might require additional CPU resources and might reduce performance. 6 tables ab_any; WARNING: Some character data was lost during transcoding in the dataset CSL.CSL1. Either the data contains characters that are not representable in the new encoding or truncation occurred during transcoding. 7 run; NOTE: There were 1 observations read from the data set CSL.CSL1. NOTE: PROCEDURE FREQ used (Total process time): real time 2.26 seconds cpu time 0.90 seconds Thanks so much for your time!

rogersaj · ‎02-02-2018

I know others have posted this same question, but even after reading the other posts and talking to SAS support, I can't solve the problem. Here's my code: libname csl 'C:\Users\rogersaj\Documents\CSL'; proc datasets library=csl nodetails nolist; modify csl1/ correctencoding='utf8'; quit; data fun; set csl.csl1; run; Here's my log: 1 libname csl 'C:\Users\rogersaj\Documents\CSL'; NOTE: Libref CSL was successfully assigned as follows: Engine: V9 Physical Name: C:\Users\rogersaj\Documents\CSL 2 proc datasets library=csl nodetails nolist; NOTE: Writing HTML Body file: sashtml.htm 3 modify csl1/ correctencoding='utf8'; WARNING: CORRECTENCODING was successful. MODIFY RUN group closed because the ENCODING value does not match the session encoding value. 4 quit; NOTE: PROCEDURE DATASETS used (Total process time): real time 0.50 seconds cpu time 0.34 seconds 5 data fun; 6 set csl.csl1; NOTE: Data file CSL.CSL1.DATA is in a format that is native to another host, or the file encoding does not match the session encoding. Cross Environment Data Access will be used, which might require additional CPU resources and might reduce performance. WARNING: Some character data was lost during transcoding in the dataset CSL.CSL1. Either the data contains characters that are not representable in the new encoding or truncation occurred during transcoding. 7 run; ERROR: Some character data was lost during transcoding in the dataset WORK.FUN. Either the data contains characters that are not representable in the new encoding or truncation occurred during transcoding. NOTE: The SAS System stopped processing this step because of errors. WARNING: The data set WORK.FUN may be incomplete. When this step was stopped there were 0 observations and 805 variables. NOTE: DATA statement used (Total process time): real time 0.05 seconds cpu time 0.04 seconds Please help? Thank you!!

rogersaj · ‎09-26-2017

I am running an analysis with a continuous outcome and a dichotomous outcome. Here's an outline of my table: I was able to get a p-value for trend in both the bivariate analyses (using PROC GLM> estimate; and PROC FREQ> trend) Now I need to determine if it is even possible to get p-value for trend in multivariate analyses. Theoretically, does that work? Here is my code: **For continuous variable; proc glm data = ....; class hivpos3 (ref = '0') marital3 (ref = '1') lmpyear (ref = '2011') agecat4 (ref = 'Aged 15 or less') distcat (ref = 'Distance less than or equal to 5 km'); model p1v1_gest = hivpos3 marital3 lmpyear agecat4 distcat /solution; estimate "Linear trend for lmpyear" lmpyear -1 0 1; estimate "Linear trend for agecat4" agecat4 -0.564165 -0.298675 0.0995585 0.7632818; estimate "Linear trend for distcat" distcat -0.447043 -0.305405 -0.084097 0.8365453; run; quit; **For dichotmous variable; proc logistic data = ....; class hivpos3 (ref = '0') marital3 (ref = '1') lmpyear (ref = '2011') agecat4 (ref = 'Aged 15 or less') distcat (ref = 'Distance less than or equal to 5 km')/param=ref; model returnee = hivpos3 marital3 lmpyear agecat4 distcat / link=logit stb; run; quit; For the continuous variable, I don't know if my "estimate" requests are valid. For the dichotmous variable, this guidance says that "In PROC LOGISTIC, the score test in the Testing Global Null Hypothesis: BETA=0 table is equivalent to the Cochran-Armitage test" so I am only getting one p-value. Is that the way it should be - just one global p-value? Thanks!

rogersaj · ‎09-26-2017

A few days ago, I asked for help interpeting a linear trend here: http://support.sas.com/kb/22/912.html It turns out my analysis was wrong because of my categories were unevenly spaced (in addition to not being sorted, which they are now). @sld pointed me to this resource to deal with unequally spaced treatments: http://support.sas.com/kb/22/912.html But I can't seem to figure out how to use PROC IML to find my coefficients for the linear trend. First, since my categories are distances, do I use the midpoint? Distance le 5km; midpoint = 2.5 Distance between 6 to 15km; midpoint = 10.5 Distance between 16 to 30km; midpoint = 23 Distance beween 31 to 100km; midpoint = 75 proc iml; distcat={2.5,10.5,23,75}; contrL=orpol(distcat,4); print contrL; quit; Output: If my code is right, how do I interpet the output to create my dose levels to test for linear trend with PROC GLM "Estimate"? Thanks!

rogersaj · ‎09-22-2017

Hi, Dependent variable = gestational age at first clinic visit (continuous) = p1v1_gest Independent variable = village distance from clinic = distcat I'd like to know if the gestational age at first clinic visit varies LINEARLY with village distance from clinic. Code: proc glm data = ...; class distcat; model p1v1_gest = distcat / solution; estimate "Linear trend for distcat" distcat -3 -1 1 3; contrast 'linear' distcat -3 -1 1 3; run; Results: (see attached) I have the following questions (If you can only answer one, the first is most important!): Can you please help me interpet the output with as regards the "Linear trend for distcat"? How does the "estimate" result differ from the "contrast" result in terms of interpetation? Was the use of -3 -1 1 3 appropriate or am I supposed to put in the median values and make them sum to zero? Once I add in covariates into the model, how will that change my interpretation? If, instead of distance, I had three marital status categories (Single, Currently Married, Previously Married) could I still test for a linear trend? If so, how would that interpretation work? Thanks SO much! AJ

rogersaj · ‎09-22-2016

I think this is fairly straightforward; my brain just can't solve this problem right now: A) I have a large SAS dataset of patients and it contains the villages that they each come from to get to clinic. B) I have a long list of distances from the village to the clinic. I need SAS to help me pull the distance from List B and put it into Dataset A. How do I do that? Thank you!

rogersaj · ‎08-04-2016

Thanks, I'll try those today andhopefully one of the strategies will work. Just as an FYI, I had trouble with the Link King because it required some elements that my database didn't have, and wasn't able to use more than one "flex" variable (ie. not one of the core identifier variables that it felt like I should be using) to identify my records, which would have been helpful in my circumstance. However, even when I was successful in getting most of my data to work within its parameters, it mysteriously deleted about 1/2 of my records. I spent half a day trouble shooting that, to no success.

rogersaj · ‎08-03-2016

I'm working with fuzzy data and trying to link prenatal care records within a dataset. See my original question post with the link below, where unfortunately none of the proposed solutions worked for me. My latest attempt is to use a working paper on "Linking records matching on at least m of n identifying keys" (http://www2.sas.com/proceedings/sugi27/p075-27.pdf). Unfortunately, I have yet to learn SQL so I'm sort of blindly submitting these statements. If someone could see what my error might be, I'd appreciate it. It creates a table with nothing in it. https://communities.sas.com/t5/Base-SAS-Programming/Assigning-subject-identifiers-to-patients-with-misspelled-names/m-p/288695 proc sql; create table pairs as select datasetX.subjectid as subjectid1, datasetY.subjectid as subjectid2 from datasetX datasetY where ( (datasetX.name1 eq datasetY.name1) + (datasetX.name2 eq datasetY.name2) + (datasetX.name3 eq datasetY.name3) + (datasetX.id eq datasetY.id) + (datasetX.dob eq datasetY.dob) + (datasetX.village eq datasetY.village) ) >= 2 and (datasetX.subjectid < datasetY.subjectid) order by subjectid1, subjectid2; quit; Thanks so much.

rogersaj · ‎08-01-2016

I'm working with a prenatal care dataset from a developing country. The data were abstracted from paper charts. Patient visits were recorded on a visit-by-visit basis, so there's no "patient file" per se. In this country, it's okay if names are spelled slightly differently, as long as it's in the ballpark phonetically. That means I now have the task of trying to link all these patient prenatal care records longitudinally while not having a consistent identifier. Here's what I do have: For identifiers, I have their first names, middle names, last names, village, age (in years, no birthdate), last menstrual period date (LMP), expected delivery date (EDD), and parity. The problem is that no one identifier is consistently right. How do I sort these patients out and assign them a subject ID? Thanks so much for your advice.

rogersaj · ‎07-20-2016

Reeza, you're my hero! It worked after I removed the equal sign and quotation marks.

rogersaj · ‎07-20-2016

I'm sure there's a simple answer to this, but how do I specify several numerical values that I'm interested in? For example, for variable "hivstatus" I want to select values 0,1,2 to mark those as "Known HIV status." I tried this but it didn't work: if hivstatus = in('0','1','2') then status = 1; Similarly, how do I select the opposite (ie, NOT 0,1,2)? Thank you!

Online Status	Offline
Date Last Visited	‎12-15-2022 08:45 AM

Polynomial regression: Fitting the curve

Re: Getting coefficient of variation and the power in the Box-Cox tran...

Getting coefficient of variation and the power in the Box-Cox transfor...

Re: Data file is in a format that is native to another host

Re: Data file is in a format that is native to another host

Data file is in a format that is native to another host

P-value for trend in multivariate analyses?

Choosing linear trend coefficients

Help interpreting linear trend

Pulling a value from a list and matching it to an observation

Re: Choosing linear trend coefficients

Re: Help interpreting linear trend

Re: Do Loop with multiple values

Re: Why this date variables are not converted into date9. format using...

Re: Count a value within group to assign new identifier

Help interpreting linear trend

Polynomial regression: Fitting the curve

Re: Getting coefficient of variation and the power in the Box-Cox tran...

Getting coefficient of variation and the power in the Box-Cox transfor...

Re: Data file is in a format that is native to another host

Re: Data file is in a format that is native to another host

Data file is in a format that is native to another host

P-value for trend in multivariate analyses?

Choosing linear trend coefficients

Help interpreting linear trend

Pulling a value from a list and matching it to an observation

Re: Proc sql code to help link fuzzy records

Proc sql code to help link fuzzy records

Assigning subject identifiers to patients with misspelled names, villa...

Re: Selecting multiple numerical values

Selecting multiple numerical values