Version as of 2004-06-07 01:04:36
Clear message
회귀분석을 통한 터널 붕괴도 및 함수 추산 프로젝트.
- 기간 ; 2004-05-26 ~ 2004-06-23(수, 예정)
- 결정사항
- 단일 어플리케이션, c++ based(maybe)
- To Do
- spss, sas 등의 기존 프로그램을 통해서 회귀분석 가능 여부 타진
; 확인 후 해당 기능이 r-project에서 가능한지, 가능하다면 c++과의 인터페이스 작성은 어떻게 해야 하는지 고려.
- 기존에 공개된 통계 관련 (회귀분석이 가능한) 라이브러리 조사.
- ui 형태 결정
- RegressionAnalyzingGrapher_관련자료
- 알려진 사항
- Leverberg-Marquardt 회귀분석 방법 사용.
- modelling functions of tunnel convergence
x : 막장거리
X : 굴착면의 영향거리에 관련된 상수
Cx : 최종 내공변위
C : 내공변위
T : 암반의 시간 의존적 변형특성에 관련된 상수
Cx(1+m) : 굴착면이 무한히 멀고 장시간이 경과한 후의 최종 내공변위
(천층_광폭터널의_내공변위_및... 에서 발췌)
오전 1:03 2004-06-07, by Kenial
tunnel <- read.table("d:tunnel.txt", TRUE)
ft1 <- nls( c ~ A * ( 1 - exp( -B * x ) ) - c0, data = tunnel,
algorithm = "plinear",
start = list( A=0, B=0, c0=0 ), trace = TRUE)
ft1 <- nls( c ~ (( 1 - exp( -B * x ) ) * A ) - c0, data = tunnel,
start = list(A=0, B=0, c0=0 ), trace = TRUE)
ft1 <- nls( c ~ (( 1 - exp( -B * t ) ) * A ) - c0, data = tunnel,
start = list(A=0, B=0, c0=0 ), trace = TRUE)
ft1 <- nls( c ~ ( A * log(1 + (B*t)) ) - c0, data = tunnel,
start = list(A=0, B=0, c0=0 ), trace = TRUE)
ft1 <- nls( c ~ ( 1 - ( ( XX / ( XX + x ) ) ** 2 ) ) - c0, data = tunnel,
start = list(XX=0, c0=0 ), trace = TRUE)
ft1 <- nls( c ~ pa * (1 - exp(-pb * x)) + pc * (1 - exp(-pd * t)) - c0, data = tunnel,
start = list(pa=0, pb=0, pc=0, pd=0, c0=0 ), trace = TRUE)
ft1 <- nls( c ~ PCX*(1- ( (PX/(PX+x)) **2 )) * (1+ PM*(1-((PT/(PT+t))**0.3))) - c0,
data = tunnel, start = list(PCX=0, PX=0, PM=0, PT=0, c0=0 ), trace = TRUE)
오전 2:20 2004-06-06, by Kenial
* A * ( 1 - EXP(-B*x) ).
* NonLinear? Regression.
MODEL PROGRAM c0=0 A=0 B=0 .
COMPUTE PRED_ = A * ( 1 - EXP(-B*x) ) - c0.
NLR c
/OUTFILE='C:\SPSSFNLR.TMP'
/PRED PRED_
/CRITERIA SSCONVERGENCE 1E-8 PCON 1E-8 .
* A * ( 1 - EXP(-B*t) ).
* NonLinear? Regression.
MODEL PROGRAM c0=0 A=0 B=0 .
COMPUTE PRED_ = A * ( 1 - EXP(-B*t) ) - c0.
NLR c
/OUTFILE='C:\SPSSFNLR.TMP'
/PRED PRED_
/CRITERIA SSCONVERGENCE 1E-8 PCON 1E-8 .
* A * ln(1 + (B*t))
* NonLinear? Regression.
MODEL PROGRAM c0=0 A=0 B=0 .
COMPUTE PRED_ = A * ln(1 + (B*t)) - c0.
NLR c
/OUTFILE='C:\SPSSFNLR.TMP'
/PRED PRED_
/CRITERIA SSCONVERGENCE 1E-8 PCON 1E-8 .
* Cx * ( 1 - ( ( XX / ( XX + x ) ) ** 2 ) ).
* NonLinear? Regression.
MODEL PROGRAM c0=0 Cx=-1.20 XX=0 .
COMPUTE PRED_ = Cx * ( 1 - ( ( XX / ( XX + x ) ) ** 2 ) ) - c0.
NLR c
/OUTFILE='C:\SPSSFNLR.TMP'
/PRED PRED_
/CRITERIA SSCONVERGENCE 1E-8 PCON 1E-8 .
* pa * (1-EXP(-pb * x)) + pc * (1-EXP(-pd * t)).
* NonLinear? Regression.
MODEL PROGRAM c0=0 PA=0 PB=0 PC=0 PD=0 .
COMPUTE PRED_ = pa * (1-EXP(-pb * x)) + pc * (1-EXP(-pd * t)) - c0.
NLR c
/OUTFILE='C:\SPSSFNLR.TMP'
/PRED PRED_
/CRITERIA SSCONVERGENCE 1E-8 PCON 1E-8 .
* PCX * ( 1 - ( ( PX / ( PX + x ) ) **2 ) ) * ( 1 + PM * ( 1- (( PT / (PT + t ) ) **0.3 ) )).
* NonLinear? Regression.
MODEL PROGRAM c0=0 PCX=0 PX=0 PT=0 PM=0 .
COMPUTE PRED_ = PCX * ( 1 - ( ( PX / ( PX + x ) ) **2 ) ) * ( 1 + PM * ( 1- (( PT / (PT + t ) ) **0.3 ) )) - c0.
NLR c
/OUTFILE='C:\SPSSFNLR.TMP'
/PRED PRED_
/CRITERIA SSCONVERGENCE 1E-8 PCON 1E-8 .
오전 1:42 2004-06-04, by Kenial
뭔가 되는 것 같기도 하고 아닌 것 같기도 하고.. 환장하겠네.~
Cm = C(x,t) - C0 ... 아 젠장 수식도 제대로 체크를 안하고 뻘짓하다니 ...
Levenberg-Marquardt Method이 아니라면... sequential quadratic programming이다!
오후 4:18 2004-06-03, by Kenial
Levenberg-Marquardt Method에 대해서...
http://groups.google.co.kr/groups?hl=ko&lr=&ie=UTF-8&newwindow=1&threadm=8kjf16%24mnf%241%40b5nntp2.channeli.net&rnum=3&prev=/groups%3Fq%3Dlevenverg-marquardt%2520algorithm%26hl%3Dko%26lr%3D%26ie%3DUTF-8%26newwindow%3D1%26sa%3DN%26tab%3Dwg
http://groups.google.co.kr/groups?q=SNLS1&btnG=%EA%B5%AC%EA%B8%80+%EA%B2%80%EC%83%89&hl=ko&lr=&ie=UTF-8&newwindow=1
http://www-fp.mcs.anl.gov/otc/Guide/OptWeb/index.html
오전 12:17 2004-06-03, kenial
- 데이터 입력.
- 모델 1~5를 통해 모수의 값을 추정.
- 추정된 모수의 값으로 비선형 회귀분석 시도.
- 회귀분석의 Squared r값을 이용해 적절한 모델 결정.
- 다중공선성 체크.
- plot 출력.
analyze - regression - curve estimation 후 함수의 타입을 테스트...!
그다음 나온 b1 값을 토대로 nonlinear regression을 쌔운다
상수항 포함/비포함을 참고
(한글 spss 10.0에 의한 알기 쉬운 다변량분석, 노형진 저, 형설출판사)
오후 11:45 2004-05-27, kenial
샘플 데이터를 이용, c = PCX * ( 1- ((PX / (PX + x ))**2)) * ( 1 + PM * (1- (PT / (PT + t))**0.3)) 방정식의 회귀분석을 시도했으나, PT + t 값의 0일 때가 있어 divide by zero 에러를 내고 사망.
PT > a 의 일정값을 지정한 후 Squared R 값의 변화 :
| a |
Squared R |
| 1E-20 |
.04616 |
| 1E-19 |
.04616 |
| 1E-18 |
.56504 |
| 1E-17 |
.56458 |
| 1E-16 |
.56465 |
| 1E-15 |
.56478 |
| 1E-14 |
.56420 |
| 1E-13 |
.56447 |
| 1E-12 |
.56293 |
| 1E-11 |
.56529 |
| 1E-10 |
.56818 |
| 1E-09 |
.59836 |
| 1E-08 |
.62288 |
| 1E-07 |
.77039 |
| 1E-06 |
.63823 |
| 1E-05 |
.74458 |
| 1E-04 |
.89697 |
| 1E-03 |
.89542 |
| 1E-02 |
.89231 |
| 1E-01 |
.89485 |
| 1 |
.89521 |
| 1E+01 |
.85480 |
| 1E+02 |
.62513 |
도대체 이 결과를 어떻게 받아들여야 할지 알 수 없음.
오후 11:45 2004-05-27, kenial
2차 미팅
- 각자 통계기초 스터디에 돌입
- R-project를 기반의 인터페이스를 통해 라이브러리를 이용해서 프로그램을 작성 : R-project에서 제공하는 프로그래밍 인터페이스에 관련된 내용을 찾아볼 것
- 라이브러리를 통해서 나온 분석 결과를 해석할 수 있는 수준까지 통계관련이론 학습 ; Leverberg-Marquardt 방식의 다중회귀분석을 사용할 것
- 분석 결과를 이용, 그래프 출력이 가능한 모듈 작성, 혹은 R-project의 package를 이용할 것.
- 라이브러리를 그대로 이용하느냐, 필요한 내용을 짜집기할 것이냐.
- 페이 관련 협의 ( ㅇㅂㅇ)/
- Deadline : 2004-06-23(수)
- 계측 데이터 수령함 ; 페이퍼 형태라 수작업 입력 필요. 샘플링 계획 세울 것.
- 논문 자료 수령 ; 관련자료 - 터널/계측 관련자료에 추가.
오전 10:00 2004-05-27, kenial
-
![[WWW]](http://kenial.net/wiki/img/moin-www.gif)
r-project 개발 환경에 대해 조사.
; r 언어 및 (개발) 환경 - 통계 처리와 그래픽에 중점을 둠. s language의 gnu 버전으로 만들어졌으며(참고:s+는 s language의 개량판이며, 일반적인 s language의 코드는 큰 수정 없이 r-project 환경에서도 사용 가능) 기본적인 s 언어와 환경을 제공하는 것 외에 contributed packages 라는 형태로 외부 패키지를 추가함으로서 기능을 확장할 수 있다. 그 외에도 서브 프로젝트들을 갖고 있으며, gui 처리 환경을 제공하거나 graphical model을 쉽게 제작할 수 있는 목적의 프로젝트, 혹은 LaTex?를 이용해 출력가능한 형태의 문서를 만들어내는 패키지 등 여러 형태의 애드인이 제공됨.
- MathML(Mathematical Markup Language) 관련 내용 정리할 것.
오후 6:35 2004-05-26, kenial
1차 미팅
- 개발 형태 결정 - 단일 어플리케이션 스타일, c++ based.
- 필요한 사항 요구
- 입력 데이터의 포맷
- 관련 자료 (회귀분석, 터널 붕괴와 관련된 논문/연구자료)
- 개발 시일 - 외부 라이브러리 사용과 관련 있음
- 페이 관련
- 개발 모듈 분리
; 입력 데이터 파서, 수식 데이터 저장/Equatation 출력/Graph 출력 클래스, Regression 프로세싱 모듈
기타분류
EditText
of this page
(last modified 2004-06-07 01:04:36)
FindPage
by browsing, searching, or an index
Or try one of these actions:
DeletePage,
DeleteUploadedFile,
LikePages,
SpellCheck,
UploadFile
