On Mon, 2018-12-10 at 22:17 +0100, Fatma Ell wrote:
> Dear all,
> I'm trying to use ks.test in order to compare two curve. I've 0 values i
> think this is why I have the follonwing warnings :impossible to calculate
> exact exact value with ex-aequos
>
> a=c(3.02040816326531, 7.95918367346939, 10.
example, the cumulative probability of reaching a point
> outside the cube (u or v or w > A) is 1 however, the bigger cube does not
> exists (because the Q is the reference space). Other words, I feel that we
> extend the space to accommodate any cube of any size! Looks a bit weird to
>
Sorry -- stupid typos in my definition below!
See at ===*** below.
On Tue, 2018-10-23 at 11:41 +0100, Ted Harding wrote:
Before the ticket finally enters the waste bin, I think it is
necessary to explicitly explain what is meant by the "domain"
of a random variable. This is not (though
Before the ticket finally enters the waste bin, I think it is
necessary to explicitly explain what is meant by the "domain"
of a random variable. This is not (though in special cases
could be) the space of possible values of the random variable.
Definition of (real-valued) Random Variable (RV):
Le
I think that one can usefully look at this question from the
point of view of what "NaN" and "NA" are abbreviations for
(at any rate, according to the understanding I have adopted
since many years -- maybe over-simplified).
NaN: Mot a Number
NA: Not Available
So NA is typically used for missing v
Pietro,
Please post this to r-help@r-project.org
not to r-help-ow...@r-project.org
which is a mailing liat concerned with list management, and
does not deal with questions regarding the use of R.
Best wishes,
Ted.
On Sat, 2018-07-14 at 13:04 +, Pietro Fabbro via R-help wrote:
> I will try to b
I've been following this thread, and wondering where it might lead.
My (possibly naive) view of these matters is basically logical,
relying on (possibly over-simplified) interpretaions of "NA" and "NaN".
These are that:
"NaN" means "Not a Number", though it can result from a
numerical calculatio
On Tue, Jul 3, 2018 at 9:25 AM, J C Nash wrote:
>
> > . . . Now, to add to the controversy, how do you set a computer on fire?
> >
> > JN
Perhaps by exploring the context of this thread,
where new values strike a match with old values???
Ted
__
R-hel
On Mon, 2018-06-25 at 09:46 +1200, Rolf Turner wrote:
> Does/should one say "the degrees of freedom is defined to be" or "the
> degrees of freedom are defined to be"?
>
> Although value of "degrees of freedom" is a single number, the first
> formulation sounds very odd to my ear.
>
> I would li
ng of Numbers", covering the
functions ceiling(), floor(), trunc(), round(), signif().
Well worth reading!
Best wishes,
Ted.
On Thu, 2018-05-31 at 08:58 +0200, Martin Maechler wrote:
> >>>>> Ted Harding
> >>>>> on Thu, 31 May 2018 07:10:32 +0100 writes:
>
Well pointed out, Jim!
It is infortunate that the documentation for options(digits=...)
does not mention that these are *significant digits* and not
*decimal places* (which is what Joshua seems to want):
"‘digits’: controls the number of digits to print when
printing numeric values."
On the fa
Apologies for disturbance! Just checking that I can
get through to r-help.
Ted.
__
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-project.org/po
Not necessarily. The R-help archives are publicly accessible,
with a "sort by date" option. So if someone sets up a web=page
monitor which reports back when new messages appear there
(at the bpottom end), then their email addresses are readily
copied (subject to " at " --> "@".
Once they have the
A. On Sat, 2018-03-31 at 15:45 +0200, Henri Moolman wrote:
> Could you please provide help with something from R that I find rather
> puzzling? In the small program below x[1]=1, . . . , x[5]=5. R also
> finds that x[1]<=5 is TRUE. Yet when you attempt to execute while, R does
> not seem to
On Sun, 2018-01-21 at 09:59 +0100, Luigi Marongiu wrote:
> Dear all,
> I have a string, let's say "testing", and I would like to extract in
> sequence each letter (character) from it. But when I use substr() I only
> properly get the first character, the rest is empty (""). What am I getting
> wron
On Thu, 2017-07-20 at 14:33 +0200, peter dalgaard wrote:
> > On 10 Jan 2013, at 15:56 , S Ellison wrote:
> >
> >
> >
> >> I am working with large numbers and identified that R looses
> >> precision for such high numbers.
> > Yes. R uses standard 32-bit double precision.
>
>
> Well, for large
[I unadvertently sent my reply below to Jeremie, instead of R-help.
Also, I havve had an additional thought which may clarify things
for R users].
[Original reply]:
The point about this is that (as Rolf wrote) FALSE & (anything)
is FALSE, provided logical NA is either TRUE ot FALSE but,
because the
gits = 53 binary places.
So this normally "almost" trivial feature can, for such a simple
calculation, lead to chaos or catastrophe (in the literal technical
sense).
For more detail, including an extension of the above, look at the
original posting in the R-help archives for Dec 22, 2
gt; 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-project.org/
>>> posting-guide.html
>>> and provide commented, minimal, self-contained, reprodu
>>
>> __
>> 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-project.org/
>> posting
check out the data.table package, as suggested.
>
> -- Bert
> Bert Gunter
>
> "The trouble with having an open mind is that people keep coming along
> and sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>
tat.uiowa.edu/~luke/R/references/weakfinex.html
>> >>
>> >> As far as I can tell, weakrefs are only available via the C API. Is
>> >> there a way to do what I want in R without resorting to C code? Is
>> >> what I want to do better achieved using something other than weakrefs?
esDemonâ
>>
>> Finally I tried
>>
>>> install.packages("devtools") library(devtools)
>>> install_github("ecbrown/LaplacesDemon")
>>
>> but I am not able to install devtools (for similar reasons). So my
>> questions are:
>
>
> saludos,
> José
---------
E-Mail: (Ted Harding)
Date: 04-Feb-2016 Time: 22:01:42
This message was sent by XFMail
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/list
defaults
> doesn't change that.
>
> I don't _think_ things have changed at my end ...
>
> Steve Ellison
-
E-Mail: (Ted Harding)
Date: 04-Feb-2016 Time: 17:03:37
This message was sent by XFMail
_
lpful replies could be sent.
>
> A milder alternative is to encourage some R-help subscribers to click the
> "Don't send" or "Save" button and think better of their replies.
>
>
> --
> Michael Friendly Email: friendly AT yorku DOT ca
> Professor, Psyc
>> >> >> Ista
>>> >> >>
>>> >> >> On Fri, Oct 30, 2015 at 9:15 PM, Val
>>wrote:
>>> >> >>> Hi all,
>>> >> >>> Iam trying to change character to numeric but have probelm
>&
p
Towards the bottom of this page is a section "Subscribing to R-help".
Follow the instructions in this section, and it should work!
Best wishes,
Ted.
---------
E-Mail: (Ted Harding)
Date: 14-Oct-2015 Time: 19:34:55
to generate 1000 numbers from N(u, a^2), however I don't
> want to include 0 and negative values. How can I use beta distribution
> approximate to N(u, a^2) in R.
>
> Thx for help
-
E-Mail: (Ted Harding)
Date: 15-Sep-2015 Time: 16:12
ewhere
in the spreadsheet? (Excel is notorious for planting things invisibly
in its spreadsheets which lead to messed-up results for no apparent
reasion ... ).
Hoping this helps,
Ted.
-
E-Mail: (Ted Harding)
Date: 09-Feb-2015 Time: 22:15:44
This me
people object to code "clutter" from parentheses that could
be more simply replaced (e.g. "var< -4" instead of "var<(-4)"),
but parentheses ensure that it's right and also make it clear
when one reads it.
Best wishes to all,
Ted.
---
Sorry, a typo in my reply below. See at "###".
On 12-Jan-2015 11:12:43 Ted Harding wrote:
> On 12-Jan-2015 10:32:41 Erik B Svensson wrote:
>> Hello
>> I've got a problem I don't know how to solve. I have got a dataset that
>> contains age intervals
t;.
Implementing the above as a procedure:
agegrp[max(which(cumsum(y1994)/sum(y1994)<0.5)+1)]
# [1] "55-64"
Note that the "obvious solution":
agegrp[max(which(cumsum(y1994)/sum(y1994) <= 0.5))]
# [1] "45-54"
give
I should have added an extra line to the code below, to complete
the picture. Here it is (see below line "##".
Ted.
On 11-Jan-2015 08:48:06 Ted Harding wrote:
> Troels, this is due to the usual tiny difference between numbers
> as computed by R and the numbers that yo
_
> 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-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
--
something special about the numbers around
>>> 1000.
>>>
>>> Back to the quesion at hand: I can avoid use of round() and speed things
>>> up
>>> a little bit by just adding a small number after multiplying by 1000:
>>>
>>>> ptm &
l data for 200
> people over one co-variate. So I was hoping instead of completely removing
> the rows, to just somehow acknowledge that the data for this particular
> co-variate is missing in the model but not completely remove the row? This
> is more what I was hoping someone would know if it
>
> [[alternative HTML version deleted]]
>
> __
> 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
>>>>>> options(latexcmd='pdflatex')
>>>>>>> options(dviExtension='pdf')
>>>>>>>
>>>>>>> ## Macintosh
>>>>>>> options(xdvicmd='open')
>>>>>>>
>>>>>>> ## Windows
wer to your question as
>> abs(qt(0.408831/2, df=1221)), but you'll get 4.117.
>>
>> Duncan Murdoch
>>
>>
>>
>
> __
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailma
17mother
107 09sibling
107 18father
107 19mother
108 16sibling
108 NAfather
108 NAmother
109 17sibling
109 NAfather
109 NAmother
That's the data. Now a litt
On 12-Aug-2014 22:22:13 Ted Harding wrote:
> On 12-Aug-2014 21:41:52 Rolf Turner wrote:
>> On 13/08/14 07:57, Ron Michael wrote:
>>> Hi,
>>>
>>> I would need to get a clarification on a quite fundamental statistics
>>> property, hope expeRts here would
cation).
The important thing when using pre-programmed functions is to know
which is being used. R uses (n-1), and this can be found from
looking at
?sd
or (with more detail) at
?cor
Ron had assumed that the denominator was n, apparently not being aware
that R
d) %*% Data_Normalized)/(dim(Data_Normalized)[1]-1)
and compare the result with
cor(Data)
And why? Look at
?sd
and note that:
Details:
Like 'var' this uses denominator n - 1.
Hoping this helps,
Ted.
-
E-Mail: (Ted Harding)
Date: 12-Aug-2014 Time: 22:32:26
This message w
4*a*b
MEAN^2 - 3*SD^2 = a*b
Hence for a >= 0 and b > a you must have MEAN^2 >= 3*SD^2.
Once you have MEAN and SD satisfying this constraint, you should
be able to solve the equations for a and b.
Hoping this helps,
Ted.
-
E-Mail: (Ted Hardi
gt; }
>>
>> It works, but words cannot describe its ugliness. And it gets really
>> slow really fast with growing x.
>>
>> So, anyway. My two questions are:
>> 1. Does R really, really, seriously lack a permutation function?
>> 2. OK, stop kidding me. So what
I tried to google and found 'tempdir()'.
>>
>> I checked my temporary file but it was empty.
>>
>> Just in case, after running 'unlink(tempdir(),recursive=TRUE)', I
>> restarted my computer, but memory space did not change. But there
>> seems
___
> 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
elapsed
# 0.028 0.000 0.029
system.time(for(i in (1:1)) (1)*2)*3)*4)*5) )
# user system elapsed
# 0.052 0.000 0.081
(though in fact the times are somwhat variable in both cases,
so I'm not sure of the value of the relationship).
Best wishes,
Ted.
---
159") * as.bigz("36662978")
> Big Integer ('bigz') :
> [1] 6195624596620653502
>
>
> On Sun, May 4, 2014 at 12:50 PM, Ted Harding
> wrote:
>> On 04-May-2014 14:13:27 Jorge I Velez wrote:
>>> Try
>>>
>>> options(digits =
t of 168988580159 * 36662978 in my algorithm and r product.
>> The Microsoft calculator confirm my number.
>>
>> Thanks.
>> --
>> *Gráfica ARTENTOR *
>>
>> de Diego L. Tentor
>> Echagüe 558
>> Tel.:0343 4310119
>> Paraná - En
ven simpler (if it is only one particular month you want,
as in your example) is:
$ cal April 2014
which yields:
April 2014
Mo Tu We Th Fr Sa Su
1 2 3 4 5 6
7 8 9 10 11 12 13
14 15 16 17 18 19 20
21 22 23 24 25 26 27
28 29 30
and now just count down the 3rd col
ailman/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-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
&g
hat the
boundary is drawn as a set of separate partial boundaries which
are in no particular order as a whole; and in some datasets the
different separate parts of the boundary do not exactly match up
at the points where they should exactly join.
Hoping this helps,
Ted.
--------
t integer when not exactly halfway between, and rounds
either always up or always down when the fractional part is exactly 1/2,
then I think (but others will probably correct me) that you may have
to write your own -- say roundup() or rounddown():
roundup <- function(x){
if((x-floor(x))==
aa0<-gsub("^[0-9]+ ","",aa)
aa0
# [1] "(472)" "(445)" "(431)" "(431)" "(415)" "(405)" "(1)"
as.numeric(gsub("[()]","",aa0))
# [1] 472 445 431 431 415 405 1
Ted.
x); so the actual byte content
of newseed is:
4b e9 76 34 41 cf 5e 17 b0 68 78 98 87 9e 8b 5f
fb 4f 52 e6 59 ef 0b 58 52 58 4a 3a df 04 c1 8d
This could be achieved via a system() call from R; and the contents
of newseed would then need to be converted into a format suitable
for use as ar
- (1-p1)*(1-p2)* ... *(1-pk)
where pj is P(Aj). Hence
punion <- function(p){1 - prod(1-p)}
should do it!
Ted.
---------
E-Mail: (Ted Harding)
Date: 18-Feb-2014 Time: 23:51:31
This message was sent by XFMail
_
sed the problem!
Ted.
> This generates a distribution with origin equal zero, but I want the
> origin to be x0
>
> How can I handle shifted gamma distribution in r?
>
> Thanks a lot!
> Rodrigo.
-
E-Mail: (Ted Harding)
Date: 1
com/
>>>> - involve R Core to give blessing for using the R logo, if necessary.
>>>> This would be similar to what Ubuntu does with AskUbuntu:
>>>> http://meta.askubuntu.com/questions/5444/is-ask-ubuntu-official-ubuntu
>
>>
>
> __
> 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, r
to give other profiling information,
> including more detail on execution position (down to the line number),
> and various measures of memory use.
>
> Duncan Murdoch
Also have a look at
?trace
which you may be able to use for what you want.
Ted.
-
E-Mail: (Ted Harding)
Date
nabble.com/Why-does-sin-pi-not-return-0-td4676963.html
>>
>> There is an example by Martin Maechler (author of Rmpfr) on how to use
>> arbitrary precision
>> with your arithmetic.
>>
>> On 22 December 2013 10:59, Ted Harding wrote:
>>> Greetings All!
>
R is off its head!
All it has to do is multiply by 2 -- and it gets it cumulatively wrong!
R just doesn't add up ...
Season's Greetings to all -- and may your calculations always
be accurate -- to within machine precision ...
Ted.
-
E-Mail
-1)] + x[2:N]
# [1] 5 13 29 23 10
Best wishes,
Ted.
-------------
E-Mail: (Ted Harding)
Date: 14-Dec-2013 Time: 10:54:00
This message was sent by XFMail
__
R-help@r-project.org mailing list
https://stat.ethz.ch/ma
4: "+XZZZU.C5BF89ZZZUBP+"
5: "+XZZZU.CZUZUBF89ZZZUBP+"
6: "+XZZZU.CZUZUBF89ZZZUBP+"
7: "+XZZZU.CZUZUBF89ZZZUBP+"
8: "+XZZZU.CZUZUBFUZZZ9ZZZUBP+"
9: "+XZZZU.CZUZUBFUZZZUZZUZZZUBP+"
A: "
_____________
> 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.
1
# [2,]112
# [3,]113
# [4,]122
# [5,]123
# [6,]133
# [7,]222
# [9,]233
#[10,]333
There may be a simpler way!
Ted.
-
E-Mail:
qnorm(0.05155075)
[1] -1.63
so maybe you mistyped "1.63" instead of "1.53"?
-----
E-Mail: (Ted Harding)
Date: 16-Oct-2013 Time: 16:12:56
This message was sent by XFMail
__
R-help@r-projec
ave missed the essential point of
your query, but worth a try ... ).
I have no idea how to achieve a similar effect in Windows ...
Ted.
-
E-Mail: (Ted Harding)
Date: 29-Sep-2013 Time: 19:31:29
This message was sent by XFMail
__
cheers,
> Rolf Turner
Though, mind you, FAQ 3.71 does also offer some consolation to R:
all.equal(0,sin(pi))
# [1] TRUE
So it depends on what you mean by "different from". Computers
have their own fuzzy concept of this ... Babak has too fussy
a concept.
Ted.
-
On 21-Aug-2013 19:08:29 David Winsemius wrote:
>
> On Aug 21, 2013, at 10:30 AM, (Ted Harding) wrote:
>
>> Greetings all.
>>
>> I suspect this question has already been asked. Apologies
>> for not having taced it ...
>>
>> In the default pairs
-Mail: (Ted Harding)
Date: 21-Aug-2013 Time: 18:30:44
This message was sent by XFMail
__
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
proportions of the populatipn:
50th to 85th = 35%; 31st to 69th = 38%; 69th to 93rd = 24%. So you are
still facing issues of what you mean, or what you want to mean.
Simpler to stick to the original "odds per unit of x" and then apply
it to whatever multiple of the unit you happen to be int
Don't worry about it. As I say, it can happen to anyone (though more
often to some than to others). If it is a proper message to R-help,
one of the moderators will approve it (though quite possible not
immediately).
Hoping this helps,
Ted (one of the mode
18 15 22 29
# [2,]29 16 23 30
# [3,]3 10 17 24 31
# [4,]4 11 18 25 32
# [5,]5 12 19 26 33
# [6,]6 13 20 27 34
# [7,]7 14 21 28 35
# To permute the rows:
t(app
eak)
choose(37,7)/choose(40,10)
# [1] 0.01214575
so the chance of all 3 being in some one of the 4 groups is
4*choose(37,7)/choose(40,10)
# [1] 0.048583
which, if you are addicted to P-values, is just significant
at the 5% (P <= 0.05) level. So this gives some indication
that the &q
gt;>
>> >From L, I wish to obtain (as directly as possible, e.g. avoiding
>> a loop) two vectors each of length N where one contains the strings
>> that are first in the pair, and the other contains the strings
>> which are second, i.e. from L (as above) I would want to extra
trings
which are second, i.e. from L (as above) I would want to extract:
V1 = c("A1","A2","A3",...)
V2 = c("B1","B2","B3",...)
Suggestions?
With thanks,
Ted.
-
E-Mail: (Ted Harding)
ng from which you can extract the
individual digits. And then on to whatever you want to do ...
Hoping this helps,
Ted.
-
E-Mail: (Ted Harding)
Date: 18-Apr-2013 Time: 10:06:43
This message was sent by XFMail
_
__
> 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
is a vector, it will be promoted to either a row or column matrix
to make the two arguments conformable. If both are vectors it
will return the inner product (as a matrix).
Usage:
x %*% y
[etc.]
Ted.
---------
E-Mail: (Ted Harding)
Date: 11
es"
will be the result of a calculation) then one useful precaution
could be to round the result:
round(0.29*100)
# [1] 29
29-round(0.29*100)
# [1] 0
length(rep(TRUE,0.29*100))
# [1] 28
length(rep(TRUE,round(0.29*100)))
# [1] 29
(The default for round() is 0 decimal places, i.e.
On 01-Apr-2013 21:26:07 Robert Baer wrote:
> On 4/1/2013 4:08 PM, Peter Ehlers wrote:
>> On 2013-04-01 13:37, Ted Harding wrote:
>>> Greetings All.
>>> This is a somewhat generic query (I'm really asking on behalf
>>> of a friend who uses R on Window
Chunk 1? (The size-change may perhaps
have to be determined empirically).
With thanks,
Ted.
-------------
E-Mail: (Ted Harding)
Date: 01-Apr-2013 Time: 21:37:17
This message was sent by XFMail
__
R-help@r-project.org mail
p
# 0.4957627
So it doesn't do the requested continuity correction in [A] because
there is no need to. But in [B1] it makes a difference (compare
with [B2]), so it does it.
Hoping this helps,
Ted.
-
E-Mail: (Ted Harding)
Date: 27-Mar-2013
-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf Of Ted Harding
> Sent: Dienstag, 26. März 2013 11:09
> To: r-help@r-project.org
> Subject: [R] edit.data() read-only?
>
> Greetings All.
>
> The function edit.data() allows a convenient spread
s, Barry.
Ted.
On 26-Mar-2013 10:20:59 Barry Rowlingson wrote:
> On Tue, Mar 26, 2013 at 10:09 AM, Ted Harding
> wrote:
>> Greetings All.
>>
>> The function edit.data() allows a convenient spreadsheet-like
>> view of a dataframe with too many rows/columns to fit on
Or some other
function which could offer similar viewing capability without
the risk of data change?
With thanks,
Ted.
-----
E-Mail: (Ted Harding)
Date: 26-Mar-2013 Time: 10:08:58
This message was sent by XFMail
_
sing order.
Is that what you mean?
Hoping this helps,
Ted.
---------
E-Mail: (Ted Harding)
Date: 22-Mar-2013 Time: 13:31:31
This message was sent by XFMail
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-hel
imilar for other arbitrary choices of first and second distribution
(so long as each has at least a second moment, hence excluding, for
example, the Cauchy distribution).
That's about as far as one can go with your question!
Hoping it helps, howevr.
Ted.
---
# [1] FALSE
(0.1 + 0.05) < (0.15 - .Machine$double.eps^0.5)
# [1] FALSE
(or similar). Observe that
.Machine$double.eps^0.5
# [1] 1.490116e-08
.Machine$double.eps
# [1] 2.220446e-16
(0.1 + 0.05) - 0.15
# [1] 2.775558e-17
Hoping this helps,
Ted.
---
The size of the first set
area2 The size of the second set
cross.area The size of the intersection between the sets
so Ben has hit the nail on the head! (Maybe "area2 = 325" was a typo?).
Best wishes to all,
Ted.
-
E-Mail: (Ted Har
> rv.t.pvalues <- array(NA, k)
>
> for(i in 1:k){
> rv.t.pvalues[i] <- t.test(rv[i, 1:(c/2)], rv[i, (c/2+1):c],
> equal.var=TRUE, alternative="two.sided")$p.value
> }
>
> hist(rv.t.pvalues)
>
> The histogram is this one:
> *http://tinyurl.com/histogram
eason why you are observing
this may become more evident!
Hoping this helps,
Ted.
-------------
E-Mail: (Ted Harding)
Date: 09-Jan-2013 Time: 10:29:21
This message was sent by XFMail
__
R-help@r-project.org mailing
Sorry, I made a mistake in re-writing your code below. See at [***]
On 18-Dec-2012 21:00:28 Ted Harding wrote:
On 18-Dec-2012 20:09:36 Beatriz González Domínguez wrote:
> Hello,
>
> I have done a scatterplot and now would like to add its regression
> line but it does not show.
>
abline(lm3)
whereas it should be
Either:
lm3 <- lm(V~U)
plot(U,V)
abline(lm3)
Or:
lm3 <- lm(U~V)
plot(V,U)
abline(lm3)
Hoping this helps,
Ted.
-
E-Mail: (Ted Harding)
Date: 18-Dec-2012 Time: 21:00:25
This message was sent by
ared = 0.0556, df = 1, p-value = 0.8136
which agrees with your calculation of
sum((x$observed-x$expected)^2/x$expected)
# [1] 0.05562457
Hoping this helps,
Ted.
-
E-Mail: (Ted Harding)
Date: 03-Dec-2012 Time: 22:44:14
This message was sent by XFMai
served table, the p-value being
the sum of such probabilities.
I hope this helps.
Ted.
-
E-Mail: (Ted Harding)
Date: 03-Dec-2012 Time: 22:24:00
This message was sent by XFMail
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo
<- 0.25+0.5*(0:20)
lines(x0,dchisq(x0,1))
Hoping this helps,
Ted.
-
E-Mail: (Ted Harding)
Date: 02-Dec-2012 Time: 15:02:45
This message was sent by XFMail
__
R-help@r-project.org mailing list
https:
1 - 100 of 1019 matches
Mail list logo
