[R] dev.off()

2007-11-21 Thread amna khan
Dear Sir

Sir when I use "png" function to save graph. At last it is written dev.off().
it does not produce any postscript file.

Please help in this regard



Thank you

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] grDevices package

2007-11-21 Thread amna khan
Dear Sir

grDevices package is not found in packages list?

Thank You

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] grDevices

2007-11-23 Thread amna khan
Dear Sir

I am using R 2.6.0. In R console grDevices are not giving results like
"png", "postcript", xlim, etc. I think they are not activated?

I request you to please provide guidance in this regard.

Regards

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] problem in closing R console

2007-11-25 Thread amna khan
Dear Sir

When I close R console then following message apear when I want to save the
workspace. If I close the R console without saving workspace then R console
is closed otherwise not.

Error in gzfile(file, "wb") : unable to open connection
In addition: Warning message:
In gzfile(file, "wb") : cannot open compressed file '.RDataTmp'
Regards

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] All files inactive in installing R

2007-11-26 Thread amna khan
Dear Sir

I was facing various problems especially in making graphs. No graph file was
displayed. I reinstall the R 2.6.0. During installation "all files" icon is
inactive. It has not been marked. Which is creating many problems
particularly grDevices are not responding.

How to solve this problem

Regards

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] gzfile

2008-01-21 Thread amna khan
Dear Sir
When I close the R console I face the following problem continuously

Error in gzfile(file, "wb") : unable to open connection
In addition: Warning message:
In gzfile(file, "wb") : cannot open compressed file '.RDataTmp'
Help is required in this regard.
-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] Kendall Test

2007-09-16 Thread amna khan
Dear Sir

In Kendall Package of R, when we use the function Kendall() or MannKendall()
then we get the value of Kendall's tau and a two-sided p-value. The question
to be asked is that is the two-sided p value for Z statistics defined by
z=(S-1)/sqrt(var(S)) where S is Kendall S statistics.

Regards


HAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] Plot for two Series

2007-10-01 Thread amna khan
Dear Sir

I want to plot the two time series having unequal number of observations,
say, one series having values from year 1960 to year 2000 and other having
values from year 1975 to year 2000. How to plot them on a single graph
sheet.

Regards


-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] plots

2007-10-01 Thread amna khan
DEAR SIR

How can I use ts.plt() function for time series of unequal length?

Regards

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


Re: [R] Plot for two Series

2007-10-02 Thread amna khan
Sir the data used in function ts.plot() is of the form


  Jan  Feb  Mar  Apr  May  Jun  Jul  Aug  Sep  Oct  Nov
1974 2134 1863 1877 1877 1492 1249 1280 1131 1209 1492 1621
1975 2103 2137 2153 1833 1403 1288 1186 1133 1053 1347 1545
1976 2020 2750 2283 1479 1189 1160 1113  970  999 1208 1467
1977 2240 1634 1722 1801 1246 1162 1087 1013  959 1179 1229
1978 2019 2284 1942 1423 1340 1187 1098 1004  970 1140 1110
1979 2263 1820 1846 1531 1215 1075 1056  975  940 1081 1294


But I have the data series of the form

   YEAR Khanpur
1   195252.1
2   195322.4
3195426.4
41955   132.8
5   195671.9
6   195710.4
7   1958 9.1
8   195949.8

because of which the sacle on x-axis are the numbers 1,2 3,..etc.

Please guid in this regard

Thank you




On 10/1/07, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
>
> See ?ts.plot
> Also ?plot.zoo and ?xyplot.zoo in the zoo package.
>
> On 10/1/07, amna khan <[EMAIL PROTECTED]> wrote:
> > Dear Sir
> >
> > I want to plot the two time series having unequal number of
> observations,
> > say, one series having values from year 1960 to year 2000 and other
> having
> > values from year 1975 to year 2000. How to plot them on a single graph
> > sheet.
> >
> > Regards
> >
> >
> > --
> > AMINA SHAHZADI
> > Department of Statistics
> > GC University Lahore, Pakistan.
> >
> >[[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.
> >
>



-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] Fwd: cluster analysis

2007-10-07 Thread amna khan
-- Forwarded message --
From: amna khan <[EMAIL PROTECTED]>
Date: Oct 6, 2007 11:02 AM
Subject: cluster analysis
To: Alberto Viglione <[EMAIL PROTECTED]>

Dear Sir

I am facing problem in using traceWminim() function for cluster analysis. A
sample of data frame that I am using is of the following form

   sitename latlong  elev map
1 Dir  35.12   71.51 4675   1472.10
2S.Sharif   34.4472.21 3155   1051.70
3Peshawar   34.01  71.35 1185  435.90
4   Kakul   34.11 73.15 4301  1330.17
5 Balakot   34.3372.21  3214 1652.20

Sir when I used the above function for 2 clusters as
> traceWminim(atsite,2)
then the following error message appears
Error in if (min.traceW.2 > traceW.1) { : missing value where TRUE/FALSE
needed

Please provide your guidance

Thank You
-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.


-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] Problem in Saving Workspace

2007-10-17 Thread amna khan
Hi Sir

When I close R console and save work space then the following message
appears

Error in gzfile(file, "wb") : unable to open connection
In addition: Warning message:
cannot open compressed file '.RDataTmp' in: gzfile(file, "wb")
How to deal with this?

Regards
-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] cluster analysis

