hello,everyone.
I tried to use cluster analysis package by reading the manual "chapter 15
cluster package"
.But when I tried to use the function "agnes",the R said :there is no
function.Can
you tell me which package I should load? Can anyone give me some sample
scripts about
clustering analysis ?T
Hello ,everyone.
I am a fresh user of R. I wrote several several R script in R with the
version 2.9.1 in Windows XP. But I only can run these scrpit via R-GUI.
Can anyone tell me how to run the scripts in background without openning
R-gui.In another word, run the R-script and get the result
Hi ,everyone ,
I draw some boxplot figure with the command "boxplot".But in the
figure,there are some bubbles at the top part of the figure.
Can anyone tell me what the correct meaning of these bubbles?and how to
remove it?
--
TANG Jie
Email: totang...@gmail.com
Tel: 0086-2154896104
Shanghai
hi everyone,
I am wondering if I can inintial all of the varibles and vectors at the
begining part of my R script to make clear all the meanning defined
in another script runned before It is just similiar with the command " close
all " in Matlab or "reinit " in grads.How can I do?
thank you.:
hi ,everyone,
I have a script shown as below:
> bj2=bjerrdata$tyerr[bjyearnum[2]+1:bjyearnum[3]]
> length(bj2)
[1] 448
while
> b=bjyearnum[3]
> a=bjyearnum[2]+1
> bj1=bjerrdata$tyerr[a:b]
> length(bj1)
[1] 169
it is different with bj2 and bj1 . and the array bjyearnum is
[1]0 279 448
interting ,
when I modified
bj2=bjerrdata$tyerr[bjyearnum[2]:bjyearnum[3]]
the length of bj2 is correct .
It seems that can not put a fomula in the dimension description .
2009/7/29 Jie TANG
> hi ,everyone,
>
> I have a script shown as below:
>
> > bj2=bjerrdata$ty
Hi ,everyone.
I want to draw a boxplot figure, main script shown as below. when I run the
command "boxplot(bjavee~yearname)",it is wrong with the error information "
model.frame.default(formula = bjavee ~ yearname) :
variation 'bjavee' list is wrong".
I want to add the yearname as the x-axis ti
Hi ,everyone.
I want to draw a boxplot figure, main script shown as below. when I
run the command "boxplot(bjavee~yearname)",it is wrong with the error
information " model.frame.default(formula = bjavee ~ yearname) :
variation 'bjavee' list is wrong".
I want to add the yearname as the x-axis t
hi ,everyone
I want subtract some dataset by subset.
>From the help running help(subset), ths information is "*subset(airquality,
Day == 1, select = -Temp)* "
while I running my script written as "*g1data<-subset(errdata, fac>12) *"
,it is wrong with the error information "*subset.default(
hello ,every one ,
I draw a figure as shown in appendix A,
The tick characters in the Y-axis is arranged in the vertical
direction.Now I want to arranged the Y-axis tick in the horizontal
direction.
How could I do ,then?
Thank you .
--
TANG Jie
Email: totang...@gmail.com
Tel: 0086-2154896104
Hello,
Now I draw a figure by boxplot ,but I am fail to add the tick information in
the axis? How could I ?
https://mail.google.com/mail/?ui=2&ik=c88089ec0a&view=att&th=1266e2c8dff41fcb&attid=0.1&disp=attd
now I want to add the tick information in the x-axis as Class1
,Class2,Class3 instead of "
I try to install a R resource package (in R.2.15 or R 3.0)
But when I run the command "configure",the Makefile does not
generated succefully as error information shown as below
--with-readline=yes (default) and headers/libs are not available
my linux sysem is "Red Hat Enterprise Linux Client rel
hi R users
I have a data set with the name AA
AA<-1:100
Now I want to get a child array from AA every 10 numbers
e.g.
ab =c(10,20,30,40,50,60,70,80,100)
How could I subtract aa from AA?
thank you .
--
TANG Jie
Email: totang...@gmail.com
Tel: 0086-2154896104
Shanghai Typhoon Institute,China
Actually, I find a method. The function reshape" OR array
data_5m<-array(data,c(5,length(data)/5))
and
data_5m[1,] is what I wanted ,but I donot think it is a good method.
2013/1/14 Jeff Newmiller
> ab <- AA[(1:10)*10]
--
TANG Jie
Email: totang...@gmail.com
Tel: 0086-2154896104
Shanghai Ty
hi Rusers
I am trying to use symbos in gtools package
symbols(data1,data3,circle=data1/data3,inches=0.1,bg="lightgreen")
Now I want to give a lengend to tell the reader the meaning or magnitude of
these circle.
How can I add these information in symbols plot just like legend in plot ?
thank y
hello ,Rusers:
I have a dataset for example
totdata=c(8.020 7.412 8.110 8.270)
plot(totdata[1:4],type='l',pch=1,lty=1,col='blue',lwd=1)
and the data could be plotted as the attachment figure.
now I want to exchange its x-axis and y-axis in the figure .How could I
modify the plot command ?
than
hi R users:
I have a dataset including u wind in x-axis and v wind in y-axis.
How can I plot the u,v wind data in vector or barb figure?
which command ?
thank you .
--
TANG Jie
[[alternative HTML version deleted]]
__
R-help@r-project.org
hi R-users,
I have a file as shown below,and the first column is time information.
"2012-08-07 00:00:00",4174830,5,8.1,34.5,9.5,32,14
"2012-08-07 00:00:01",4174831,4.7,8.6,34.5,9.9,29,14
"2012-08-07 00:00:02",4174832,4.7,8.6,34.5,9.9,29,14
"2012-08-07 00:00:03",4174833,5,8.5,34.5,9.8,30,14
"2012
Hi R users:
I have a big data and want to calculate the running mean of this data .
How can I get this kind of result ? I have check the command "mean"
it seems "mean" could not get the running mean?
--
TANG Jie
[[alternative HTML version deleted]]
forecast', dependencies = TRUE)
> library(forecast)
>
> ?ma
>
>
> Hope this helps,
>
> Rui Barradas
>
>
> Em 12-06-2013 08:21, Jie Tang escreveu:
>
> Hi R users:
>>I have a big data and want to calculate the running mean of this data .
>>
hi r users:
I have a datadset and want to write into a file .
when i use :
write.table(u_bar, file = "u_test.txt")
the data in the outpuf file is shown as below which included the
row number from "1" to "14" and colname .e.g."X32N"
How could I output the numbers and excluded the colname and ro
hi R users
I have a large 1D dataset and some of them is NA value .
I found I cound get the spectrum by such a command.
ua=c£¨10£¬30 £¬40£¬50£¬NA£©
spectrum(ua)
and I could not use na.rm just like mean or sd function
How could I get the spectrum of ua ?
thank you .
--
TANG Jie
[[a
hi
I used a two nested cycle by if and for by such code
for(ity in 1:4)
{
if (ity==1£©
{
print(ity)
}
}
when I run the code it failed and R tell me that
"error: unrespected '}' in "}""
and when I reduce a }
for(ity in 1:4)
{
if (ity==2£©
{
print(ity)
}
R will print 4 but not 2 as what I repec
hi
I have a wind dataset and I want to analze its spectrum with
the x-axis and y-axis in log .But I found that
when I set log="dB" in spectrum as shown below,what I can get
is a figure by log y-axis.
How could I get a log-axis figure both in x-axis and y-axis ?
thank you .
yylab<-c("Wind Spectrum(
Futher more
When I write as below :,the result is still not satisfied .
asp=spectrum(ww_spec)
windspec=data.frame(asp)
plot(windspec,log="dB")
2014/1/7 Jie Tang
> hi
> I have a wind dataset and I want to analze its spectrum with
> the x-axis and y-axis in log .But I foun
I use
windlog=log10(windspec)
plot(windlog)
it seems good but stupid
is there any better method ?
2014/1/7 Jie Tang
> Futher more
>
> When I write as below :,the result is still not satisfied .
> asp=spectrum(ww_spec)
> windspec=data.frame(asp)
> plot(windspec,log=&q
Hello,everyone .
Now I write some R scripts to calculate some variables and I run them in C
shell script command in linux..
Now my some new R scripts had to be moved in someother's computer with
windows system.
How to run these R scripts easily without installing many soft package?
--
TANG Ji
here ,I want to plot two lines in one figure.But I have two problems
1) how to move one of the y-axis to be the right ? I tried to the
command"axis(2)",But I failed.
2) how to add the axis information correctly.Since I have use the cmommand
"axis(1,at=1:6,labels=gradeinfo$gradenam)"
but it seems
list ;-). Try searching the help archiwe (www.rseek.org) and you will
> find solutions. I would guess that you need to use something like:
>
> axis(4)
>
> as the sides of the plot are always numbered from bottom,left,top,right
>
> HTH
> Jannis
> Jie TANG schrieb:
>
>> he
hi ,R user folks .
Nowadays I read a paper which draw a probability ellipse circle figure
shown in the appendix.
I wonder how to draw this figure by R ?
the x-axis and y-axis both express the error but in different direction .
--
TANG Jie
Email: totang...@gmail.com
Tel: 0086-2154896104
Shanghai
hi ,R USERS:
I am a new R user. I have study some basic R syntax and get some
interesting results.
Now I want to know how much R can do .where can i get some examples or demo
resource to study advanced function of R ?
--
thank your
TANG Jie
Email: totang...@gmail.com
Shanghai Typhoon Institu
hi,R users,
Now I have some scatter figures, is there some method can plot equalprobable
error ellipse ?
--
TANG Jie
Email: totang...@gmail.com
Tel: 0086-2154896104
Shanghai Typhoon Institute,China
[[alternative HTML version deleted]]
__
R-
for example:
a data.frame
data
V1 V2 V3
1 2 3
4 5 6
If I want to calculate the V2
I can get data$V2
But now I want to calculate the dimension name and then analysis the data.
for example
No=3*x ;x=1,2,3,...n
anadata=data$VNo
How could I do?
--
TANG Jie
[[alterna
my qustion is if my miss value in my dataset is -
for example :
1 3 21 33 - 23 33 -
how can I plot or analysis this dataset ignoring the missing value ? thanks
.
--
TANG Jie
[[alternative HTML version deleted]]
__
R-help@r-pr
hi ,
Now I have two dataset and want to compare them with same sample.
Dataset A:
year mon day obs
2010 03 1212
2010 03 1822
2010 04 1262
2010 07 24 29
Dataset B:
year mon day obs
2010 03 1215
2010 04 1257
2010 07 2432
2010 08
) .OO#. .OO#. rocks...1k
> ---
>
> Sent from my phone. Please excuse my brevity.
>
> Jie TANG wrote:
>
>> hi ,
>> Now I have two dataset and want to compare them with same sample.
>> Dataset A:
>> year mon day obs
>> 2010 03 1212
&g
for example : I have files with the name
"ma01.dat","ma02.dat","ma03.dat","ma04.dat",I want to read the data in
these files into one data.frame
flnm<-paste("obs",101:114,"_err.dat",sep="")
newdata<-read.table(flnm,skip=2)
data<-(flnm,skip=2)
but the data only contains data from the flnm[1]
I als
hi, R-users
I have a data.frame for example test$newdataday24 and test$newdataday48
I can plot them by
plot(test$newdataday24)
but now i want to plot different data by define a variable to describe them
dayno<-c(24,48)
newnam<-paste("test$newdataday",dayno,sep="")
plot(newnam[1])
but i failed,
='')
>
> plot(test[[newnam[1]]])
>
>
> On Mon, Aug 29, 2011 at 12:29 PM, Jie TANG wrote:
>
>> hi, R-users
>> I have a data.frame for example test$newdataday24 and test$newdataday48
>> I can plot them by
>> plot(test$newdataday24)
>> but
hi R user
mtdno<-paste("data",1:3,sep="")
tyno<-paste("obs",1:5,sep="")
flnm<-paste(mtdno,tyno,"_err.dat",sep="")
flnm is
[1] "data1obs1_err.dat" "data2obs2_err.dat" "data3obs3_err.dat"
[4] "data1obs4_err.dat" "data2obs5_err.dat"
but actually what i want is data from 1 to 3 and obs from 1 to 5.
uot;)
flnm would be a 1-dimension with 15 elements?
thankyou
2011/9/1 Jorge I Velez
> Hi Jie,
>
> Try
>
> c(outer(mtdno, tyno, FUN = paste, "_err.dat", sep = ""))
>
> HTH,
> Jorge
>
>
> On Thu, Sep 1, 2011 at 3:11 AM, Jie TANG <> wrote:
>
nice. thank you all.
R is a great tool and R users are also so kind .
2011/9/1 David Winsemius
>
> On Sep 1, 2011, at 3:45 AM, Jie TANG wrote:
>
> thank you. it works.
>> but further question is that if we can let the " flnm" to be a
>> 2-dimension
>>
hi
i have a dataframe with the name "obsdata"
V1 V2 V3 V4V5V6V7V8 V9 V10 V11 V12 V13
11001 3 24 12 24.7 44.4 70.1 49.3 33.7 3.0 6.8 2.7NA
21001 3 25 0 70.1 49.3 33.7 138.2 152.5NA 4.2 6.9 17.5
31001 3 25 12 33.7 187.7
1) how to modify the the tickment of x-axis or y-axis.
boxplot(data[,1:5])
the tickment in x-axis in V1 V2 V3 V4 V5 ,I want to be some name for
example
name<-c("1day","2day","3day","4day","5day")
2) how to overlap two plot into one figure?
plot(data[1:5])
boxplot(newdata[,1:5])
?
--
TAN
hi ,when i read a lots of files
for (i in 1:totnum)
{
cop_x_data<-read.table(flnm[i],skip=2)
if(i==1) {cop_data=cop_x_data}
else {cop_data=rbind(cop_data,cop_x_data)}
}
some of the files are missing . so this loop can not go on .I wonder how can
I go on
the loop cycle while reading the files jus
hi ,R users
I read a series of file by the command shown as below.
cop_x_data<-read.table(flnm,skip=2)
the first two line of the files are headfile and I skip them.
and sometimes the original data file is null and there is no any data
in the file except for the head information.
At this situation
hi ,
I am reading a series of files by the command shown as below.
cop_x_data<-read.table(flnm(i)
,skip=2,allowEscapes=TRUE,blank.lines.skip=TRUE)
the first two line of the files are headfile and I skip them by skip=2.
and sometimes the data file is null and there is no any data
in the file except
hi,R-users:
Now I plot some data with the name(aveobsdata) in column , How can I add
the
some number(e,g. the sample number) in each of the column?
plot(aveobsdata,type='h',lwd=line_width,col=line_col,main=titleinfo,xlab=xxlab,ylab=yylab,xaxt
= "n")
axis(1, at = 1:nums, label = name)
--
_
hi R-users:
I have some data related with geographical information.
I wonder how to plot some data (in contour,shaded or scatter) on a map.
Could anyone tell me which package or reference or examples I can learn
?
thank you .
--
TANG Jie
Email: totang...@gmail.com
[[alternativ
hi R-users
Now I have a figure in emf or png or tiff format that have been drawn
by other tool and I want to insert this figure into
my new figure by R script. I wonder if is possible ?
--
TANG Jie
Email: totang...@gmail.com
Tel: 0086-2154896104
Shanghai Typhoon Institute,China
[[al
hi R-users:
I want to draw three plot into one figure by layout and the script has
been shown below.
But I find R does not show the xlab and ylab information completely as
shown the figure attached.
How can I midify the script.? thank you .
xxlab<-paste(cpmd," (",ro,"%)",sep=" ")
yylab<-pas
I defined a data.frame by a two-dimensional array.
> aa = data.frame(rbind( 11:20, 1:10))
> aa
X1 X2 X3 X4 X5 X6 X7 X8 X9 X10
1 11 12 13 14 15 16 17 18 19 20
2 1 2 3 45 6 78 9 10
Now I want to use the data in the second line of aa which from 1 to 10 to
output.But I fin
hi everyone .
Here I have a text where there are some integer and string variables.But I
can not read them by readLines and scan
the text is :
weight ;30;130
food:2;1;12
color:white;black
the first column is the names of the variables and others are the value of
them.
the column in different l
HI R-users
I downloaded the latest version of R with the name R-latest.tar.gz from
the website
But when I tar the package and configure the package,
some error message showed that
configure: error: --with-readline=yes (default) and headers/libs are not
available
and I run make command:
error
-- Forwarded message --
From: Jie TANG
Date: 2011/10/25
Subject: how can I install the latest version of r in linux?
To: r-help@r-project.org
HI R-users
I downloaded the latest version of R with the name R-latest.tar.gz from
the website
But when I tar the package and
hi everyone .
I have a group of files that needs to be read in .
And I put their filename into one file with the name "file.txt"
And in the file.txt,the filename are shown as below :
H20100902.TXT
H20100903.TXT
H20100904.TXT
H20100905.TXT
H20100906.TXT
H20100907.TXT
H20100908.TXT
Now I want to re
;
> >
> > files <- scan('file.txt', what = '') # read names from the file
> Read 7 items
> > files # display the names
> [1] "H20100902.TXT" "H20100903.TXT" "H20100904.TXT" "H20100905.TXT"
> "H20100906.TX
here is my file.txt that generated by other command "ls "
2012/3/7 Jie Tang
> there is something wrong :
>
> scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, :
> scan() needs 'a real', but not 'H20100903.TXT'
>
> 2012/3/7 j
0904.TXT"
> ,"H20100905.TXT", "H20100906.TXT",
> "H20100907.TXT","H20100908.TXT"),function(x) scan(x, other options))
>
> Kees
>
> On Wed, Mar 7, 2012 at 3:44 PM, Jie Tang wrote:
> > there is something wrong :
> >
> > scan(
hi
I plot a series of observation data every minutes in a day as the
attachment below
plot(wnd,type='l',lty=1,col='red',lwd=1,xlab=xxlab,ylab=yylab,ylim=YY)
In the figure, the x-axis tickment is the number of data How can I change
it
fore example 1h 2h 3h 4h and so on ?
--
TANG Jie
Email: to
hi everyone .
Now I want to draw several lines in one frame.And it seems needs more
colors except for blue red,black .Where can i found these color name or
define some new color ?thank you .
--
TANG Jie
Email: totang...@gmail.com
Tel: 0086-2154896104
Shanghai Typhoon Institute,China
[[
hi ,I am trying to using wavelet package.
require(wavelets)
wt <- dwt(aa, n.levels=4, boundary="reflection", fast=TRUE)
plot(wt)
and we can get the figure as the attachment.
aa is a data array that came from the one day obvervation in 1HZ and has
the length of 86400 .
Can anyone tell me what i
hi everyone .
I am trying to use some packages but there are some functions have the same
name in different package.
for example
dwt function
both in packages wavelets and waveslim
dwt(X, filter="la8", n.levels, boundary="periodic")
How can I avoid mixing them up ?
--
TANG Jie
Email: totang
hi
I plot a series of observation data every minutes in a day as the
attachment below
plot(wnd,type='l',lty=1,col='red',lwd=1,xlab=xxlab,ylab=yylab,ylim=YY)
In the figure, the x-axis tickment is the number of data How can I change
it
fore example 1h 2h 3h 4h and so on ?
someone has told me that
t;>
>> op <- par()
>> par(ask = TRUE)
>> plot(1:10) # Unhappy axis
>>
>> plot(1:10, xaxt = "n")
>> axis(1, at = 1:10, label = paste(1:10, "h", sep = "")) # Happy axis!
>> par(ask = op$ask)
>>
>> The (good) a
hi
I plot a series of observation data every minutes in a day as the
attachment below
plot(wnd,type='l',lty=1,col='red',lwd=1,xlab=xxlab,ylab=yylab,ylim=YY)
In the figure, the x-axis tickment is the number of data How can I change
it
fore example 1h 2h 3h 4h and so on ?
someone has told me that
;)
> axis(1, at = 1:10, label = paste(1:10, "h", sep = "")) # Happy axis!
> par(ask = op$ask)
>
> The (good) advice to use dput() was to help us help you. Reread what was
> said.
>
> Michael
>
> On Sat, Mar 17, 2012 at 10:22 AM, Jie Tang wrote:
>
67 matches
Mail list logo