[R] multi-regression with more than 50 independent variables

R DF rdfmodel at gmail.com
Mon Feb 13 16:19:06 CET 2012


Hi R Users,

I am going to run a multiple linear regression with around 57 independent
variables. Each time I run the model with just 11 variables, the results
are reasonable. With increasing the number of independent variables more
than 11, the coefficients will get “NA” in the output.  Is there any
limitation for the number of independent variables in multiple linear
regressions in R? I attached my dataset as well as R codes below:



mlr.data<- read.table("./multiple.txt",header=T)

mlr.output<- lm(formula = CaV ~ SHG +  TrD+  CrH+ SPAD+ FlN+ FrN+   YT+
LA+ LDMP+    B+Cu+  Zn+   Mn +   Fe+   K +  P+   N +Clay30 +Silt30 +Sand30
+Clay60 +Silt60 +Sand60 +ESP30 +NaEx30+ CEC30+Cl30+ SAR30 +KSol30+ NaSol30
+CaMgSol3 +ZnAv30 +FeAv30 +OC30 +PAv30 +KAv30 +TNV30+ pH30+ EC30 +SP30
+ESP60 +NaEx60 +CEC60  +Cl60 +SAR60 +KSol60 +NaSol60 +CaMgSol6
+ZnAv60+FeAv60 +OC60 +PAv60 +KAv60 +TNV60 +pH60 + EC60 +SP60, data=mlr.data)

summary (mlr.output)




Regards,

