Statistics

sas_file1.docx

Home >Mathematics homework help >Statistics

Study about the effect of fungi, cultivars, and countries on seed germination and number of days it took the seeds from plant them until the shooting system appear in the soil surface.

All of that was done by SAS

Table of Country by Cultivar
Country(Country)	Cultivar(Cultivar)
	ba	ha	th	Total
sa	100	100	100	300
us	100	100	100	300
Total	200	200	200	600

Fungi
Fungi	Frequency	Percent
0	238	39.67
1	362	60.33

Germination
Germination	Frequency	Percent
0	231	38.50
1	369	61.50

Table of Fungi by Germination
Fungi(Fungi)	Germination(Germination)
	0	1	Total
0	75	163	238
1	156	206	362
Total	231	369	600

Type
Type	Frequency	Percent
Alternaria	5	0.83
Alternaria+Aspergillus	1	0.17
Alternaria+Penicillium	1	0.17
Aspergillus	137	22.83
Aspergillus+Chaetomium	5	0.83
Chaetomium	53	8.83
None	238	39.67
Penicillium	105	17.50
Penicillium+Aspergillus	12	2.00
Penicillium+Chaetomium	15	2.50
Ulocladium	26	4.33
Ulocladium+Penicillium	1	0.17
Unknow	1	0.17

We have two countries: sa and us

3 Cultivar: ba, ha, th

Fungi: 0 means absent of fungi, 1 means present of fungi

Here we tested the germination

GERMINATION

Model Information
Data Set	WORK.GERMINATION
Response Variable (Events)	COUNT
Response Variable (Trials)	total
Response Distribution	Binomial
Link Function	Logit
Variance Function	Default
Variance Matrix	Diagonal
Estimation Technique	Maximum Likelihood
Degrees of Freedom Method	Residual

Class Level Information
Class	Levels	Values
Country	2	sa us
Cultivar	3	ba ha th
Fungi	2	0 1

Type III Tests of Fixed Effects
Effect	Num DF	Den DF	F Value	Pr > F
Country	1	146	14.40	0.0002
Cultivar	2	146	26.83	<.0001
Fungi	1	146	3.01	0.0851
*CountryCultivar**	2	146	4.39	0.0141
*CountryFungi**	1	146	14.79	0.0002
*CultivarFungi**	2	146	22.14	<.0001

*CultivarFungi Least Squares Means**
Cultivar	Fungi	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean
ba	0	-5.1080	0.2532	146	-20.17	<.0001	0.05	-5.6084	-4.6077	0.006012	0.001513
ba	1	-4.0225	0.1586	146	-25.37	<.0001	0.05	-4.3358	-3.7091	0.01759	0.002741
ha	0	-3.4502	0.1136	146	-30.38	<.0001	0.05	-3.6747	-3.2257	0.03076	0.003387
ha	1	-4.4022	0.1797	146	-24.50	<.0001	0.05	-4.7574	-4.0470	0.01210	0.002149
th	0	-3.6548	0.1258	146	-29.06	<.0001	0.05	-3.9034	-3.4063	0.02521	0.003091
th	1	-3.0993	0.09743	146	-31.81	<.0001	0.05	-3.2919	-2.9067	0.04314	0.004021

*CountryFungi Least Squares Means**
Country	Fungi	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	0	-4.0783	0.1378	146	-29.60	<.0001	0.05	-4.3505	-3.8060	0.01665	0.002256	0.01274	0.02175
sa	1	-4.3071	0.1446	146	-29.79	<.0001	0.05	-4.5928	-4.0213	0.01329	0.001896	0.01002	0.01761
us	0	-4.0638	0.1316	146	-30.88	<.0001	0.05	-4.3239	-3.8036	0.01689	0.002186	0.01308	0.02180
us	1	-3.3756	0.09228	146	-36.58	<.0001	0.05	-3.5580	-3.1932	0.03307	0.002951	0.02771	0.03942

