On Oct 21, 2010, at 12:52 PM, Etn wrote:

Hi All,

Apologies for the simplicity of my question, but I would be grateful for any
advice. Thanks

I'm trying to put the output from a for loop into a data frame, however I
have not been successful.

The steps I have taken are:



*R-code:*

for (k in 1:(nt-1-n0) ){
  n<- n0-1+k
  lam=n/nt
  Q=x[n]

so n=n0+1+(1:(nt-1-n0))
and lam=(n0+1+(1:(nt-1-n0))/nt
and Q=x[0+1+(1:(nt-1-n0))]

Using the parameters below:

nt=20
n0=5
x=c(1:20)

> output3 = data.frame(k = 1:(nt-1-n0),
+                      n = n0+1+(1:(nt-1-n0)),
+                      lam = (n0+1+(1:(nt-1-n0)))/nt ,
+                      Q = x[0+1+(1:(nt-1-n0))] )


(I don't get what you calculated but you should be able to work toward a solution.)
> output3
    k  n  lam  Q
1   1  7 0.35  2
2   2  8 0.40  3
3   3  9 0.45  4
4   4 10 0.50  5
5   5 11 0.55  6
6   6 12 0.60  7
7   7 13 0.65  8
8   8 14 0.70  9
9   9 15 0.75 10
10 10 16 0.80 11
11 11 17 0.85 12
12 12 18 0.90 13
13 13 19 0.95 14
14 14 20 1.00 15

output1<-data.frame(cbind(k,n,lam,Q))
output1
}

output1


*R-Output *
    k    n  lam   Q
1 14  18  0.9   18


I would like the output in this format, but for all the values of k (i.e.
1-14 as opposed to just the last value)


I have also tried

*R-code:*

nt=20
n0=5
x=c(1:20)
for (k in 1:(nt-1-n0) ){
    n<- n0-1+k
    lam=n/nt
    Q=x[n]

output1<-data.frame(cbind(k,n,lam,Q))
print(output1)
}


*R-Output  *

   k n  lam Q
1 1 5 0.25 5
  k n lam Q
1 2 6 0.3 6
  k n  lam Q
1 3 7 0.35 7
  k n lam Q
1 4 8 0.4 8
 k n  lam Q
1 5 9 0.45 9
 k  n lam  Q
1 6 10 0.5 10
 k  n  lam  Q
1 7 11 0.55 11

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

David Winsemius, MD
West Hartford, CT

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to