About gretaolsson

gretaolsson · ‎05-10-2017

Hi, I want to execute Exact logistic regression on a sample with a continuous variable. The sample size is 100 and below I show an extract from the data: "","num","y","x" "1",1,1,0.242106882663965 "2",1,0,-1.0215795935404 "3",1,1,-0.166510155016807 "4",1,0,-0.979868182931502 "5",1,1,0.078228780249738 "6",1,0,-0.37946718709135 "7",1,1,2.22519466020747 "8",1,0,0.0303031993682403 "9",1,0,-1.82711480734735 "10",1,0,1.1355902983012 I use the following code: proc import out = filename file = "/folders/myfolders/sasuser.v94/Btest.csv" dbms = CSV replace; run; data want; set filename; sampleid=ceil(_n_/100); /* Assign sampleid as record blocks of 10 */ run; proc logistic data=want desc; by sampleid; model y = x; exact x / estimate = both; run; (The beginning of the code I have because later on I want to have a csv file with several samples to be analyzed.) And I get the following message: An exception was thrown while sending a packet to the peer. What am I doing wrong?

gretaolsson · ‎05-05-2017

I would like to save all the parameter estimates from each sample to a vector or list, is that possible?

gretaolsson · ‎05-05-2017

What do you mean with <dataset to store>? Should I create a new file where I can save, for example, estimates? How do I do that then?

gretaolsson · ‎05-04-2017

Hi! I changed the code and now it works! Thank you very much! PROC IMPORT OUT = filename FILE = "/folders/myshortcuts/sf_myfolder/filename.csv" DBMS = CSV REPLACE; RUN; data want; set filename; sampleid=ceil(_n_/10); /* Assign sampleid as record blocks of 10 */ run; proc logistic data=want; by sampleid; class x y; model y = x; run; Is it possible to extract specific values (estimates and p-values) from the result of each regression and put them together into a vector or something similar?

gretaolsson · ‎05-04-2017

Hi, This was really helpful, thankyou! I have some questions, I think Im missing something obvious. Now I have the following code: PROC IMPORT OUT = filename FILE = "/folders/myshortcuts/sf_myfolder/filename.csv" DBMS = CSV REPLACE; RUN; data want; set have; sampleid=ceil(_n_/10); /* Assign sampleid as record blocks of 10 */ run; proc logistic data=want; by sampleid; class x y; model y = x; run; Of course, this does not work and my log looks like this: ... NOTE: WORK.FILENAME data set was successfully created. NOTE: The data set WORK.FILENAME has 100 observations and 4 variables. NOTE: PROCEDURE IMPORT used (Total process time): real time 0.06 seconds cpu time 0.04 seconds 95 96 data want; 97 set have; ERROR: File WORK.HAVE.DATA does not exist. 98 sampleid=ceil(_n_/10); /* Assign sampleid as record blocks of 10 */ 99 run; NOTE: The SAS System stopped processing this step because of errors. WARNING: The data set WORK.WANT 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.00 seconds cpu time 0.01 seconds 100 NOTE: PROCEDURE LOGISTIC used (Total process time): real time 0.00 seconds cpu time 0.00 seconds ERROR: Variable X not found. ERROR: Variable Y not found. NOTE: The SAS System stopped processing this step because of errors. 101 proc logistic data=want; 102 by sampleid; 103 class x y; 104 model y = x; 105 run; 106 107 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 120 I thought I would replace "want" with "filename", but it did not work. How should I do?

gretaolsson · ‎05-04-2017

Okay, I did not really catch any question, sorry. I have a CSV file containing 100 samples. Each sample consists of 10 observations. (So overall I have 1000 observations in one file.) I now want to let SAS read the CSV file and perform logistic regression on the first sample, which is observation "1" to "10" in the CSV-file, then I want SAS to perform logistic regression on the second sample, which is observation "11" to "20" in the CSV-file, and so on and on. Above I have sent an excerpt on how the data looks in the CSV file. "1" - "10" is the first sample, "11" - "20" is the second sample, ...., "91" - "100" is the tenth sample and so on up to the hundred sample. If you have any more questions that you need answered, please let me know.