Fungi Least Squares Means
Fungi	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
0	-4.0710	0.1013	146	-40.20	<.0001	0.05	-4.2711	-3.8709	0.01677	0.001670	0.01377	0.02041
1	-3.8413	0.08793	146	-43.69	<.0001	0.05	-4.0151	-3.6676	0.02101	0.001809	0.01772	0.02490

Country Least Squares Means
Country	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	-4.1927	0.1029	146	-40.76	<.0001	0.05	-4.3960	-3.9894	0.01488	0.001508	0.01218	0.01818
us	-3.7197	0.07996	146	-46.52	<.0001	0.05	-3.8777	-3.5617	0.02367	0.001848	0.02028	0.02761

Cultivar Least Squares Means
Cultivar	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
ba	-4.5653	0.1527	146	-29.90	<.0001	0.05	-4.8670	-4.2635	0.01030	0.001556	0.007637	0.01388
ha	-3.9262	0.1064	146	-36.89	<.0001	0.05	-4.1365	-3.7159	0.01934	0.002018	0.01573	0.02376
th	-3.3771	0.07918	146	-42.65	<.0001	0.05	-3.5335	-3.2206	0.03302	0.002528	0.02837	0.03840

And here we tested the days

Days

Model Information
Data Set	WORK.DAYS
Response Variable	Days
Response Distribution	Poisson
Link Function	Log
Variance Function	Default
Variance Matrix	Diagonal
Estimation Technique	Maximum Likelihood
Degrees of Freedom Method	Residual

Class Level Information
Class	Levels	Values
Country	2	sa us
Cultivar	3	ba ha th
Fungi	2	0 1

Type III Tests of Fixed Effects
Effect	Num DF	Den DF	F Value	Pr > F
Country	1	359	8.18	0.0045
Cultivar	2	359	1.50	0.2236
Fungi	1	359	14.08	0.0002
*CountryCultivar**	2	359	3.06	0.0482
*CountryFungi**	1	359	1.28	0.2592
*CultivarFungi**	2	359	1.92	0.1487
*CultivarFungi Least Squares Means**
Cultivar	Fungi	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
ba	0	4.1483	0.03154	359	131.54	<.0001	0.05	4.0863	4.2103	63.3245	1.9971	59.5164	67.3763
ba	1	4.2209	0.02054	359	205.50	<.0001	0.05	4.1805	4.2613	68.0943	1.3987	65.3985	70.9012
ha	0	4.1254	0.01669	359	247.20	<.0001	0.05	4.0926	4.1583	61.8948	1.0329	59.8964	63.9599
ha	1	4.1591	0.02688	359	154.75	<.0001	0.05	4.1062	4.2119	64.0132	1.7204	60.7177	67.4876
th	0	4.1094	0.02130	359	192.96	<.0001	0.05	4.0675	4.1513	60.9084	1.2971	58.4101	63.5134
th	1	4.2219	0.01311	359	322.05	<.0001	0.05	4.1961	4.2477	68.1649	0.8936	66.4300	69.9451

*CountryFungi Least Squares Means**
Country	Fungi	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	0	4.1420	0.01717	359	241.22	<.0001	0.05	4.1082	4.1757	62.9272	1.0805	60.8377	65.0885
sa	1	4.2367	0.01650	359	256.72	<.0001	0.05	4.2043	4.2692	69.1811	1.1417	66.9718	71.4632
us	0	4.1134	0.02162	359	190.28	<.0001	0.05	4.0709	4.1559	61.1548	1.3220	58.6094	63.8107
us	1	4.1645	0.01710	359	243.60	<.0001	0.05	4.1309	4.1982	64.3635	1.1003	62.2356	66.5643

Fungi Least Squares Means
Fungi	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
0	4.1277	0.01379	359	299.24	<.0001	0.05	4.1006	4.1548	62.0347	0.8557	60.3745	63.7405
1	4.2006	0.01149	359	365.71	<.0001	0.05	4.1780	4.2232	66.7289	0.7665	65.2384	68.2533

