03-01-2023
jacksonan123
Lapis Lazuli | Level 10
Member since
05-11-2015
- 283 Posts
- 8 Likes Given
- 4 Solutions
- 28 Likes Received
-
Latest posts by jacksonan123
Subject Views Posted 928 01-28-2023 02:55 PM 1004 01-28-2023 07:19 AM 377 03-20-2022 11:30 PM 1235 01-10-2021 09:49 AM 1302 01-09-2021 11:47 AM 1374 01-09-2021 09:32 AM 870 11-04-2020 05:21 PM 912 11-04-2020 02:36 PM 972 11-04-2020 05:17 AM 1059 11-03-2020 01:55 PM -
Activity Feed for jacksonan123
- Posted Re: Data Step- populate data from prior time to current time on SAS Programming. 01-28-2023 02:55 PM
- Posted Data Step- populate data from prior time to current time on SAS Programming. 01-28-2023 07:19 AM
- Got a Like for In file statement SAS studio. 03-20-2022 11:31 PM
- Posted In file statement SAS studio on SAS Programming. 03-20-2022 11:30 PM
- Liked Re: Sorting complex folder names for Ksharp. 01-11-2021 07:27 AM
- Posted Re: Sorting complex folder names on SAS Programming. 01-10-2021 09:49 AM
- Posted Re: Sorting complex folder names on SAS Programming. 01-09-2021 11:47 AM
- Liked Re: Sorting complex folder names for novinosrin. 01-09-2021 09:58 AM
- Posted Sorting complex folder names on SAS Programming. 01-09-2021 09:32 AM
- Posted Re: Replacing missing column variables on SAS Programming. 11-04-2020 05:21 PM
- Posted Replacing missing column variables on SAS Programming. 11-04-2020 02:36 PM
- Posted Re: SAS Macro on SAS Programming. 11-04-2020 05:17 AM
- Posted SAS Macro on SAS Programming. 11-03-2020 01:55 PM
- Posted Re: SAS Macro on SAS Programming. 10-31-2020 04:59 PM
- Posted Re: SAS Macro on SAS Programming. 10-30-2020 04:39 PM
- Posted Re: SAS Macro on SAS Programming. 10-30-2020 04:29 PM
- Posted Re: SAS Macro on SAS Programming. 10-30-2020 01:00 PM
- Posted Re: SAS Macro on SAS Programming. 10-30-2020 12:46 PM
- Posted Re: SAS Macro on SAS Programming. 10-30-2020 12:31 PM
- Posted Re: SAS Macro on SAS Programming. 10-30-2020 11:22 AM
-
Posts I Liked
Subject Likes Author Latest Post 3 1 1 1 1 -
My Liked Posts
Subject Likes Posted 1 03-20-2022 11:30 PM 2 07-23-2019 11:14 AM 1 08-01-2017 10:20 AM
08-26-2016
06:39 AM
I conducted a BOX-COX analysis using the following code: ods rtf file='/folders/myfolders/regr4BCoxNresuls/CLnrGFR1.rtf' style=journal; ods TRACE ON; ods output FitStatistics=FitStat ParameterEstimates=param; proc transreg data=boot2 plots=all; /*plots(only)=FITPLOT(stats=none); */ model boxcox(CLnr/convenient lambda=-3 to 3 by 0.125)=identity(GFR1)/cl; /*SPEC;*/ output out=pred ; proc print data=fitstat;run; proc print data=param;run; proc print data=pred;run; run; quit; ods rtf close; The best lamda was =0 which means that the log transformation for CLnr was best. The data is posted below (Nstudyr4.csv)set which if I use log transformation for CLnr would give negative values. On line it was suggested that I can do log transformation if I add a constant to each value. My question is how do I interpret the value for the intercept if I do a regression when I take the antilog with the added constant? I have attached the data set and the BOX-Cox output. Subject CL CLr CLnr GFR1 GFR2 1 2.02 1.73 0.29 125.5 124.1 2 2.18 1.63 0.55 114.3 113.9 3 1.95 1.28 0.67 102.2 102 4 1.83 1.16 0.67 101.6 100.7 5 1.64 1.47 0.17 96.6 95.6 6 1.94 1.72 0.22 90.7 89.2 7 1.98 1.47 0.51 87.8 86.2 8 1.67 1.07 0.6 84.5 83.6 9 1.25 1.19 0.06 84.2 83.4 10 0.9 0.74 0.16 77.7 76.5 11 1.53 1.06 0.47 75.6 75.2 12 1.52 0.82 0.7 68 67.5 13 1.58 1.23 0.35 66.5 65.3 14 0.93 0.75 0.18 64.4 64.2 15 0.93 0.95 . 60 59.6 17 0.61 0.48 1.2 52.3 52.5 18 0.76 0.46 0.15 52 52.3 19 0.57 0.44 0.32 51.7 52.3 20 0.66 0.4 0.17 49.3 49.7 21 0.51 0.15 0.51 37.5 37.7 22 0.41 0.32 0.19 35.5 35.7 23 0.65 0.57 . 34.3 34.4 24 0.9 0.58 0.07 31 31.4 25 0.5 0.29 0.61 28.1 29.3 26 0.31 0.22 0.28 26.6 26.8 27 0.3 0.21 0.1 24.8 24.7 28 0.48 0.38 . 18 18.6 29 0.42 0.19 0.29 16.3 16.6 30 0.31 0.15 0.27 16 16.6 31 0.2 0.12 0.19 12.7 12.7 32 0.18 0.2 . 12.7 12.7 33 0.25 0.18 . 12.7 12.7 34 0.35 0.13 0.12 10.3 12.7 35 0.2 0.12 0.23 9.8 12.7 36 0.15 0.08 0.12 8 11.5 37 0.15 0.21 . 7.7 10.2 38 0.11 0.092 0.058 6.5 10.2 39 0.046 0.092 0.018 4.7 8.9 40 0.076 0.061 . 4.1 6.4 41 0.14 0.034 0.042 4.1 3.8 42 0.16 0.018 0.122 4.1 3.8 43 0.18 0.015 0.145 1.8 1.5 44 0.15 0.018 0.162 1.5 1.5
... View more
08-15-2016
11:01 AM
I will try each of the suggested options and see which works best. Thanks for the advice.
... View more
08-14-2016
12:22 PM
I have data for which I did a regression and the White test for constant variance had a p value =0.0016 indicating heterosdedasticity of the variance. See attached graph for residuals vs predicted value. The data was normally distributed see atached distribution graph which had a p=0.79 for the Shapiro Wilk test. These results indicate that I need to do a weighted regression. proc reg; /* weighted linear regression */ model y = x; weight w; In the literature I read, " If however we know the noise variance σ 2 i at each measurement i, and set wi = 1/σ2 i , we get the heteroskedastic MLE, and recover efficiency." My question is how do we know this weight value and based upon my data what would be an appropriate weight i.e., 1/y or something else?
... View more
08-04-2016
02:25 PM
Where in the present code would I set page size and what is the correct syntax?
... View more
08-04-2016
01:21 PM
I have the following code which runs fine:
ods rtf file='/folders/myfolders/regr1resuls/CLCLr.rtf' style=journal;
ods TRACE ON;
ods output FitStatistics=FitStat ParameterEstimates=param;
proc reg data=boot1 plots=all;
/*plots(only)=FITPLOT(stats=none);*/
model CL=CLr /clb SPEC;
output out=pred predicted=p lcl=lcl lclm=lclm ucl=ucl uclm=uclm residual=r;
proc print data=fitstat;
run;
proc print data=param;
run;
proc print data=pred;
run;
run;
quit;
ods rtf close;
The problem is that the outputted graphs are truncated (i.e., the X axis is not visible). I have attached an output for your perusal.
Can someone tell me why this occurs and how to correct it?
... View more
08-02-2016
09:45 AM
This is my current code for which I get no errors but I also get no stat output in my Table CLuGfru. The output is only the ANOVA and the graphs as before. ods rtf file='/...../CLUGFRu.rtf' style=journal; ods trace on; ods output FitStatistics=FitStatistics OutputStatistics=OutputStatistics; proc reg data=boot1 plots=all; model CLu=GFru; run; quit; ods rtf close; I have attached my file boot1 so that you can run the data.
... View more
08-02-2016
09:23 AM
I put it in but got no output although it works with base SAS but I guess it may not be the proper syntax for University SAS. I have trace on for my ODS and there is an output called FitStatistics which I have tried to add to the code as: proc reg data =boot1 plots=all Fit Statistics; However I get a syntax error. Can someone tell me how to correct the statement to get the output of Fit Statistics and the plots?
... View more
08-02-2016
08:43 AM
I am running the following code as previously discussed: ods rtf file='/folders/myfolders/peter/sample1.rtf' style=journal; ods trace on; proc reg data=boot plots=all; model Clu=GFru; run; quit; ods rtf close; It ran okay and I got all of the plots and the analysis of variance results. My question is how does one get the observed and predicted values (slope, intercept etc) output as well as their 90% confidence intervals in a table?
... View more
08-01-2016
08:25 AM
I am using the following code within a macro to export files in SAS University and the resulting file is a word document. proc export data=normal&i outfile="folders/mufolders/regr/nnormal&i" dbms=csv replace;run; How can I edit this code to have the output be a *.csv file?
... View more
07-29-2016
11:53 AM
I will use SGplot which will give me more control similar to when one uses gplot in base SAS. Thanks for the help.
... View more
07-29-2016
11:37 AM
Yes, that worked quite well. One final question if I wanted to use the graph for a presentation is there any way for me to edit the output to remove legends etc if it were output as a word document instead of an rtf file?
... View more
07-29-2016
11:25 AM
With your suggested code I did get the names for the outputs. The only one that I really want is the regression with confidence intervals which is part of the FitPlot. Unfortunately I get residual plot, and fit diagnostics which I don't want outputted. Since all of these are part of the FitPlot is there any way to get just the desired plot for the regression?
... View more
07-29-2016
06:30 AM
The code used was : ods rtf file='/folders/myfolders/peter_regr/sample.rtf' style=journal; proc reg data=er plots=(all); model change=conc; run; quit; ods rtf close; This code ran okay. The question that remains is there a way to selectively output only selected graphs using a different syntax for the plots(all) command?
... View more
- « Previous
- Next »