Based on very fast converging series for special transcendental
numbers. Is there some R code available?
See details about the RNG at
http://www.analyticbridge.com/profiles/blogs/new-state-of-the-art-random-number-generator-simple-strong-and-fa
__
R-hel
I am trying to put an
R2 value with R2 formatted with a superscript 2 followed by "=" and the
value :
the first mtext prints the R2 correctly formatted but follows it with
"=round(summary(mylm)$r.squared,3)))" as text
the second prints "R^2 =" followed by the value of
round(summary(mylm)$r.squ
Hi all,
Not sure if this is the best list, please point me to a more appropriate list
if necessary.
Running Mac OSX 10.7.1
R version 2.13.1 Patched (2011-08-14 r56741)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-b
Hi
I have a data comprised of repeated measures of growth (5-15 records per
individual) for 580 fish (similar to Orange dataset from nlme library). I would
like to standardise these growth measures (yi â Å·/sd) using mean and
standard deviation unique to each fish. Can someone suggest a functi
Vince,
When that happens, one possible reason is that your DAE is of index > 1,
which cannot be solved by daspk.
The solver radau, also from deSolve can handle DAEs up to index 3, but
you need to rewrite the problem in the form M*y' = f(x,y), where M is a
mass matrix.
If you do that for your pr
Dear John,
You need to combine scale with a grouping function.
data(Orange)
library(plyr)
Orange <- ddply(Orange, .(Tree), function(x){
x$ddplyAge <- scale(x$age)[, 1]
x
})
Orange$aveAge <- ave(Orange$age, by = Orange$Tree, FUN = scale)
all.equal(Orange$ddplyAge, Orange$aveAge)
On Fri, Sep 30, 2011 at 2:41 PM, Duncan Murdoch
wrote:
> As an aside, I don't recommend the workflow you describe: it's very slow
> and cumbersome. It's much better to tell your text editor how to run both
> Sweave and Latex in one command. In the upcoming release of R 2.14.0, this
>
Another ap
Dear all
How does print.htest display the p-value in scientific notation?
> (x <- cor.test(iris[[1]], iris[[3]]))
Pearson's product-moment correlation
data: iris[[1]] and iris[[3]]
t = 21.65, df = 148, p-value < 2.2e-16
alternative hypothesis: true correlation is not equal to 0
95 percen
Hi Sandeep,
You need to provide way more detail for us to be able to help you:
- Which version of gstat do you want to install, ?
- where did you get the sources, ?
- is it the version inside R or the standalone version?
- what version of R are you running?
- *what is the exact error that gstat g
Isn't it true that 0 < 2.2e-16?
cheers,
Rolf Turner
On 03/10/11 20:53, Liviu Andronic wrote:
Dear all
How does print.htest display the p-value in scientific notation?
(x<- cor.test(iris[[1]], iris[[3]]))
Pearson's product-moment correlation
data: iris[[1]] and iris[[3]
> is(x)
[1] "htest"
> # take a look at stats:::print.htest
> format.pval(x$p.value)
[1] "< 2.22e-16"
Does that answer your question?
KK
On Mon, Oct 3, 2011 at 10:53 AM, Liviu Andronic wrote:
> Dear all
> How does print.htest display the p-value in scientific notation?
>> (x <- cor.test(iris
On 10/02/2011 08:00 PM, Nadine Melhem wrote:
> please send me the "patients.txt" data.
>
> thanks.
>
> __
> 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/post
One point to note, for information, in this discussion is
that cor.test() has apparently returned the P-value as
an exact zero:
x$p.value == 0
# [1] TRUE
identical(x$p.value, 0)
# [1] TRUE
(which, by the way, I was led to after trying log10(x$p.value)
and getting -Inf). Perhaps a more in
On 09/30/2011 04:39 PM, dM/ wrote:
> Suppose that I'm working on Hadley's diamond dataset and I want to
> review the relationship between price, colour and carat.
>
> I might run the following:
>
> library(ggplot2)
>
> #plot scatter and add some hex binning
> q<-qplot(carat,price,data=diamonds
On 09/30/2011 04:39 PM, dM/ wrote:
> Suppose that I'm working on Hadley's diamond dataset and I want to
> review the relationship between price, colour and carat.
>
> I might run the following:
>
> library(ggplot2)
>
> #plot scatter and add some hex binning
> q<-qplot(carat,price,data=diamonds
dear users,
sorry for the late reply! thank you all very much for your input.
i think for now i will try rstudio which was recommended most and later on i
will start learning vim or emacs.
it was very interesting to read about your experiences. thanx!
marion
2011/9/20 Ken Hutchison
> Rstudio and
Thanks all for your pointers. The following does trick:
> base::format.pval(x$p.value) ##Hmisc also has such a function
[1] "<2e-16"
On Mon, Oct 3, 2011 at 10:25 AM, Rolf Turner wrote:
> Isn't it true that 0 < 2.2e-16?
>
Yes, but it doesn't mean that the p-value actually hits absolute zero.
An
Hi
>
> I realize that this is terribly basic, but I just don't seem to see it
at
> this moment, so I would very much appreciate your help.
>
>
> How shall I transform this dataframe:
>
> > df1
> Name Index Value
> 1a 1 0.1
> 2a 2 0.2
> 3a 3 0.3
> 4a 4
On Fri, 2011-09-30 at 12:33 -0400, David Winsemius wrote:
> On Sep 30, 2011, at 2:16 AM, Petr PIKAL wrote:
>
> >>
> >> David,
> >> thank you for your reply
> >>
> >> I tried this
> >> attach(mtcars)
> >> interaction.plot(cyl, gear, mpg, type="o", pch=5:8, lty=1 )
> >>
> >> but I got this error:
>
format.pval is documented and accessible from outside of base.
So you do not have to qualify it as base::format.pval
On 10/3/2011 11:24 AM, Liviu Andronic wrote:
> Thanks all for your pointers. The following does trick:
>> base::format.pval(x$p.value) ##Hmisc also has such a function
> [1] "<2e-1
dataguru gmail.com> writes:
> Based on very fast converging series for special transcendental
> numbers. Is there some R code available?
>
> See details about the RNG at
> http://www.analyticbridge.com/profiles/blogs/new-state-of-the-art-
> random-number-generator-simple-strong-and-fa
(Broken
Noel Faux unimelb.edu.au> writes:
>
> Hi all,
>
> Not sure if this is the best list, please point me to a more appropriate list
if necessary.
>
> Running Mac OSX 10.7.1
> R version 2.13.1 Patched (2011-08-14 r56741)
> Copyright (C) 2011 The R Foundation for Statistical Computing
> ISBN 3-90005
I've re-looked at survexp with the question of efficiency. As it
stands, the code will have 3-4 (I think it's 4) active copies of the X
matrix at one point; this is likely the reason it takes so much memory
when you have a large data set.
Some of this is history; key parts of the code were writt
Hello,
I have exactly the same problem that bibtex is not being called and so the
bibliography is not being processed...
Did you find any solution for that?
Many thanks
syrvn
--
View this message in context:
http://r.789695.n4.nabble.com/How-to-run-Bibtex-with-pdfLatex-in-StatEt-MikTex-on-Wind
Hi everyone,
My (simplified) problem is the following one: I have two tables. The first
table contains 5 columns with 5 values and the second table contains a
single value for each vector name of the first table (see tables below):
Table 1:
A B C D
Hi Nevil,
the function bquote() should do what you want:
(found here: http://r.789695.n4.nabble.com/expression-td904189.html)
mtext(bquote(R^2==.(round(summary(mylm)$r.squared,3))),1)
Regards,
Sina Rüeger
--
View this message in context:
http://r.789695.n4.nabble.com/How-to-format-Rsuperscript
Hi everyone,
I have the following problem:
I have three variables, 'group', 'city' and 'pressure'
There is an interaction effect between group and city and I'd like to show
this in an interaction plot:
interaction.plot(group, city, pressure, type="b",
col= c(1:2),
Dear all,
I have a time series with a frequency of 10 days (so 36 yearly). one year is
completely NA. Now I want to do a stl-decomposition, but using e.g.
na.action= na.approx makes no sense for a whole year, of course. Is there a
way of simulating this single year or to just make stl not using t
On 03/10/2011 10:12 AM, syrvn wrote:
Hello,
I have exactly the same problem that bibtex is not being called and so the
bibliography is not being processed...
Did you find any solution for that?
That sounds like a StatET question. You should get bibtex if you run R
CMD texi2dvi --pdf (or R
Dear all,
I'd like to have a dataframe store information about the units of
the data it contains.
You'll find below a minimal exemple of the way I do, so far. I add a
"units" attribute to the dataframe. But I dont' like the long syntax
needed to access to the unit of a given variable (namely
Hi guys,
Have the following problem computing vectors with pure vector algebra and end
up reverting to recursion or for-looping.
Function my_cumsum calculates a weighted average (W) of ratios (R), but only up
to the given size/volume (v). Now I recurse into the vector (from left to
right) wit
It worked, thank you Jim
Erik
--
View this message in context:
http://r.789695.n4.nabble.com/Keep-ALL-duplicate-records-tp3865136p3867709.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.et
It works, thanks a lot Gabor
Erik
--
View this message in context:
http://r.789695.n4.nabble.com/Find-all-duplicate-records-tp3865139p3867724.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://sta
On Oct 3, 2011, at 9:35 AM, bruno Piguet wrote:
> Dear all,
>
> I'd like to have a dataframe store information about the units of
> the data it contains.
>
> You'll find below a minimal exemple of the way I do, so far. I add a
> "units" attribute to the dataframe. But I dont' like the long s
Try this:
X = 1:5; names(X) = letters[sample(5)]
Y = matrix(1:25, 5); colnames(Y) = letters[1:5]
Y[ , names(X)]
Hope this helps,
Michael Weylandt
On Oct 3, 2011, at 6:20 AM, "Samir Benzerfa" wrote:
> Hi everyone,
>
>
>
> My (simplified) problem is the following one: I have two tables. The
Hi Max,
Thanks for the note. In your last paragraph, did you mean "in
createDataPartition"? I'm a little vague about what returnTrain option
does.
Bonnie
Quoting Max Kuhn :
Basically, createDataPartition is used when you need to make one or
more simple two-way splits of your data. For ex
Hi Bruno,
It sounds like what you want is really a separate class, one that has
stores information about units for each variable. This is far from an
elegant example, but depending on your situation may be useful. I
create a new class inheriting from the data frame class. This is
likely fraught
No, it is an argument to createFolds. Type ?createFolds to see the
appropriate syntax: "returnTrain a logical. When true, the values
returned are the sample positions corresponding to the data used
during training. This argument only works in conjunction with list =
TRUE"
On Mon, Oct 3, 20
sum(ifelse(cumsum(W)<=v, W, 0) * R)
HTH,
David L. Reiner
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of fernando.cabr...@nordea.com
Sent: Monday, October 03, 2011 9:50 AM
To: r-help@r-project.org
Subject: [SPAM] - [R] Matrix/Vect
That works great. Thanks, Rich.
From: Richard M. Heiberger [mailto:r...@temple.edu]
Sent: Wednesday, September 28, 2011 1:13 PM
To: Joseph Boyer
Cc: r-help@r-project.org
Subject: Re: [R] how to get old packages to work on R 2.12.1
Joe,
Most firewall issues can be resolved by entering
setIntern
Alternatively, you can make a separate S3 class for a data frame with this
attribute -- and then have appropriate, print, plot, etc. methods that
access it. That way, your messy syntax has to be written only once in the
method functions.
-- Bert
On Mon, Oct 3, 2011 at 8:15 AM, Joshua Wiley wrot
I know you already got a good answer from Joshua Wiley.
Here is a function I find useful.
Anyone who wants to suggest improvements, please do so!
> my.ls
function(pos=1, sorted=FALSE, mode, class){
.result <- sapply(ls(pos=pos, all.names=TRUE),
function(..x)object.size(eval(as.symbol(..x
Hi:
A small toy example:
fakedata <- data.frame(group = factor(rep(1:3, each = 10), labels =
paste('Therapy', 1:3)),
city = factor(rep(c('Amsterdam', 'Rotterdam'), each = 5)),
pressure = rnorm(30))
with(fakedata, interaction.plot(group, city, pressure, t
Hi,
If you want to take advantage of Josh's example below (using an S4
subclass of data.frame), perhaps you might be interested in taking
advantage of the multitude of useful objects/classes defined in the
bioconductor IRanges package:
http://www.bioconductor.org/packages/release/bioc/html/IRange
Hi Duncan,
you were right. texi2dvi does latex + bibtex. Unfortunately I cannot get it
running.
When I run texi2dvi(file = "path/to/tex/file", pdf=TRUE, quiet=FALSE) then I
get the following error message:
___
Error in texi2
On Mon, Oct 3, 2011 at 10:35 AM, bruno Piguet wrote:
> Dear all,
>
> I'd like to have a dataframe store information about the units of
> the data it contains.
>
> You'll find below a minimal exemple of the way I do, so far. I add a
> "units" attribute to the dataframe. But I dont' like the long
Dear R-list,
I would like to generate an intensity map based on a x * y matrix.
Each point in the matrix should get plotted at the coordinate: x = column / y =
row with
a color-intensity (for example gray-value) based on the actual value of this
point.
Is there a convenient package / functio
I am a bit confused...
I wrote a little shell script called tex2pdf which just calls: texi2pdf
body.tex and when I execute it
in a terminal it runs smoothly. If I type into R: system("/Path/to/tex2pdf")
I get the same error message:
You don't have a working TeX binary (tex) installed anywhere in
See ?image
Uwe Ligges
On 03.10.2011 19:05, Martin Batholdy wrote:
Dear R-list,
I would like to generate an intensity map based on a x * y matrix.
Each point in the matrix should get plotted at the coordinate: x = column / y =
row with
a color-intensity (for example gray-value) based on the
On 03.10.2011 03:29, venerealdisease wrote:
Hi, all,
I am a newbie for [R]
Would anyone help me how to transpose a 3x3x3 array for 1:27
Eg.
A<-array(1:27, c(3,3,3)
What is the logic to transpose it to B<-aperm(A, c(3,2,1))
It simply says third dimension first, second second, and first thir
It seems you don't really know how predict works. If you don't supply
new data, it will only return the least squares fit to the old data,
which is the large data block you saw. Check the first example given
in ?predict to see how this works for new (out of sample) data.
More importantly, use of l
Mr. Disease,
As Uwe points out, the syntax is pretty clear, but it is perhaps worth
mulling over why:
R> identical(A[,,1], t(B[1,,]))
TRUE
to confirm that you understand the function.
Michael Weylandt
PS -- Might I suggest, if you insist on anonymity, a different handle?
I'm answering on my lu
I'm interested in creating a graphic -like- this:
c <- ggplot(mtcars, aes(qsec, wt))
c + geom_point() + stat_smooth(fill="blue", colour="darkblue", size=2,
alpha = 0.2)
but I need to show 2 sets of bands (with different shading) using 5%,
25%, 75%, 95% limits that I specify and where the hea
On Oct 2, 2011, at 18:25 , Mao Jianfeng wrote:
> But, vennerable can not be installed on my Mac book.
Works for me. What are the symptoms?
--
Peter Dalgaard, Professor,
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email: p
Hi,
I know now why it did not work. See here:
http://www.mail-archive.com/r-help@r-project.org/msg19682.html
I started R via the terminal with open -a R and then called the script again
as well as the texi2dvi function.
Both worked. The problem now is I use R within Eclipse and I don't know how
Hi:
I would think that, at least in principle, this should work:
a <- ggplot(mtcars, aes(qsec, wt))
a + geom_point() + stat_smooth(fill="blue", colour="darkblue", size=2,
level = 0.9, alpha = 0.2) +
stat_smooth(fill = 'blue', colour = 'darkblue', size = 2,
Your subject line says Windows, but your error message suggests MacOS.
I think you need to post questions that are answerable if you want an
answer.
Duncan Murdoch
On 03/10/2011 12:48 PM, syrvn wrote:
Hi Duncan,
you were right. texi2dvi does latex + bibtex. Unfortunately I cannot get it
run
Hi,
I need to run a PCoA (PCO) for a data set wich has both positive and negative
values for variables. I could not find any distancecoefficient other than
euclidean distace running for the data set. Are there any other coefficient
works with negtive values.Also I cannot get summary out put (th
Hello,
I'm a complete newbie to R. Spent this past weekend reading The Art of R
Programming, The R Cookbook, the language spec, Wikis and FAQs. I sort-of have
my head around R; the dizzying selection of libraries, packages, etc? Not
really. I've probably missed or failed to understand somethin
Hello,
I am interested in generating a compact letter display for the results
of Tukey HSD tests that contain interaction effects. The 'cld' method in
the 'multcomp' package seems only to work for main effects. Does such a
thing exist already? Thank you for any thoughts,
Josh
--
Joshua Capl
Dear All,
Thank you for the replies to my first thread here:
http://r.789695.n4.nabble.com/global-optimisation-with-inequality-constraints-td3799258.html.
So far the best result is achieved via a penalised objective function. This
was suggested by someone on this list privately. I am still l
Hi,
I have a simple question about importing data, I would be very grateful if
you could help me out.
I have used read.csv(file name, header=T, sep=",") to bring in a csv file I
saved in MS Excel.The problem is I have white spaces in the middle of values
(not in the column names), and this messes
Hi -
I am conducting a meta-analysis and I have a matrix of f-statistics, Ns and
dfs from a series of studies that tested for an interaction in a 2x2 anova.
I'd like to test whether the 2x2 interaction is significant in the
aggregate.
Similarly, I have a matrix of chi-square statistics that I'
The Georgia R School (http://georgia-r-school.org) is a non-profit
educational organization. With faculty from Virginia Commonwealth
University, we are conducting live, interactive, synchronous month-long
online courses on PLS path modeling and R statistical software during
October and November. A
I am using readLines to read a fairly large ASCII file. readLines reads
a fixed number of lines, then other R code processes the data, then
readLines reads the same number of lines again, then other R code
processes the data, then
Sort of like:
conn<-file('filename','r')
for (chunk i
thanks!
On 03.10.2011, at 19:11, Uwe Ligges wrote:
> See ?image
>
> Uwe Ligges
>
>
> On 03.10.2011 19:05, Martin Batholdy wrote:
>> Dear R-list,
>>
>>
>> I would like to generate an intensity map based on a x * y matrix.
>>
>> Each point in the matrix should get plotted at the coordinate:
I'm struggling to find a way of parsing a vector of data in this sort of form:
A,B,C
B,B
A,AA,C
A,B,BB,BBB,B,B
into a matrix (or data frame). The catch is that I don't know a priori how many
entries there will be in each element, nor how many characters there will be.
strsplit(vec,",") gets me
Hello everyone
I have manually compiled directory of gstat in a particular folder of my
Unix system.
I want to install this and am unable to use either of the following two
commands
1. R CMD INSTALL
2. Install.packages
I do not understand how to coax above commands to locate the directory that
i
I am very new to R and have been struggling trying to read a basic ".csv" file
from a password protected site with the following code:
myURL
="http://www.frontierweather.com/degreedays/L15N15PowerRegionAverages_10weeks.txt";
test2=read.table(url(myURL),header=TRUE,sep=",")
A 'data.frame' is re
Hi,
I am trying to find a function to match two data frames of different lengths
for one field only.
So, for example,
df1 is:
Name Position location
francesca A 75
cristina B 36
And df2 is:
location Country
75 UK
56 Austria
And I would like to match on "Location" and the output to be something
More on my previous question ...
I have put in timing statements to try to get a better idea of where the
problem is, like so:
conn<-file('filename','r')
for (chunk in 1:10) {
print(paste('begin read at',date()))
Lines<-readLines(conn,n=25)
print(paste('begin processing at',date()))
On 10/03/2011 01:07 PM, davidyeager wrote:
Hi -
I am conducting a meta-analysis and I have a matrix of f-statistics, Ns and
dfs from a series of studies that tested for an interaction in a 2x2 anova.
I'd like to test whether the 2x2 interaction is significant in the
aggregate.
Similarly, I have
Dear Peter,
I am glad to hearing your reply. That is really nice. Thanks a lot.
###
# (1) the problem of the plot venneuler generated me is sets (A,B,C,D,E,F)
should shared 69604 elements.
# But, it illustrated nothing for me for this 6 sets sharing.
>
Start out with merge():
> df <- merge(df1, df2, all.x=TRUE) # could add by="location" for emphasis
> df
location Name Position Country
1 36 cristinaB
2 75 francescaA UK
You could make make your 'Match' column from is.na(df$Country) if you
knew
Hi,
On Mon, Oct 3, 2011 at 1:54 PM, francy wrote:
> Hi,
>
> I am trying to find a function to match two data frames of different lengths
> for one field only.
> So, for example,
> df1 is:
>
> Name Position location
> francesca A 75
> cristina B 36
>
> And df2 is:
>
> location Country
> 75 UK
> 5
Hi Mike,
On Mon, Oct 3, 2011 at 12:31 PM, Mike Pfeiff wrote:
> I am very new to R and have been struggling trying to read a basic ".csv"
> file from a password protected site with the following code:
>
> myURL
> ="http://www.frontierweather.com/degreedays/L15N15PowerRegionAverages_10weeks.txt"
Hi,
On Mon, Oct 3, 2011 at 11:14 AM, francy wrote:
> Hi,
>
> I have a simple question about importing data, I would be very grateful if
> you could help me out.
>
> I have used read.csv(file name, header=T, sep=",") to bring in a csv file I
> saved in MS Excel.The problem is I have white spaces i
Andrés,
Thank you for your help, but that does not capture what I'm looking for. I need
to be able to control the
shaded bound limits and they need to be coincident.
Tom
On 10/3/11 3:37 PM, Andrés Aragón wrote:
Hi,
Try some like this:
c<- ggplot(mtcars, aes(qsec, mpg, colour=factor(cyl)))
library(sos)
cl <- ???'compact letter display'
cl
This opened a table for me in a web browser identifying 3 help
pages in 'multcomp' and 1 in 'multcompView' that mention 'compact letter
display'. Functions in multcompView will work anything that looks like
a distance or similarity
Enrico,
The idea of looking at the code never occurred to me. That's a great idea which
will come in useful
I'm sure.
In this particular case the using list() or as.matrix() does not solve the
problem, but
those are both good ideas.
Thanks for the reply.
-Original Message-
From: Enric
Hi,
I've assumed that you meant to send this to the R-help list, and not just me.
On Mon, Oct 3, 2011 at 3:54 PM, Mike Pfeiff wrote:
> Sarah, Thanks for the suggestion. Although,
> read.table("http://userid:passw...@my.url/file.csv";) did not work as it
> returned the following:
>
> E
On 10/1/2011 3:03 AM, jim holtman wrote:
try this:
vec1<-
c(4,5,6,7,8,9,10,11,12,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81)
vec2<- c
(1,2,3,4,5,6,7,8,9,10,11,12,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,
If you are using rbind() at each iteration, that can slow things down
greatly. Look up a document called the R Inferno which discusses this
in great detail in circle 2.
Michael Weylandt
On Mon, Oct 3, 2011 at 2:26 PM, Cable, Sam B Civ USAF AFMC AFRL/RVBXI
wrote:
> More on my previous question ..
Well how do you want it be made into a matrix if the rows are all
different lengths? Methinks you are finding this tricky for a
reason...
Michael
On Mon, Oct 3, 2011 at 11:40 AM, Benjamin Wright wrote:
>
> I'm struggling to find a way of parsing a vector of data in this sort of form:
>
> A,B,C
>
One order of the usual coming right up!
1 course of "Why does XXX not work for you?" a la francaise, where XXX
is, in your case, the Euclidean distance. Specifically, any metric
worth its salt (in a normed space) satisfies dist(a,b) = dist(a+c,b+c)
so why are negative values a problem?...
2 side
Hi Doug,
Thanks for taking the time to write a great question.
On Mon, Oct 3, 2011 at 12:23 PM, Douglas Philips wrote:
> Hello,
> I'm a complete newbie to R. Spent this past weekend reading The Art of R
> Programming, The R Cookbook, the language spec, Wikis and FAQs. I sort-of
> have my head
That works a treat Thierry, thanks! I wasn't aware of the plyr package but I
like what it does- I'll put it to use work in the future.
Regards
John
-Original Message-
From: ONKELINX, Thierry [mailto:thierry.onkel...@inbo.be]
Sent: Monday, 3 October 2011 6:36 PM
To: Morrongiello, John
More concisely,
ddply(Orange, .(Tree), transform, scaled = scale(age))
HTH,
baptiste
On 4 October 2011 11:24, wrote:
> That works a treat Thierry, thanks! I wasn't aware of the plyr package but I
> like what it does- I'll put it to use work in the future.
>
> Regards
>
> John
>
> -Origin
I like that one too Baptiste, thanks
-Original Message-
From: baptiste auguie [mailto:baptiste.aug...@googlemail.com]
Sent: Tuesday, 4 October 2011 9:56 AM
To: Morrongiello, John (CMAR, Hobart)
Cc: thierry.onkel...@inbo.be; r-help@r-project.org
Subject: Re: [R] new standardised variable b
Hi,
I have two questions want to ask.
1. If I have a matrix like this, and I want to figure out the rows whose
value in the 3rd column are less than 0.05. How can I do it with R.
hsa-let-7a--MBTD10.5282391972.41E-05
hsa-let-7a--APOBEC10.5078694095.51E-05
hsa-let-7a--PAPOLA0.47
thanks
On Fri, Sep 30, 2011 at 5:25 PM, Francois Pepin <
francois.pe...@sequentainc.com> wrote:
> Hi Jiang,
>
> where did you get that definition of the Benjamini-Hochberg correction?
> That is simply not how it works. You can take a look at the original paper,
> it is available online (
> http:/
You still haven't explained what's wrong with *almost every metric
there is*, but if you want other distance metrics have you considered
those in the package you are using, via the function dsvdis().
Consider, for example:
library(labdsv)
X <- get(data(bryceveg));
X[, sample(NROW(X))] <- (-1)*X[
Comments inline:
On Mon, Oct 3, 2011 at 11:27 PM, dilshan benaragama
wrote:
> Yes I think you did not get my problem.
No, you did not state your problem. I have replied to everything you
have actually included to this point. Admittedly, I have failed to
reply to things you did not say...
> Act
Qian Liu gmail.com> writes:
>
> Hi I need some help with ploting the ROC for K-nearest neighbors. Since KNN
> is a non-parametric classification methods, the predicted value will be
> either 0 or 1.
> It will not be able to test for different cutoff to plot ROC. What is the
> package or functio
I use the following function I stole somewhere. There's probably better ways.
white <- function(x){
x <- as.data.frame(x)
W <- function(x) gsub(" +", "", x)
sapply(x,W)
}
#EXAPLE
dat <- paste(letters," ", " ", LETTERS)
(DAT <- data.frame(dat, dat)) #nasty white spaces
white(DAT) #white spa
Hi,
Try some like this:
c <- ggplot(mtcars, aes(qsec, mpg, colour=factor(cyl)))
c + stat_smooth(aes(group=cyl))+stat_smooth(aes(fill=factor(cyl)))+geom_point()
Andrés AM
2011/10/3, Thomas Adams :
> I'm interested in creating a graphic -like- this:
>
> c <- ggplot(mtcars, aes(qsec, wt))
> c
Dear all,
I am new in R and I have been faced with the following problem, that slows
me down a lot. I am short of ideas to circumvent it. So, any help would be
highly appreciated:
I have 2 dataframes x and y. x is very big (70 million observations),
whereas y is smaller (30 observations).
Al
Yes, I meant to reply to all (sorry still new at asking for help)
1. No, I am not able to open the file when I insert "userID:password@" between
"http://"; and "www"
http://userid:passw...@www.frontierweather.com/degreedays/L15N15PowerRegionAverages_10weeks.txt
(replacing useri
Hi Michael,
Many thanks for your reply.
I am handicapped in 2 different areas; firstly R and most importantly
statistical concepts!
May
I ask you to expand on the 2nd paragraph of your reply below. I was
originally using dynlm but then changed back to lm to avoid
complication.
Also, wha
Hi
I want to do simulation on quasi-binomial distribution with some covariates.
Does anyone have an idea how to do that?
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEA
1 - 100 of 111 matches
Mail list logo
