R help - Mar 2004 - (no subject)

How to use "contour" function?
I type "contour(th1,th2,SumofSquares,levels=c(seq(1000,4000,200)))"

But I didn't see contour in the plot, I guess it is because my data is so 
big that contour is out of range.
Does anybody know how to fix this problem?


here th1 is a 1x101 array(158741.8-298529.6)
th2 is another 1X101 array(-0.0058352269-0.2093140935)

SumofSquares is a 101X101 matrix
here is 1st row of this matrix:
SumofSquares[1,]
   [1] 58276823733 45108959063 36097139203 29640961600 24852628699 
21205691351 18370693380 16132073683 14343415033
  [10] 12902043760 11733956375 10784525542 
10012575486  9386493133  8881608427  8478388277  8161165261  7917225781
  [19]  7736144795  7609292848  7529465556  7490601474  7487564677 
7515975345  7572076406  7652627560  7754820315
  [28]  7876209323  8014656457  8168284935  8335441457  8514664731 
8704659189  8904272893  9112478876  9328359317
  [37]  9551092042  9779938969 10014236185 10253385367 10496846372 
10744130780 10994796282 11248441758 11504702968
  [46] 11763248764 12023777744 12286015297 12549710985 12814636220 
13080582187 13347358003 13614789066 13882715573
  [55] 14150991196 14419481894 14688064834 14956627427 15225066452 
15493287270 15761203103 16028734388 16295808187
  [64] 16562357650 16828321525 17093643719 17358272883 17622162049 
17885268285 18147552386 18408978592 18669514322
  [73] 18929129935 19187798513 19445495656 19702199292 19957889512 
20212548406 20466159917 20718709708 20970185038
  [82] 21220574642 21469868628 21718058379 21965136456 22211096520 
22455933247 22699642261 22942220059 23183663956
  [91] 23423972019 23663143017 23901176371 24138072104 24373830797 
24608453551 24841941947 25074298011 25305524181
[100] 25535623276 25764598466



thanks,
Aimin Yan

Aimin Yan wrote:
> How to use "contour" function?
> I type
"contour(th1,th2,SumofSquares,levels=c(seq(1000,4000,200)))"
> 
> But I didn't see contour in the plot, I guess it is because my data is 
> so big that contour is out of range.
  It looks like it, given the data you showed. If you leave out the 
'levels' argument then the contour function tries to pick good values 
for the contour lines. Why dont you try that?

  contour(th1,th2,SumOfSquares)

Barry

Aimin Yan wrote:
> How to use "contour" function?
> I type
"contour(th1,th2,SumofSquares,levels=c(seq(1000,4000,200)))"
> 
> But I didn't see contour in the plot, I guess it is because my data is 
> so big that contour is out of range.
> Does anybody know how to fix this problem?

By either not specifying "levels", or by specifying it reasonable,
i.e.
values within the range of your matrix "SumofSquares".
Are there values < 4000 in that matrix? You should see a line in that case.

Uwe Ligges


> 
> here th1 is a 1x101 array(158741.8-298529.6)
> th2 is another 1X101 array(-0.0058352269-0.2093140935)
> 
> SumofSquares is a 101X101 matrix
> here is 1st row of this matrix:
> SumofSquares[1,]
>   [1] 58276823733 45108959063 36097139203 29640961600 24852628699 
> 21205691351 18370693380 16132073683 14343415033
>  [10] 12902043760 11733956375 10784525542 10012575486  9386493133  
> 8881608427  8478388277  8161165261  7917225781
>  [19]  7736144795  7609292848  7529465556  7490601474  7487564677 
> 7515975345  7572076406  7652627560  7754820315
>  [28]  7876209323  8014656457  8168284935  8335441457  8514664731 
> 8704659189  8904272893  9112478876  9328359317
>  [37]  9551092042  9779938969 10014236185 10253385367 10496846372 
> 10744130780 10994796282 11248441758 11504702968
>  [46] 11763248764 12023777744 12286015297 12549710985 12814636220 
> 13080582187 13347358003 13614789066 13882715573
>  [55] 14150991196 14419481894 14688064834 14956627427 15225066452 
> 15493287270 15761203103 16028734388 16295808187
>  [64] 16562357650 16828321525 17093643719 17358272883 17622162049 
> 17885268285 18147552386 18408978592 18669514322
>  [73] 18929129935 19187798513 19445495656 19702199292 19957889512 
> 20212548406 20466159917 20718709708 20970185038
>  [82] 21220574642 21469868628 21718058379 21965136456 22211096520 
> 22455933247 22699642261 22942220059 23183663956
>  [91] 23423972019 23663143017 23901176371 24138072104 24373830797 
> 24608453551 24841941947 25074298011 25305524181
> [100] 25535623276 25764598466
> 
> 
> 
> thanks,
> Aimin Yan
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! 
> http://www.R-project.org/posting-guide.html