SASSTATHW1ANOVA (SAS-3task)
Описание файла
Файл "SASSTATHW1ANOVA" внутри архива находится в папке "STATHW1_Kazachuk". Документ из архива "SAS-3task", который расположен в категории "". Всё это находится в предмете "(ппп соиад) (sas) пакеты прикладных программ для статистической обработки и анализа данных" из 10 семестр (2 семестр магистратуры), которые можно найти в файловом архиве МГУ им. Ломоносова. Не смотря на прямую связь этого архива с МГУ им. Ломоносова, его также можно найти и в других разделах. .
Онлайн просмотр документа "SASSTATHW1ANOVA"
Текст из документа "SASSTATHW1ANOVA"
2
SAS/STAT. HOME WORK 1. ANOVA
| ФИО | Казачук Мария Андреевна |
ФОРМАТ ОТВЕТА
Архив STATHW1_<ваша фамилия>, содержащий три файла:
-
Этот документ с заполненными таблицами «ОТВЕТ», где каждая строка соответствует вопросу [P1], [P2] и т.д. (см. TASK 1. ниже).
-
Код, который вы использовали для выполнения заданий – один файл *.sas (!!! Добавляйте комментарии с номерами заданий и вопросов /* 1 [P1]*/, /* 2 [P3] */ и т.д. !!!)
-
Результат запуска кода в формате RTF (в окне редактирования кода нажать кнопку Properties (вверху справа) Results RTF, Graph format = JPEG – см. картинку ниже):
DATA
Рассмотрим эксперимент, в котором мы хотим протестировать три торговых марки (brands) бетона (concrete) и узнать делает ли добавка (additive) бетон каждой марки прочнее. Несколько участков земли было забетонировано. Следующие характеристики зафиксированы (набор данных CONCRETE_S3 в библиотеке libname mydata "/courses/u_cmc.msu.ru1/i_889205/c_6043/saslib" access=readonly; ):
Strength – прочность.
Additive – была ли использована добавка.
Brand – марка бетона, подвергшаяся тестированию.
TASKS
-
Используя процедуру MEANS, познакомиться с данными. [P1] Вывести средние в новый набор данных. [P2] Построить график процедурой SGPLOT, используя полученный набор данных: strength но оси Y, Additive по оси X, группировать по переменной Brand. [P3] Что вы можете сказать о данных? [P4] Основываясь на графике, нужно ли использовать пересечение факторов Additive и Brand в модели?
| ОТВЕТ [статистика=значение и/или график, ваш вывод] |
| [P1] |
| [P2] |
| [P3] Во всех трех случаях добавка (Reinforsed) существенно увеличивает прочность бетона. |
| [P4] Пересечение факторов Additive и Brand стоит учитывать, потому что в зависимости от марки бетона, мы получаем различные значения разницы прочности бетона с использованием и без использования добавки. |
-
[P1] Проверьте гипотезу о том, что средняя прочность одинакова для всех марок. Проверить предположения. Если возможно, сравните все марки с маркой Graystone. [P2] Добавьте оставшийся фактор – Additive. Какие выводы можно сделать сейчас? [P3] Если графики говорят из п.1 говорят, что нужно использовать пересечение, то добавьте его. Какие выводы вы можете сделать на данном шаге анализа.
| ОТВЕТ [код, статистика=значение и/или график, ваш вывод] |
| [P1] proc glm data=tmp1 plots(only)=diagnostics; class Brand; model Strength_Mean=Brand; means Brand/ hovtest welch; title 'Testing for Equality of Means with PROC GLM'; run; quit; Значение Pr > F равно 0.5184. Это значит, что гипотеза не верна, средняя прочность не одинакова для всех марок. У Graystone наибольшая прочность. |
| [P2] proc glm data=tmp1 plots(only)=diagnostics; class Brand Additive; model Strength_Mean=Brand Additive; title 'ANOVA for Randomized Block Design'; run; quit; Значение Pr > F равно 0.0028. Это значит, что гипотеза верна. Средняя прочность приблизительно равна для брендов с одинаковыми добавками. |
| [P3] proc glm data=tmp1 order=internal; class Brand Additive; model Strength_Mean=Brand Additive Brand*Additive; lsmeans Brand*Additive/ slice= Additive; title 'title1'; title2 'title2'; run; quit;
Графики не пересекаются, что нам и необходимо. |
-
Выполните подходящие множественные сравнения для статистически значимых переменных.
| ОТВЕТ [код, статистика=значение и/или график, ваш вывод] |
| proc glm data=tmp1 plots(only)=(controlplot diffplot(center)); class Brand Additive; model Strength_Mean=Brand Additive; lsmeans Brand / pdiff=all adjust=tukey; lsmeans Brand / pdiff=control('Graystone') adjust=dunnett; title 'Brand Data: Multiple Comparisons'; run; quit; На данных таблицах мы видим разницу между прочностями разных марок бетона (Least Squares Means). |