Country Least Squares Means
Country	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	4.1894	0.01226	359	341.77	<.0001	0.05	4.1652	4.2135	65.9801	0.8088	64.4086	67.5899
us	4.1390	0.01175	359	352.16	<.0001	0.05	4.1159	4.1621	62.7386	0.7374	61.3052	64.2056

Cultivar Least Squares Means
Cultivar	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
ba	4.1846	0.01849	359	226.31	<.0001	0.05	4.1482	4.2209	65.6661	1.2142	63.3212	68.0979
ha	4.1423	0.01460	359	283.65	<.0001	0.05	4.1135	4.1710	62.9451	0.9192	61.1631	64.7790
th	4.1656	0.01193	359	349.15	<.0001	0.05	4.1422	4.1891	64.4346	0.7688	62.9403	65.9642

Here we tested each fungus separately

Aspergillus

Model Information
Data Set	WORK.DAYS
Response Variable	Days
Response Distribution	Poisson
Link Function	Log
Variance Function	Default
Variance Matrix	Diagonal
Estimation Technique	Maximum Likelihood
Degrees of Freedom Method	Residual

Class Level Information
Class	Levels	Values
Country	2	sa us
Cultivar	3	ba ha th
Fungi	1	1

Fit Statistics
-2 Log Likelihood	219.37
AIC (smaller is better)	231.37
AICC (smaller is better)	234.87
BIC (smaller is better)	239.98
CAIC (smaller is better)	245.98
HQIC (smaller is better)	234.18
Pearson Chi-Square	29.86
Pearson Chi-Square / DF	1.19

Type III Tests of Fixed Effects
Effect	Num DF	Den DF	F Value	Pr > F
Country	1	25	1.78	0.1940
Cultivar	2	25	2.04	0.1513
*CountryCultivar**	2	25	2.10	0.1436

*CountryCultivar Least Squares Means**
Country	Cultivar	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	ba	4.2965	0.03890	25	110.46	<.0001	0.05	4.2164	4.3766	73.4444	2.8567	67.7905	79.5699
sa	ha	4.2691	0.03567	25	119.69	<.0001	0.05	4.1956	4.3425	71.4545	2.5487	66.3936	76.9013
sa	th	4.2836	0.05872	25	72.95	<.0001	0.05	4.1626	4.4045	72.5000	4.2573	64.2413	81.8204
us	ba	4.1537	0.07236	25	57.40	<.0001	0.05	4.0046	4.3027	63.6667	4.6068	54.8520	73.8979
us	ha	4.0943	0.1291	25	31.71	<.0001	0.05	3.8285	4.3602	60.0000	7.7460	45.9916	78.2751
us	th	4.3610	0.06523	25	66.85	<.0001	0.05	4.2266	4.4953	78.3333	5.1099	68.4856	89.5971

Cultivar Least Squares Means
Cultivar	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
ba	4.2251	0.04107	25	102.86	<.0001	0.05	4.1405	4.3097	68.3810	2.8087	62.8343	74.4174
ha	4.1817	0.06697	25	62.44	<.0001	0.05	4.0438	4.3196	65.4773	4.3849	57.0415	75.1606
th	4.3223	0.04389	25	98.49	<.0001	0.05	4.2319	4.4127	75.3602	3.3072	68.8477	82.4888

Country Least Squares Means
Country	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	4.2831	0.02632	25	162.75	<.0001	0.05	4.2289	4.3373	72.4618	1.9070	68.6388	76.4977
us	4.2030	0.05391	25	77.96	<.0001	0.05	4.0920	4.3140	66.8862	3.6059	59.8572	74.7407

Model Information
Data Set	WORK.DAYS
Response Variable	Days
Response Distribution	Poisson
Link Function	Log
Variance Function	Default
Variance Matrix	Diagonal
Estimation Technique	Maximum Likelihood
Degrees of Freedom Method	Residual

Chaetomium

Class Level Information
Class	Levels	Values
Country	2	sa us
Cultivar	3	ba ha th
Fungi	1	1

