I appreciate your kind guidance! I did not read the manual carefully
(it's my fault).
Thank you so much, Prof. John Fox!
Chel Hee Lee
On 01/21/2016 12:52 AM, Fox, John wrote:
Dear Chel Hee Lee,
With the formula method, the default na.action is na.omit; thus,
aggregate(y~grp, data=tmp
What a surprise! Is this a bug? I would appreciate if you share the
results after testing the codes. Thank you so much for your helps in
advance!
Chel Hee Lee
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.
of translating R manuals. Comments and corrections is of
course most welcome.
Chel Hee Lee
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R
Thank you so much, David! I also confirmed that the link moved is not
found.
Chel Hee Lee
On 15-06-24 04:15 PM, David Winsemius wrote:
On Jun 24, 2015, at 10:52 AM, Chel Hee Lee wrote:
Could you also kindly check the following links in the rw-FAQ manual at
http://cran.r-project.org/bin
Thank you so much, Jim. Yes, I also see that the wiki
(rwiki.sciviews.org) is not found.
Chel Hee Lee
On 15-06-25 04:54 AM, Jim Lemon wrote:
Hi Chel Hee,
The last link (for FAQ 7.39) is dead and the domain
(rwiki.sciviews.org) is not there either.
Jim
On Thu, Jun 25, 2015 at 3:45 AM, Chel
://windowsvistablog.com/blogs/windowsvista/archive/2007/01/23/security-features-vs-convenience.aspx
I appreciate your help!
Chel Hee Lee
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
,
* Bioconductor software packages
http://www.bioconductor.org/packages/bioc/
Under the section 7.39 How do I create a plot with two y-axes?,
* http://rwiki.sciviews.org/doku.php?id=tips:graphics-base:2yaxes
I appreciate your helps!
Chel Hee Lee
__
R
A~ha~!! Thank you, Prof. Peter Dalgaard, so much for your wonderful
lesson!!! Learning new things everyday from this R-help mailing list!
Chel Hee Lee
On 2/11/2015 10:37 AM, peter dalgaard wrote:
On 11 Feb 2015, at 17:11 , Chel Hee Lee chl...@mail.usask.ca wrote:
The functional form
.
Chel Hee Lee
On 2/10/2015 2:29 AM, Rolf Turner wrote:
On 10/02/15 14:04, Ssuhanchen wrote:
Hi!
I want to use R to calculate the variable x which is in a complex
equation
in below:
2
Σ[exp(-x/2)*(x^k)/(2^k*k!)]=0.05
k=0
how to solve this equation to get the exact x in R
lapply(1:2, function(x) t(A[rev(1:3),,x]))
[[1]]
[,1] [,2] [,3]
[1,] g d a
[2,] h e b
[3,] i f c
[[2]]
[,1] [,2] [,3]
[1,] p m j
[2,] q n k
[3,] r o l
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 02/05/2015 03:17 PM, Karim Mezhoud wrote:
Dear
hope this helps.
Chel Hee Lee
On 1/31/2015 3:06 AM, Camilla Timo wrote:
Hi there,
I have a problem, i have a dataset, in which there are time series from 2010-
to 2014, like this:
-2014-01-30 15:39:46
-2012-04-20 14:49:02
And so on .
I want to have a situation in which there are days of week
r_arch
;x64
Chel Hee Lee
On 1/30/2015 9:55 PM, Jim Lemon wrote:
Hi Ahmed,
Hmmm, this seems to work for me (R-3.1.2, Linux)
legend(0,2100, legend=c(2009,2010,2012,2013,2014),
col = 1,cex=1,lty=NA,pch=c(1,2,6,7,8),lwd=2,bty=n,pt.cex=2)
Jim
On Sat, Jan 31, 2015 at 8:49 AM
I hope this helps.
Chel Hee Lee
On 1/30/2015 11:52 AM, JS Huang wrote:
Hi,
Here is my implementation. Hope this helps.
b
[1] 1 2 3 4 5
c
[1] 1 2 1 3 5 4
sapply(b,function(x)ifelse(x==c,1,0))
[,1] [,2] [,3] [,4] [,5]
[1,]10000
[2,]01000
[3
Why don't try this?
nams-c(P2O5 ['%'], D50 [mu*m])
pdf(test.pdf)
for(i in 1:2)
+ plot(1,1,xlab=parse(text=(bquote(.(nams[i])
dev.off()
Is this what you are looking for? I hope this helps.
On 1/30/2015 7:54 AM, PIKAL Petr wrote:
Thanks
So the key is to put expression(s) directly into
paste(rep(v1, each=3), 1:3, sep=_)
[1] a_1 a_2 a_3 b_1 b_2 b_3 c_1 c_2 c_3
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 1/30/2015 7:34 AM, Knut Hansen wrote:
I have a vector of several character elements:
v1 - c(a, b, c)
I want each of these elements to be repeated
(parse(text=DFName))[,1])
[1] 5
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 1/29/2015 12:34 AM, Jeff Newmiller wrote:
This approach is fraught with dangers.
I recommend that you put all of those data frames into a list and have your
function accept the list
Or, you can also do the same job using 'colSums()' as shown in the below:
colSums(status==I)
2010 2011 2012
344
I hope this helps.
Chel Hee Lee
On 1/28/2015 7:31 PM, JSHuang wrote:
Hi,
I think you need quotation around I like the following:
status
2010 2011 2012
1
Hi Bert! yes, you are VERY correct!!! Why am I making this simple thing
so complicated??? ;) Thank you so much for your nice lesson!
Chel Hee Lee
On 01/28/2015 09:59 PM, Bert Gunter wrote:
eek!
Chel Hee,anything that complicated should engender fear and trembling.
Much simpler and more
.
Chel Hee Lee
On 01/28/2015 07:38 PM, William Dunlap wrote:
with(dat1, ave(integer(length(Date)), Date, FUN=seq_along))
[1] 1 1 2 1 1 2 1 2 1 2 1
Bill Dunlap
TIBCO Software
wdunlap tibco.com
On Wed, Jan 28, 2015 at 4:54 PM, Morway, Eric emor...@usgs.gov wrote:
The two datasets below
I am using just the first row of your data (i.e. ID1).
ID1 - c(A, A, T, G, C, T, G, C, G, T, C, G,
T, A)
do.call(c,lapply(tapply(ID1, gl(7,2), c), paste, collapse=))
1234567
AA TG CT GC GT CG TA
Is this what you are looking for? I hope this helps.
Chel Hee Lee
do.call(rbind, lapply(split(data, data$Name), function(x)
x[order(x$CheckInDate),][nrow(x),]))
Name CheckInDate Temp
John John 2014-04-01 99.0
Mary Mary 2014-03-01 98.1
Sam Sam 2014-04-01 97.5
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 01/23/2015 05:43
=a, INDEX=grp.id, FUN=sum)
1 2 3 4
3 3 15 15
aggregate(x=a, by=list(grp.id), FUN=sum)
Group.1 x
1 1 3
2 2 3
3 3 15
4 4 15
xtabs(formula=a~grp.id)
grp.id
1 2 3 4
3 3 15 15
I hope this helps.
Chel Hee Lee
On 1/23/2015 3:32 AM, Bert Gunter wrote:
If Jim's
x - c(hola mundo mundo);
table(unlist(strsplit(x, )))
hola mundo
1 2
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 1/22/2015 8:25 AM, bgnumis bgnum wrote:
Hi all,
I want to cout the different words in a text.
You see if the text is: hola mundo mundo
the error as below if 'is.null()' is used for
the fourth component.
is.null(A[[4]])
Error in A[[4]] : subscript out of bounds
# Why? It is because the list 'A' has only three components.
length(A)
[1] 3
I hope this helps. I would appreciate if you DO READ the posting guide.
Chel Hee Lee
On 1
are looking for? I hope this helps.
Chel Hee Lee
On 1/20/2015 9:23 AM, Shuhua Zhan wrote:
Hi Chel,
Thank you very much for your help!
I'm sorry I did not post my wanted output correctly. I only want the colnames
of the data frame occur once as a table:
trait locus chr pos ci.low ci.high lod
bp
), col=red)
panel.text(1, -1, bquote(R^2 == .(summary(lm(y~x))$r.squared)))
},
grid = TRUE
)
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 01/17/2015 07:08 AM, Naresh Gurbuxani wrote:
In a conditional xyplot, I would like to add some numerical results in every
panel
.loc45 7 47.7 36.7156.7 6.11
2bp c15.loc8 15 12.0 3.9622.8 5.29
[[2]]
trait locus chr pos ci.low ci.high lod
1hr c2.loc54 2 59.8 14.887.8 4.19
2hr D15MIT184 15 22.8 12.036.0 3.15
Is this what you are looking for? I hope this helps.
Chel Hee Lee
222324
Hm this function 'extr.1()' seems to be (much) slower than the
function 'extr()'.
Chel Hee Lee
On 01/16/2015 11:34 AM, Bert Gunter wrote:
Chee Hee's approach is both simpler and almost surely more efficient,
but I wanted to show another that walks the tree (i.e
This approach may not be fancy as what you are looking for.
xl - unlist(x)
xl[grep(A, names(xl))]
f1.x1.A f1.x2.A f2.x3.A f2.x4.A
11 12 13 14
I hope this helps.
Chel Hee Lee
On 01/16/2015 04:40 AM, Rainer M Krug wrote:
Hi
Consider the following variable:
--8
for? I hope this helps.
Chel Hee Lee
On 01/16/2015 02:16 AM, philippe massicotte wrote:
Hi all.
How we evaluate a formula in R?
Ex.:
params - list(a = 2, b = 3)
x - seq(1,10, length.out = 100)
func1 - as.formula(y ~ a*x^2 + b*x)
##How to evaluate func1 using x and the params list
???
Thank
Thank you so much. Your explanation is very clear and helpful to
understand. Thank you, again.
Chel Hee Lee
On 1/12/2015 5:56 AM, Prof Brian Ripley wrote:
You don't want to change the locale, rather the message language. That
is not a locale, and e.g.
LANGUAGE=en
is the correct form
en_CA.utf8
Sys.getlocale(LC_MESSAGES)
[1] en_CA.utf8
a
エラー: オブジェクト 'a' がありません
Is this a bug? Any advice would be greatly appreciated.
Chel Hee Lee
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman
are looking for? I hope this helps.
Chel Hee Lee
On 12/22/2014 7:45 AM, najuzz wrote:
#Hi guys,
#I would like to count the number of individuals that receive X=0 troughout
their observational period.
#example dataset:
ID-c(1,1,1,1,2,2,3,3,3)
X-c(0,1,2,1,0,0,0,0,0)
Time-c(1,2,3,4,1,2,1,2,3)
Test
you are looking for? I hope this
helps.
Chel Hee Lee
On 12/21/2014 01:25 PM, eliza botto wrote:
Dear UseRs,
A point was plotted by the following command
plot(2,4,ylim=c(0,10),xlim=c(0,5))
how to divide the space around the plotted point into six regions each of 60
degree as shown
You are welcome, I am glad that I was able to help.
Chel Hee Lee, PhD
Biostatistician and Manager
Clinical Research Support Unit
College of Medicine
University of Saskatchewan
On 12/20/2014 03:53 AM, Kristi Glover wrote:
Thank you Prof. Lee for your code. I am sorry that I noticed that I sent
0 B 0
2015-04-07
x1 0
x2 0
x3 0
x4 0
x5 0
x6 A
x7 0
x8 0
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 12/19/2014 10:45 PM, Kristi Glover wrote:
Hi R User, Would you suggest me on how I can build a pivot table using two variables? I want to put
] 1,20001 5001,20001 10001,20001 15001,20001 20001,20001
[16] 1,25001 5001,25001 10001,25001 15001,25001 20001,25001
[21] 25001,25001
I hope this helps. Thank you, David, for your good lesson.
Chel Hee Lee, PhD.
Biostatistian and Manager
Clinical Research Support Unit
College of Medicine
RK.Records_CL V.Other
1005317NA NA
1007183NA NA
1008833NA NA
1012281NA NA
1015285 1 NA
1015315NA NA
1015322NA NA
Chel Hee Lee
On 12/18/2014 10:02 AM, Jeff Newmiller wrote:
No guarantees
it if you would provide
reproducible examples next time.
Chel Hee Lee
On 12/18/2014 11:48 AM, Jeff Newmiller wrote:
Read the posting guide. The solution is likely to depend on your operating
system and graphics devices
. I hope this helps.
Chel Hee Lee
On 12/18/2014 08:57 PM, Crombie, Burnette N wrote:
Chel, your solution is fantastic on the dataset I submitted in my question but
it is not working when I import my real dataset into R. Do I need to vectorize
the columns in my real dataset after importing? I
$ID
tmp
$MCZ4325
[1] C23.2 T43.2
$GDR2343
[1] M20.64 Y32.1 M20.64 T44.2
$BZD2643
[1] B83.2 T43.2 B83.2 Y32.1 B83.2 T44.2
$BCM3455
[1] B83.2 T43.2
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 12/18/2014 07:41 AM, Michael Dewey wrote:
Not sure how much help
1.988448
CID_10.32 CID_10 Actinobacteria1.644548 6 1.644548
CID_10.33 CID_10 Others1.582823 7 5.546493
I hope this helps.
Chel Hee Lee
On 12/07/2014 08:21 AM, Morway, Eric wrote:
Using the dataset dat (found below), I'm seeking a way to condense down
the data.frame
Anyway, move to next!
(continue)
I hope this helps, and thank you so much, Martin Morgan, for your good
lesson.
Chel Hee Lee
On 12/06/2014 07:22 PM, Martin Morgan wrote:
On 12/06/2014 02:53 PM, ce wrote:
Dear all,
Let's say I have this script , below. tryCatch indeed catches the
error
It seems that you would like to make a spaghetti plot (in a longitudinal
data analysis). You can use the function 'interaction.plot()'.
with(my.df, interaction.plot(TIME, ID, X))
I hope this helps.
Chel Hee Lee
On 12/06/2014 02:24 AM, arun wrote:
Not sure whether it is a scatterplot
for? I hope this helps.
Chel Hee Lee
On 12/05/2014 02:30 AM, Adrien Bonvin wrote:
Bonjour
Hi everybody,
Firs of all, sorry for my terrible English,
I would like to know if it’s possible to create an “empty plot” in which i
could add two ablines I created on two different plots earlyer in my
))
[1] 00200 00700 01200
I like to read the documentation by typing
help(formatC)
help(sprintf)
You may find answers what you wish to get. Documentation has been my
best friend when using R. I hope this helps.
Chel Hee Lee
On 12/05/2014 09:35 AM, Zilefac Elvis wrote:
Hi Chel,
How can I
list.files(pattern=*.dat)
[1] new_example_01.dat new_example_02.dat new_example_03.dat
[4] new_example_04.dat new_example_05.dat new_example_06.dat
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 12/04/2014 09:44 PM, Zilefac Elvis via R-help wrote:
Hello,
I would like to rename
I see that a function 'format()' is used in your code.
format(c(1,5,32,100), width=3, flag=0)
[1] 1 5 32 100
formatC(c(1,5,32,100), width=3, flag=0)
[1] 001 005 032 100
I hope this helps.
Chel Hee Lee
On 12/04/2014 10:54 PM, Zilefac Elvis wrote:
Hi Chel,
Thanks for the timely reply
haldate oiddate
$row.names
[1] 1 2 3 4 5 6 7
$class
[1] data.frame
I hope this helps.
Chel Hee Lee
On 12/03/2014 04:10 PM, Muhuri, Pradip (SAMHSA/CBHSQ) wrote:
Hello,
Two alternative approaches - mutate() vs. sapply() - were used to get the
desired results (i.e., creating a new column
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 12/03/2014 03:14 PM, Brock Huntsman wrote:
I apologize if this is a relatively easy problem, but have been stuck on
this issue for a few days. I am attempting to combine values from 2
separate dataframes. Each dataframe
to
get final bar plots. I hope this helps.
Chel Hee Lee
On 12/03/2014 12:29 PM, Charles R Parker wrote:
I am trying to create groups of barplots from data that have different number
of records in the groups, in such a way that all of the plots will have the
same numbers and sizes of bars
This question seems to be the problem specific to Ubuntu. What if you
post the message to r-sig-deb...@r-project.org?? I hope you get
answers from that mailing list.
Chel Hee Lee
On 12/02/2014 11:10 AM, VG wrote:
Hi everyone,
I was having trouble with R i installed some time ago on my
,between)
table(achtergrond$spits)
between morning evening
1636 579 142
You can also use function 'findInterval()' instead of using 'cut()'. I
hope this helps.
Chel Hee Lee
On 12/02/2014 02:04 PM, William Dunlap wrote:
You can do this in 2 steps - have cut() make a factor
xx - as.factor(c(AL, AK, CA, FL))
xx
[1] AL AK CA FL
Levels: AK AL CA FL
as.character(xx)
[1] AL AK CA FL
I hope this helps.
Chel Hee Lee
On 14-11-22 01:12 AM, Aditya Singh wrote:
Dear Boris and R-Experts,
I have a variable my_state which is a 2-letter character string telling which
Another (simpler) way that I can think is that
y * matrix(rep(z,3), ncol=ncol(y), byrow=TRUE)
[,1] [,2]
[1,]00
[2,]6 -3
[3,] 12 -6
I hope this helps.
Chel Hee Lee
On 14-11-19 08:43 AM, Ruima E. wrote:
Thank you Chel Hee.
Isn't there a simpler way to do so?
On Wed
y = matrix(cbind(c(0, 0.5, 1),c(0, 0.5, 1)),ncol=2)
z = matrix(c(12, -6),ncol=2)
t(apply(y, 1, function(x) x*z))
[,1] [,2]
[1,]00
[2,]6 -3
[3,] 12 -6
I hope this helps.
Chel Hee Lee
On 14-11-19 08:22 AM, Ruima E. wrote:
Hi,
I have this:
y = matrix(cbind(c(0, 0.5
/BatteriesO.O#. #.O#. with
/Software/Embedded Controllers) .OO#. .OO#. rocks...1k
---
Sent from my phone. Please excuse my brevity.
CHEL HEE LEE gnust...@gmail.com wrote:
I
I have some trouble to deal the value of 'NaN'. For example,
exp(1e3)
[1] Inf
exp(1e3)*0
[1] NaN
The correct answer should be 0 rather than NaN. I will very appreciate
if anyone can share some technique to get a correct answer.
__
Thank you, David. I found some good stuffs for which I have been
looking.
On Sun, 2012-08-26 at 09:01 -0700, David Winsemius wrote:
On Aug 25, 2012, at 9:00 PM, Chel Hee Lee wrote:
I will very appreciate if anyone can provide some materials to draw a
simplex plot of a Dirichlet
I will very appreciate if anyone can provide some materials to draw a
simplex plot of a Dirichlet distribution in R as shown in the page at
http://en.wikipedia.org/wiki/File:Dirichlet_distributions.png .
__
R-help@r-project.org mailing list
Oh... I see!!! I very appreciate for your clear explanation! Thank you, Uwe
Ligges.
--
View this message in context:
http://r.789695.n4.nabble.com/General-question-on-the-message-with-non-zero-exit-status-tp4551438p4568839.html
Sent from the R help mailing list archive at Nabble.com.
I have a general question how to understand the message with non-zero exit
status when new package is installed. Based on my experience, this message
implies a package dependence. Am I correct to understand this message? Are
there anyone who can provide some reference or explain details about the
62 matches
Mail list logo