yesterday
TomHsiung
Pyrite | Level 9
Member since
01-04-2018
- 142 Posts
- 52 Likes Given
- 4 Solutions
- 7 Likes Received
-
Latest posts by TomHsiung
Subject Views Posted 559 a week ago 228 a week ago 229 a week ago 283 a week ago 679 a week ago 680 a week ago 690 a week ago 702 a week ago 710 a week ago 725 a week ago -
Activity Feed for TomHsiung
- Liked Re: PROC Tabulate omit one observation but PROC Univariate not for Kathryn_SAS. a week ago
- Posted Re: PROC Tabulate omit one observation but PROC Univariate not on Statistical Procedures. a week ago
- Posted Re: Two table statements in PROC Tabulate due to different statistics and we want merge them into on on SAS Programming. a week ago
- Posted Re: Two table statements in PROC Tabulate due to different statistics and we want merge them into on on SAS Programming. a week ago
- Liked Re: Two table statements in PROC Tabulate due to different statistics and we want merge them into on for Cynthia_sas. a week ago
- Posted Two table statements in PROC Tabulate due to different statistics and we want merge them into one on SAS Programming. a week ago
- Posted Re: PROC Tabulate omit one observation but PROC Univariate not on Statistical Procedures. a week ago
- Posted Re: PROC Tabulate omit one observation but PROC Univariate not on Statistical Procedures. a week ago
- Posted Re: PROC Tabulate omit one observation but PROC Univariate not on Statistical Procedures. a week ago
- Posted Re: PROC Tabulate omit one observation but PROC Univariate not on Statistical Procedures. a week ago
- Posted Re: PROC Tabulate omit one observation but PROC Univariate not on Statistical Procedures. a week ago
- Posted PROC Tabulate omit one observation but PROC Univariate not on Statistical Procedures. a week ago
- Posted Re: Did SAS Viya for Learners update the rule for naming the column when importing from xlsx on SAS Programming. 2 weeks ago
- Posted Did SAS Viya for Learners update the rule for naming the column when importing from xlsx on SAS Programming. 2 weeks ago
- Liked Re: Can we do the approximately unbiased estimator for the probability density function of cont. var for Season. 04-07-2025 08:12 AM
- Posted Re: Can we do the approximately unbiased estimator for the probability density function of cont. var on Statistical Procedures. 04-07-2025 06:32 AM
- Liked Re: Can we do the approximately unbiased estimator for the probability density function of cont. var for FreelanceReinh. 04-07-2025 05:10 AM
- Liked Re: Can we do the approximately unbiased estimator for the probability density function of cont. var for Season. 04-07-2025 02:10 AM
- Posted Re: Can we do the approximately unbiased estimator for the probability density function of cont. var on Statistical Procedures. 04-06-2025 11:31 AM
- Liked Re: Can we do the approximately unbiased estimator for the probability density function of cont. var for Season. 04-06-2025 07:45 AM
-
Posts I Liked
Subject Likes Author Latest Post 3 1 1 3 1 -
My Liked Posts
Subject Likes Posted 1 04-05-2025 03:22 AM 1 04-04-2025 10:20 AM 1 04-04-2025 09:36 AM 2 04-23-2024 02:51 AM 2 04-19-2024 11:38 AM
04-04-2025
11:31 AM
Hi, Season. Thank you for your answer. Yep, the bandwidth determines the degree of the approximation of unbiased estimation. The smaller, the better, I think. Can we set it manually in PROC UNIVARIATE or PROC SGPLOT? Thx
... View more
04-04-2025
10:20 AM
1 Like
Ticket created and thank you, pal.
... View more
04-04-2025
09:36 AM
1 Like
Hello, guys
I am reading a textbooks of mathematical statistics and it is funny and precise. I wonder:
Is there any procedure to get f_hat(x) in SAS? Thanks
... View more
03-30-2025
04:22 AM
Dear, guys
Today a new publication on NEJM says their data is sharable and would be available on "a specialized SAS data platform". Therefore, I would like to download the data and make statistical inferences by myself. And, I don't where is the SAS data platform? Does it belong to SAS?
Thanks
Source: https://www.nejm.org/doi/suppl/10.1056/NEJMoa2501006/suppl_file/nejmoa2501006_data-sharing.pdf
... View more
03-09-2025
04:04 AM
Thanks for the feedback.
We think we have a two-step approach.
First, we use the procedure to transfer continuous to spl variables, e.g., from the variable of age to spl variables age1, age2, and age3. Later, the regression procedure will estimate the coefficients of age1, age2, and age3.
Next, we want to map each age with its corresponding spl variables age1, age2, and age3. If there are n observations for age, there should be n for age1, age2, and age3 too.
The second step is that we use the data of age1, age2, and age3, along with their coefficients to estimate the outcome (e.g., HRs) for each observation by hand. Finally, we plot the relationship of the variable age and the outcome via proc sgplot.
So far, we don't know how to get the age1, age2, and age3 for each observation for the proc phreg.
... View more
03-08-2025
06:02 AM
Thanks. Better to know each observation's spline format of the continuous variable.
Stata achieves this goal by (another dataset but it use the Cox model too):
mkspline ages = ageY, cubic nknots(4) displayknots
stcox ages* BSA genderC AF Hypertension CHF serumAlb akiC ddiC indicationC CYP2C9code VKORC1code AFandWarfarinHistory
list ageY ages* if ageY == 60
predictnl xb = _b[ages1]*(ages1 - 60) + _b[ages2]*(ages2 - 13.19037) + _b[ages3]*(ages3 - 0.3942672) + _b[CYP2C9code]*(1 - 0), ci(lo hi)
twoway (line uci ageY, sort lc(black black) lp(- -)) (line lci ageY, sort lp(dash_dot)) (line hr ageY, sort lc(black) lp(l)), scheme(sj) ytitle(Hazard ratio (CYP2C9)) xtitle(Age (years)) yline(1.00, lpattern(dot))
... View more
03-08-2025
01:41 AM
proc phreg data=work.d202_raw plot(cl) = surv;
class genderC AF Hypertension CHF akiC ddiC indicationC CYP2C9code VKORC1Code / desc;
effect spl = spline(serumAlb / naturalcubic basis = bspline knotmethod = percentilelist(25 50 75));
model timeToFirstInrAttain*firstInrAttain(0) = bsa ageY genderC AF Hypertension CHF
spl akiC ddiC indicationC CYP2C9code VKORC1Code / rl = both;
hazardratio AF;
hazardratio Hypertension;
run;
We know the method to integrate the spline effect of the variable serumAlb into the Cox model and we get the regression coefficients of the spl variables.
But we don't know how we plot the spline plot, which has a y-axis for hazard ratios and x-axis for the original variable of serumAlb.
I know there is an example but it's for another procedure (which outputs the spl variable of each observation to a new dataset) instead of proc phreg.
I'm looking forward to your reply and I appreciate your help.
... View more
03-07-2025
11:03 AM
Hello, Rick Thank you for your respond. We intend to plot the spline (HR vs a continuous variable). Not sure about the strategy but Stata can output each observation with the spline data and plot the graph. So, how about in SAS Viya by code? Thanks.
... View more
03-07-2025
08:44 AM
Can we have an equivalent output in PROC PHREG? Thanks.
... View more
12-02-2024
07:22 AM
Thanks for the feedback. OK, if R-side model does not have the Zu matrix, then what is the difference between the mixed effect model and a regular linear model?
PS: R-side model: Y = X*beta + error is the same to the regular linear model: Y = X*beta + error. What makes the difference?
... View more
11-30-2024
07:48 AM
Does R-side PROC GLIMMIX neglect the random effect (Z*u)?
PS: y = X*Beta + Z*u + error
... View more
11-07-2024
09:30 AM
Thanks for sharing. I wonder how do we get the variance-covariance matrix of regression coefficients? Thanks.
... View more
- Tags:
- .
10-30-2024
09:21 AM
Thanks for the suggestion, sir. But I think I must first understand the G- and R-side models first. I know understand the formula of a general linear mixed model but the textbook I read does not discuss this topic in detail.
... View more
10-24-2024
08:08 AM
Hello, @Ksharp
I spent a couple of weeks learning the general linear mixed model. The model is in the form of the formula in the screenshot.
We can use the maximum likelihood method or matrix of variance-covariance to solve the point and interval estimation of the regression coefficients. What I don't understand is: What is the difference between the G- and R-side models? Both belong to the general linear mixed model, so they have the same model formula.
... View more
- Tags:
- to
09-30-2024
10:52 AM
Thanks for the update and I have additional questions.
What is the function of this line?
*random intercept / subject = Patient_ID;
With this line un-asterisked, I had the error:
NOTE: An R-side variance component is confounded with the profiled variance.
NOTE: The GLIMMIX procedure is modeling the probability that inr_attain='1'.
NOTE: Did not converge.
NOTE: The GLIMMIX procedure deleted the model item store WORK.MIXEDMODEL because of incomplete information for a subsequent
analysis.
In addition, in proc mixed, the random and repeated statements are used to represent between- and within-individual differences. What about the proc glimmix?
... View more