Fit Statistics
-2 Log Likelihood	299.61
AIC (smaller is better)	309.61
AICC (smaller is better)	311.19
BIC (smaller is better)	318.54
CAIC (smaller is better)	323.54
HQIC (smaller is better)	312.92
Pearson Chi-Square	32.80
Pearson Chi-Square / DF	0.84

Type III Tests of Fixed Effects
Effect	Num DF	Den DF	F Value	Pr > F
Country	1	39	2.49	0.1225
Cultivar	2	39	1.53	0.2289
*CountryCultivar**	1	39	1.16	0.2881

*CountryCultivar Least Squares Means**
Country	Cultivar	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	ha	4.1301	0.03824	39	108.02	<.0001	0.05	4.0527	4.2074	62.1818	2.3776	57.5540	67.1818
sa	th	4.3277	0.03186	39	135.82	<.0001	0.05	4.2632	4.3921	75.7692	2.4142	71.0401	80.8132
us	ba	4.2370	0.05376	39	78.81	<.0001	0.05	4.1283	4.3457	69.2000	3.7202	62.0698	77.1492
us	ha	4.0943	0.1291	39	31.71	<.0001	0.05	3.8332	4.3555	60.0000	7.7460	46.2109	77.9037
us	th	4.1386	0.03375	39	122.63	<.0001	0.05	4.0703	4.2069	62.7143	2.1165	58.5761	67.1448

Cultivar Least Squares Means
Cultivar	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
ba	Non-est	.	.	.	.	.	.	.	Non-est	.	.	.
ha	4.1122	0.06732	39	61.08	<.0001	0.05	3.9760	4.2484	61.0812	4.1121	53.3052	69.9915
th	4.2331	0.02321	39	182.41	<.0001	0.05	4.1862	4.2801	68.9334	1.5997	65.7724	72.2463

Country Least Squares Means
Country	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	Non-est	.	.	.	.	.	.	.	Non-est	.	.	.
us	4.1566	0.04795	39	86.68	<.0001	0.05	4.0596	4.2536	63.8569	3.0622	57.9540	70.3611

Model Information
Data Set	WORK.DAYS
Response Variable	Days
Response Distribution	Poisson
Link Function	Log
Variance Function	Default
Variance Matrix	Diagonal
Estimation Technique	Maximum Likelihood
Degrees of Freedom Method	Residual

None

Class Level Information
Class	Levels	Values
Country	2	sa us
Cultivar	3	ba ha th
Fungi	2	0 1

Fit Statistics
-2 Log Likelihood	1110.17
AIC (smaller is better)	1122.17
AICC (smaller is better)	1122.71
BIC (smaller is better)	1140.74
CAIC (smaller is better)	1146.74
HQIC (smaller is better)	1129.71
Pearson Chi-Square	136.20
Pearson Chi-Square / DF	0.87

Type III Tests of Fixed Effects
Effect	Num DF	Den DF	F Value	Pr > F
Country	1	157	0.61	0.4370
Cultivar	2	157	1.03	0.3602
*CountryCultivar**	2	157	2.22	0.1116

*CountryCultivar Least Squares Means**
Country	Cultivar	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	ba	4.1589	0.04167	157	99.81	<.0001	0.05	4.0766	4.2412	64.0000	2.6667	58.9437	69.4900
sa	ha	4.1100	0.02939	157	139.86	<.0001	0.05	4.0520	4.1681	60.9474	1.7910	57.5105	64.5897
sa	th	4.1418	0.01655	157	250.19	<.0001	0.05	4.1091	4.1745	62.9138	1.0415	60.8899	65.0050
us	ba	4.1454	0.04757	157	87.15	<.0001	0.05	4.0514	4.2393	63.1429	3.0034	57.4807	69.3627
us	ha	4.1457	0.01598	157	259.43	<.0001	0.05	4.1141	4.1773	63.1613	1.0093	61.1988	65.1867
us	th	4.0518	0.04663	157	86.90	<.0001	0.05	3.9597	4.1439	57.5000	2.6810	52.4411	63.0469

