[R] for loop output
stats12
skarmv at gmail.com
Wed Oct 17 20:16:53 CEST 2012
Dear R users,
In the code below, I am trying to print the result of my loop function. The
output first gives me the result for k=1, and then for k=1 and k=2. I only
want the last output which is
[,1] [,2]
[1,] 0.1700065 0.5002659
[2,] 0.3080273 0.4954731
[3,] 0.4844886 0.4544306
[4,] 0.5062987 0.1868154
[5,] 0.5846982 0.4353522
[6,] 0.4332621 0.2202922
[7,] 0.4391985 0.2424147
[8,] 0.4242643 0.2372497
[9,] 0.3274367 0.3664741
[10,] 0.1937940 0.4192495
How can I extract only this part of my output. Thank you very much in
advance.
d<-matrix(c(1,1,0,0,0,0,0,0,2,1,0,0,1,1,0,1,2,2,1,0),nrow=10,ncol=2)
h<-matrix(runif(20,0,1),10)
int<-matrix(c(0),nrow=10, ncol=2)
for (k in 1:2){
for(s in 1:10){
integrand<-function(x) x^d[s,k]*exp(-x*gamma(1+1/1.6))^1.6*exp(-x*h[s,k])
integ[s,k]<-quadgk(integrand,0,1000)
}
print(int)
}
[,1] [,2]
[1,] 0.1700065 0
[2,] 0.3080273 0
[3,] 0.4844886 0
[4,] 0.5062987 0
[5,] 0.5846982 0
[6,] 0.4332621 0
[7,] 0.4391985 0
[8,] 0.4242643 0
[9,] 0.3274367 0
[10,] 0.1937940 0
[,1] [,2]
[1,] 0.1700065 0.5002659
[2,] 0.3080273 0.4954731
[3,] 0.4844886 0.4544306
[4,] 0.5062987 0.1868154
[5,] 0.5846982 0.4353522
[6,] 0.4332621 0.2202922
[7,] 0.4391985 0.2424147
[8,] 0.4242643 0.2372497
[9,] 0.3274367 0.3664741
[10,] 0.1937940 0.4192495
--
View this message in context: http://r.789695.n4.nabble.com/for-loop-output-tp4646496.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list