2007-10-18 Thread amna khan
Hi Sir

How to perform cluster analysis using Ward's method and K- means clustering?

Regards

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] Stoping rule for hierarchical clustering

2007-10-29 Thread amna khan
Hi Sir

I am new user of 'cluster' and 'amap'  packages for hierarchical clustering.

I want to find the optimum number of clusters described by Mojena (1977)
'Hierarchical
grouping methods and stopping rules: An evaluation' The Computer
Journal, 20(4).

I did not find the such function for optimal choice of clusters.

Some suggestions or hints are required.

Thank you


-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] cluster analysis

2007-11-01 Thread amna khan
Hi Sir

How can we select the optimum number of clusters?

Best Regards

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] ctest package

2007-11-01 Thread amna khan
Hi Sir

I did find ctest package from CRAN.
Where it can be found?

Regards

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] RFA package

2007-11-01 Thread amna khan
Hi Sir

In RFA package, using reganalysis() function, How can we select few sites
from of our own choice? Mean to say that how can I select the first two and
then the last two sites from given 10 sites?

Regards

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] R code for L-moment digram

2007-11-19 Thread amna khan
Dear Sir I am getting errors in using following R code for L-moment ratio
diagram
Help in this regard

> "plotlmrdia" <-
+ function(lmr,
+  nopoints=FALSE,
+ nolines=FALSE,
+ nolimits=FALSE,
+ nogev=FALSE,
+ noglo=FALSE,
+ nogpa=FALSE,
+ nope3=FALSE,
+ nogno=FALSE,
+ noexp=FALSE,
+ nonor=FALSE,
+ nogum=FALSE,
+ nouni=FALSE,
+  xlab="L-SKEWNESS",
+  ylab="L-KURTOSIS") {
+plot(lmr$limits, xlab = xlab, ylab = ylab, type = "n",font.lab=2)
+if(nolimits == FALSE) {
+ lines(lmr$limits,lwd=2.5,col=8)
+}
+if(nolines == FALSE) {
+  if(nogev == FALSE) lines(lmr$gev, col=1,lty=2,lwd=2)
+  if(noglo == FALSE) lines(lmr$glo, col=6,lwd=2)
+  if(nogno == FALSE) lines(lmr$gno, col=4, lty=2,lwd=2)
+  if(nogpa == FALSE) lines(lmr$gpa, col=2,lty=2,lwd=2)
+  if(nope3 == FALSE) lines(lmr$pe3, col=3,lwd=2)
+}
+if(nopoints == FALSE) {
+  if(noexp == FALSE) points(lmr$exp,pch=16,col=1)
+  if(nonor == FALSE) points(lmr$nor,pch=15,col=1)
+  if(nogum == FALSE) points(lmr$gum,pch=17,col=1)
+  if(nouni == FALSE) points(lmr$uniform,pch=18,cex=1.5,col=1)}}
>
> png("plot.png",width=480,height=400)
> par(mar=c(5,4,4,8))
>
> plotlmrdia(lmr, nopoints=FALSE, nolines=FALSE, nolimits=FALSE,
+ nogev=FALSE, noglo=FALSE, nogpa=FALSE, nope3=FALSE,
+ nogno=FALSE, noexp=FALSE, nonor=FALSE, nogum=FALSE,
+ nouni=FALSE)
Error in plot(lmr$limits, xlab = xlab, ylab = ylab, type = "n", font.lab =
2) :
  object "lmr" not found
>  grid()
Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) :
  plot.new has not been called yet
> par(xpd=TRUE)
> legend(1.15,1.05,c("LmL*","GEV","GLO","GNO","GPA","PE3","EXP","GUM","UNI","N"),
lty=c(1,2,1,2,2,1,NA,NA,NA,NA),col=c(8,1,6,4,2,3))
Error in rect(left, top, r, b, angle = angle, density = density, ...) :
  plot.new has not been called yet
> points(rep(1.3,4),c(0.32,0.22,0.12,0.02),pch=c(16,17,18,15),col=1)
Error in plot.xy(xy.coords(x, y), type = type, ...) :
  plot.new has not been called yet
> text(1.6,-0.2,"*L-moment Limits")
Error in text.default(1.6, -0.2, "*L-moment Limits") :
  plot.new has not been called yet
> par(xpd=FALSE,mar=c(5,4,4,2))
> dev.off()
windows
  2
THANK YOU

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] Scatter Plot

2009-07-30 Thread amna khan
Dear Sir

I want to write the numbers 1,2,3,on a scatter plot instead of points,
like 1 corresponding to first point on plot, 2 corresponding second point
etc.
Help in this regard.

Regards

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.


[R] Wakeby Curve

2009-08-04 Thread amna khan
Hi Sir
How to get Wakeby distribution curve on L-moment ratio diagram?

Regards

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

[[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.