Cultivar Least Squares Means
Cultivar	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
ba	4.1521	0.03162	157	131.33	<.0001	0.05	4.0897	4.2146	63.5700	2.0099	59.7215	67.6665
ha	4.1279	0.01673	157	246.81	<.0001	0.05	4.0948	4.1609	62.0445	1.0377	60.0283	64.1283
th	4.0968	0.02474	157	165.60	<.0001	0.05	4.0479	4.1456	60.1460	1.4879	57.2777	63.1579

Country Least Squares Means
Country	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	4.1369	0.01787	157	231.51	<.0001	0.05	4.1016	4.1722	62.6076	1.1187	60.4364	64.8568
us	4.1143	0.02283	157	180.20	<.0001	0.05	4.0692	4.1594	61.2089	1.3975	58.5098	64.0324

Penicillium

Model Information
Data Set	WORK.DAYS
Response Variable	Days
Response Distribution	Poisson
Link Function	Log
Variance Function	Default
Variance Matrix	Diagonal
Estimation Technique	Maximum Likelihood
Degrees of Freedom Method	Residual

Class Level Information
Class	Levels	Values
Country	2	sa us
Cultivar	3	ba ha th
Fungi	2	0 1

Fit Statistics
-2 Log Likelihood	640.80
AIC (smaller is better)	650.80
AICC (smaller is better)	651.55
BIC (smaller is better)	663.13
CAIC (smaller is better)	668.13
HQIC (smaller is better)	655.77
Pearson Chi-Square	120.36
Pearson Chi-Square / DF	1.47

Type III Tests of Fixed Effects
Effect	Num DF	Den DF	F Value	Pr > F
Country	1	82	2.62	0.1093
Cultivar	2	82	0.50	0.6082
*CountryCultivar**	1	82	2.48	0.1191

*CountryCultivar Least Squares Means**
Country	Cultivar	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	ha	4.0456	0.05000	82	80.91	<.0001	0.05	3.9461	4.1450	57.1429	2.8571	51.7326	63.1189
sa	th	4.1554	0.04174	82	99.56	<.0001	0.05	4.0724	4.2384	63.7778	2.6620	58.6960	69.2995
us	ba	4.1603	0.03766	82	110.46	<.0001	0.05	4.0854	4.2352	64.0909	2.4138	59.4646	69.0772
us	ha	4.2195	0.08575	82	49.21	<.0001	0.05	4.0489	4.3901	68.0000	5.8310	57.3358	80.6477
us	th	4.1578	0.01642	82	253.18	<.0001	0.05	4.1251	4.1905	63.9310	1.0499	61.8762	66.0541

Cultivar Least Squares Means
Cultivar	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
ba	Non-est	.	.	.	.	.	.	.	Non-est	.	.	.
ha	4.1325	0.04963	82	83.27	<.0001	0.05	4.0338	4.2313	62.3355	3.0938	56.4751	68.8041
th	4.1566	0.02243	82	185.34	<.0001	0.05	4.1120	4.2012	63.8544	1.4320	61.0682	66.7677

Country Least Squares Means
Country	Estimate	Standard Error	DF	t Value	Pr > \|t\|	Alpha	Lower	Upper	Mean	Standard Error Mean	Lower Mean	Upper Mean
sa	Non-est	.	.	.	.	.	.	.	Non-est	.	.	.
us	4.1792	0.03169	82	131.86	<.0001	0.05	4.1162	4.2423	65.3139	2.0701	61.3229	69.5646

Here is the Survival Analysis by using R studio- for each fungus

Study

about the

effect of

fungi, cultivars, and countries on seed germination

and number of days it took

the seeds from plant them until the shooting

system

appear

in the soil surface.

of that was done by SAS

Table

Country

Cultivar

Country(Country)

Cultivar(Cultivar)

Total

100

300

100

300

Total

200

600

Fungi

Frequency

Percent

238

39.67

362

60.33

Germination

Frequency

Percent

231

38.50

369

61.50

Table

Fungi

Germination

Fungi(Fungi)

Germination(Germination)

Total

163

238