gretaolsson · ‎05-04-2017

How does SAS know how many observations it is in each sample? This did not work for me ....

gretaolsson · ‎05-03-2017

Hello, I have a CSV file containing 100 samples a 10 observations each. How do I perform logistic regression on each individual sample from the one csv-file? How do I extract the values I'm interested in to a vector, such as estimatates and p-values? proc logistic data = WORK.IMPORT; class x y; model y = x; run; Here is an extract from the CSV-file: "","num","y","x" "1",1,1,0 "2",1,1,1 "3",1,1,1 "4",1,1,1 "5",1,1,0 "6",1,1,0 "7",1,1,1 "8",1,1,1 "9",1,1,1 "10",1,1,1 "11",1,1,1 "12",1,1,1 "13",1,1,0 "14",1,0,0 "15",1,1,0 "16",1,1,0 "17",1,1,0 "18",1,1,1 "19",1,1,1 "20",1,1,0 "21",1,0,0 "22",1,1,1 "23",1,0,0 "24",1,1,1 "25",1,1,1 "26",1,1,1 "27",1,1,0 "28",1,1,1 "29",1,0,0 "30",1,1,0 "31",1,1,1 "32",1,1,0 "33",1,1,1 "34",1,1,0 "35",1,0,0 "36",1,1,0 "37",1,1,0 "38",1,0,1 "39",1,1,1 "40",1,1,1 "41",1,1,0 Grateful for all the help I can get!

gretaolsson · ‎04-13-2017

Hello, I want to analyze simulated data of different sizes in both R and SAS. I have simulated all the samples in R and where I've done all my analysis, now I want to analyze all of those samples in SAS. All samples are in a the data frame. I want to do logistic regression and exact logistic regression, and I'm interested in the estimated coefficients. (university edition) 1. How do I transfer all samples to SAS smoothly? (So that point 2 gets easy) 2. How do I make a loop so that I can analyze all the data directly? Note that the sample size goes from 5-1500 with an increment of 5 and for now I only have one explanatory variable. Plz help :)

gretaolsson · ‎04-12-2017

My Data: y x1 x2 1. 1 1.489611900786800 -0.486983894512530 2. 1 0.887638190472230 -0.899961461187430 3. 1 -0.328400349680380 0.320480850960210 4. 0 -1.283346136073470 0.314729922388780 5. 1 -0.014384666024895 -1.040793737862780 6. 1 1.005337941612940 0.385444205622100 7. 1 0.403112850999760 0.797554772638080 8. 1 1.432508077938930 -0.553701810045310 9. 1 0.137341139238340 0.177313212434980 10. 0 -1.341507064615280 0.042985039337917 Logg: 1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 61 62 PROC IMPORT DATAFILE=REFFILE 63 DBMS=DBF 64 OUT=WORK.IMPORT1; 65 RUN; NOTE: Import cancelled. Output dataset WORK.IMPORT1 already exists. Specify REPLACE option to overwrite it. NOTE: The SAS System stopped processing this step because of errors. NOTE: PROCEDURE IMPORT used (Total process time): real time 0.00 seconds cpu time 0.00 seconds 66 67 68 proc logistic data = WORK.IMPORT1; 69 model y = x1 x2; 70 run; NOTE: PROC LOGISTIC is modeling the probability that y=0. One way to change this to model the probability that y=1 is to specify the response variable option EVENT='1'. WARNING: There is a complete separation of data points. The maximum likelihood estimate does not exist. WARNING: The LOGISTIC procedure continues in spite of the above warning. Results shown are based on the last maximum likelihood iteration. Validity of the model fit is questionable. NOTE: There were 10 observations read from the data set WORK.IMPORT1. NOTE: PROCEDURE LOGISTIC used (Total process time): real time 0.12 seconds cpu time 0.12 seconds 71 72 proc logistic data = WORK.IMPORT1 desc; 73 model y = x1 x2; 74 exact x1 x2 /estimate=both; 75 run; NOTE: PROC LOGISTIC is modeling the probability that y=1. WARNING: There is a complete separation of data points. The maximum likelihood estimate does not exist. WARNING: The LOGISTIC procedure continues in spite of the above warning. Results shown are based on the last maximum likelihood iteration. Validity of the model fit is questionable. NOTE: There were 10 observations read from the data set WORK.IMPORT1. NOTE: PROCEDURE LOGISTIC used (Total process time): real time 0.12 seconds cpu time 0.12 seconds 76 77 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 90

