File PRIN.RM

Directory of image this file is from
This file as a plain text file

	SUBROUTINE PRIN1
        COMMON IBUFF(255),DUMI(85),MANU,IBLK,CP,CV,CF
        COMMON/BLK1/VOLN(250),VOLX(250),FLON(260),FLOX(260),PRSN(250)
        COMMON/BLK2/ACCN(250),ACCX(250),PRSX(250)
        COMMON/BLK3/INEG,IPOS,IDIM,ITMP,IDIM1,IPAS,ISTOP,IREC
        COMMON/BLK4/HOLD(30),K,X,Y,APOS,J
	COMMON/BLK5/AMAT(5,5),COEF(4),COV(5,5),COR(5,5),A21(4)
	COMMON/BLK6/SD(5),A12(4),CEPT,DET,EPS,SX(5),SXY(5),SY
	DIMENSION SE(4),FV(4)
	AM=X
	A22=APOS
	DF=AM-4.
	X=0.
	DO 20 L=1,4
	COEF(L)=0.
	DO 10 I=1,4
  10	COEF(L)=COEF(L)+AMAT(I,L)*A12(I)
	X=X+COEF(L)*A21(L)
  20	CONTINUE
	BDD=A22-X
	RSS=A22-BDD
	AMS=BDD/DF
	F=RMS/AMS
	R=SQRT(RSS/A22)
	SEE=SQRT(AMS)
	DO 30 I=1,4
	SE(I)=SEE*AMAT(I,I)**0.5
  30	FV(I)=(COEF(I)/SE(I))**2
	AN=AM-1.
	Q=4.
	X=AM-Q
	WRITE(3,40)AN,A22,Q,RSS,RMS,F,X,BDD,AMS
  40	FORMAT(/////////1X,T30,'A N A L Y S I S     O F     V A',
     1	' R I A N C E     S U M M A R Y     T A B L E'///
     2	1X,T10,'SOURCE',T35,'DF',T50,'SS',T65,'MS',T80,'F',
     3	//,1X,T10,'TOTAL',T28,E15.6,T43,E15.6,//,1X,T13,'REGR',
     4	'ESSION',T28,E15.6,T43,E15.6,T58,E15.6,T73,E15.6,//,
     5	1X,T13,'DEVIATIONS',T28,E15.6,T43,E15.6,T58,E15.6,//)
	X=R**2
	WRITE(3,50)R,X,SEE,(COEF(I),SE(I),FV(I),I=1,4)
  50	FORMAT(1X,T10,'MULTIPLE CORRELATION COEFFICIENTS',//,
     1	1X,T15,'R =',T18,E15.6,T39,'R2 =',T43,E15.6,T62,'SE-EST =',
     2	T70,E15.6,///,1X,T10,'MULTIPLE REGRESSION COEFFI',
     3	'CIENTS',//,1X,T15,'COEF',T38,'SE',T58,'F',/,
     4	4(/1X,T10,E15.6,T30,E15.6,T50,E15.6))
	RETURN
	END



Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with.  I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation. 

PDP-8 Home Page   PDP-8 Site Map   PDP-8 Site Search