Reza
-------------- next part --------------
CaV	SHG	TrD	CrH	SPAD	FlN	FrN	YT	LA	LDMP	B	Cu	Zn	Mn	Fe	K	P	N	Clay30	Silt30	Sand30	Clay60	Silt60	Sand60	ESP30	NaEx30	CEC30	Cl30	SAR30	KSol30	NaSol30	CaMgSol3	ZnAv30	FeAv30	OC30	PAv30	KAv30	TNV30	pH30	EC30	SP30	ESP60	NaEx60	CEC60	Cl60	SAR60	KSol60	NaSol60	CaMgSol6	ZnAv60	FeAv60	OC60	PAv60	KAv60	TNV60	pH60	EC60	SP60
49.83	15.8	22.32	45.8	82.7	126	5	55.8	34.7	57.7	17.9	8.4	14.7	50.5	144.7	0.9	0.1	1.7	14	50	36	10	26	64	3.2	0.3	11.9	87.7	5.6	2.5	25.7	69.8	2.9	3.1	0.66	64.4	360	14.8	9.94	0.6	39.6	4.8	0.2	12.2	51.4	6.9	1.4	33.4	69.9	1.8	3.4	0.2	3.5	290	11.5	7.19	3.21	40.1
37.85	18.9	20.89	53.2	75.9	169	4	67.4	40.1	58.6	15.8	7.9	12.1	53.2	141.6	0.7	0.1	1.9	13	51	34	11	28	64.0	1.9	0.2	12.5	45.7	7.8	3.8	22.9	77.9	3.1	2.5	0.69	45.6	390	14.8	7.2	1.1	40.1	5.4	0.1	13.3	59.9	7.9	1.7	56.6	58.9	1.5	5.5	0.6	1.9	350	10.2	7.1	3.7	42.1
64.12	20.7	19.04	43.9	74.3	133	5	55.9	38.7	60.0	12.1	6.3	12.6	47.4	159.5	1.0	0.1	1.6	12	48	33	13	26	62.0	6.4	1.1	12.6	33.2	10.1	1.8	43.8	99.3	1.5	2.8	0.57	20.5	470	15.1	7.1	0.7	39.8	4.9	0.3	10.9	65.5	7.1	2.3	42.1	66.4	1.1	3.9	0.7	2.9	288.9	12.2	7.2	4.1	35.6
90.28	14.6	19.52	56.9	61.9	145	7	66.5	33.2	59.1	13.7	4.7	10.0	52.1	241.1	0.8	0.1	1.9	13	47	32	10	30	57.0	5.5	0.3	11.7	31.2	3.5	2.4	50.4	65.7	1.9	2.5	0.88	41.9	398	14.3	7	2.4	38.7	3.1	0.5	14.1	71.9	8.4	1.9	36.7	59.1	0.9	2.6	0.5	2.7	290.7	13.6	7	2.6	43.7
111.18	13.2	16.53	61.3	78	127	6	49.9	41.7	51.6	14.7	4.7	10.0	55.8	148.9	0.7	0.1	1.9	14	46	37	11	28	60.0	3.9	0.7	11	21.7	4.3	0.9	33.7	73.6	2.3	2.7	0.78	33.4	349	15.2	7.4	0.9	39.1	2.8	0.1	12.5	76.6	8.9	3.1	32.1	71.4	0.5	6.9	0.5	2.5	256.9	15.1	7.5	1.7	36.7
59.11	12.9	21.34	45.3	78.9	178	6	63.3	39.8	52.0	19.5	4.2	8.9	54.7	229.5	0.7	0.1	1.7	13	46	36	12	30	62.0	2.7	0.7	12.9	19.5	2.6	2.8	61.2	86.9	2.2	3.7	0.86	27.8	400.5	17.1	7.1	2.1	39.9	3.9	0.3	11.9	57.4	6.7	1.6	30.1	89.9	1.8	5.8	0.3	3.7	224.8	12.9	7.3	5.5	34.9
80.89	17.9	15.86	40.3	66.8	154	7	45.6	36.8	47.8	21.6	3.7	12.6	54.7	162.1	0.7	0.1	1.9	11	50	35	13	31	61.0	2.9	0.4	10.9	37.9	7.1	1.9	19.5	55.8	2.8	2.9	0.66	45.1	459	15.6	7.2	0.8	36.1	5.1	0.4	13.1	85.5	5.7	2.1	29.1	92.1	1.9	7.8	0.7	2.8	278.9	11.8	7.2	6.1	32.7
122.74	16.6	17.29	43	77.8	140	6	32.9	37.7	55.6	20.0	4.2	12.1	47.4	152.6	0.9	0.1	1.7	14	49	36	11	25	58.0	3.2	1.5	11.5	24.5	3.9	3.7	20.8	56.9	1.9	2.6	0.72	40.7	398	16.3	7.3	0.7	37.7	2.2	0.1	11.8	55.5	6.1	1.8	39.9	69.9	2.1	2.4	0.3	1.8	312.8	10.6	7.4	7.3	31.2
20.08	7.6	18.28	33.8	60.5	81	1	11.3	37.0	59.8	21.3	3.2	4.9	67.5	153.4	0.5	0.1	0.7	16	46	38	15	45	40	7.7	0.8	11.2	139.8	10.7	1.3	80.7	106.6	2.4	2.4	0.6	10.9	345.8	14.8	7.5	12.2	36.3	6.3	0.7	10.5	108.2	9.8	0.5	63.7	84.4	1.3	2.5	0.4	3.3	213.0	15.7	7.6	10.5	38.2
22.09	11.8	14.44	40.3	61.4	67	2	13.5	34.6	58.7	21.6	2.6	7.6	57.4	186.5	0.5	0.1	0.6	13	48	35	12	30	59	6.8	0.8	12.1	45.8	5.6	1.0	26.1	37.1	1.3	2.5	0.6	4.3	407.2	17.5	7.7	4.7	34.5	6.9	0.8	12.0	65.9	7.4	0.7	38.2	50.9	0.9	2.9	0.4	1.9	265.3	18.3	7.6	6.8	35.2
25.04	8.4	15.29	46.3	65.5	60	1	11.0	33.2	56.2	21.4	2.2	10.1	54.6	167.8	0.4	0.1	0.8	13	48	35	12	30	58	11.3	1.2	10.4	79.7	10.4	0.8	60.6	61.9	0.9	2.2	0.5	3.6	341.3	16.3	7.6	8.3	33.3	9.0	1.0	10.8	82.0	11.8	0.5	65.7	60.5	0.5	2.1	0.3	0.7	205.0	17.9	7.6	8.8	35.1
27.77	23.9	17.80	58.6	80.1	77	2	16.7	37.1	54.9	17.6	3.8	2.7	43.3	108.3	0.4	0.1	1.1	13	48	35	12	31	57	5.3	0.4	8.3	117.0	14.5	2.0	96.7	71.8	0.9	2.4	0.5	3.9	395.0	15.0	7.7	10.0	32.8	11.8	0.9	8.1	73.4	10.7	1.6	57.9	58.1	0.4	3.1	0.3	1.1	193.8	15.8	7.7	8.0	30.1
34.03	8.7	21.02	50.3	68.7	58	1	15.5	36.8	55.9	18.4	2.5	4.6	41.3	141.2	0.4	0.1	0.9	13	48	36	12	31	57	4.4	0.5	10.2	34.1	4.7	1.3	22.4	30.8	0.7	2.5	0.8	4.8	446.4	14.8	7.7	4.3	33.4	4.0	0.4	10.5	58.8	7.6	0.8	38.9	46.6	0.4	2.3	0.5	1.7	316.7	15.7	7.7	6.4	34.2


More information about the R-help mailing list