gretaolsson · ‎04-12-2017

Do I need to give more information to get any help?

gretaolsson · ‎04-11-2017

Hi, I want to perform exact logistic regression in SAS. I've found the following code that I want to apply to different samples of varying size. (I use the university edition. ) PROC IMPORT DATAFILE=REFFILE DBMS=DBF OUT=WORK.IMPORT; RUN; proc logistic data = WORK.IMPORT desc; model y = x1 x2; exact x1 x2 / estimate = both; run; When I run this code I get empty tables with no estimates... Must the data be written in a specific way, in that case, how? I can perform ordinary logistic regression on the samples, and my goal is to compare the results. I have attached the three files, log, results and data - that contains 20 observations. Because the files did not have the valid extension they are all in paint, sorry for that. I'm grateful for all the help I can get.

gretaolsson · ‎04-07-2017

I am completely new to SAS. You can see my code below. When I run the last part, I get the error message: 94 95 96 Proc sgplot data = pred; 97 plot p * x_1 = x_2 / vaxis = Axis1 haxis = Axis2; ____ 180 ERROR 180-322: Statement is not valid or it is used out of proper order. Why is that? How do I solve it? I also wonder how I can run the same code for several different samples, lets say 300 different samples. Can I make a loop for this and how do I do that? (SAS® Studio university edition) PROC IMPORT DATAFILE=REFFILE DBMS=DBF OUT=WORK.IMPORT1; RUN; proc freq data = WORK.IMPORT1; tables x_1*(x_2 ytest); tables x_2*ytest; weight num; run; proc tabulate data = WORK.IMPORT1; class x_1 x_2 ytest; tables x_1='x_1', ytest*x_2='x_2'*F=6. / rts=13.; freq num; run; proc logistic data = WORK.IMPORT1 desc; class x_1 x_2 ytest; freq num; model ytest = x_1 x_2; exact x_1 x_2 / estimate = both; run; proc logistic data = WORK.IMPORT1 desc; class x_1 x_2 ytest; freq num; model ytest = x_1 x_2; exact x_1 x_2 / estimate = both; output out = pred predicted = p; run; symbol1 c=blue v=circle i=join; symbol2 c=red v=plus i=join; symbol3 c=black v=square i=join; axis1 label=(r=0 a=90) minor=none; axis2 minor=none order=(0 1); proc sgplot data= pred; plot p*x_1=x_2 / vaxis=axis1 haxis=axis2; run; quit;

Online Status	Offline
Date Last Visited	‎05-15-2017 03:24 PM

Exact logistic regression continous variable

Re: Do loop

Re: Do loop

Re: Do loop

Re: Do loop

Re: Do loop

Re: Do loop

Do loop

Transfer data and loops

Re: Exact logistic regression

Re: Exact logistic regression continous variable

Exact logistic regression continous variable

Re: Do loop

Re: Do loop

Re: Do loop

Re: Do loop

Re: Do loop

Re: Do loop

Do loop

Transfer data and loops

Re: Exact logistic regression

Re: Exact logistic regression

Exact logistic regression

How do I solve Error: Statement is not valid or it is used out of prop...