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
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
commandaxis(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
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:
here ,I want to plot two lines in one figure.But I have two problems
1) how to move one of the y-axis
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
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=2ik=c88089ec0aview=attth=1266e2c8dff41fcbattid=0.1disp=attd
now I want to add the tick information in the x-axis as Class1
,Class2,Class3 instead of
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.:)
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 totang...@gmail.com
hi ,everyone,
I have a script shown as below:
bj2=bjerrdata$tyerr[bjyearnum
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 tick.
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
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, fac12) *
,it is wrong with the error information
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
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
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]]
__
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
---
Sent from my phone. Please excuse my brevity.
Jie TANG totang...@gmail.com wrote:
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
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 also tried as below
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,the
]]])
On Mon, Aug 29, 2011 at 12:29 PM, Jie TANG totang...@gmail.com 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 now i want to plot different data by define a variable to describe
them
dayno-c(24,48
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. thus ,I
can read 15
...@gmail.com
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:
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
nice. thank you all.
R is a great tool and R users are also so kind .
2011/9/1 David Winsemius dwinsem...@comcast.net
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
matrix [3,5]?
since
mtdno
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])
?
--
TANG Jie
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
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
.TXT H20100903.TXT H20100904.TXT H20100905.TXT
H20100906.TXT
[6] H20100907.TXT H20100908.TXT
On Wed, Mar 7, 2012 at 8:49 AM, Jie Tang totang...@gmail.com wrote:
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
here is my file.txt that generated by other command ls
2012/3/7 Jie Tang totang...@gmail.com
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 jim holtman jholt...@gmail.com
,H20100908.TXT),function(x) scan(x, other options))
Kees
On Wed, Mar 7, 2012 at 3:44 PM, Jie Tang totang...@gmail.com wrote:
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 jim
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:
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 is
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:
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
()
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) 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 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
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 totang...@gmail.com wrote:
hi
I plot a series of observation data every minutes in a day as the
attachment below
plot(wnd,type
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 totang...@gmail.com
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
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
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 jdnew...@dcn.davis.ca.us
ab - AA[(1:10)*10]
--
TANG Jie
Email: totang...@gmail.com
Tel:
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
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 ?
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
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
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
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= )
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-08-07
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]]
', 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 .
How can I get this kind of result ? I have check the command mean
it seems mean
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 row
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
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 find
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
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]]
__
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]]
__
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
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 repect
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(db)
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 totang...@gmail.com
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
I use
windlog=log10(windspec)
plot(windlog)
it seems good but stupid
is there any better method ?
2014/1/7 Jie Tang totang...@gmail.com
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
65 matches
Mail list logo