Hi,
it seems that the directions in
http://wiki.r-project.org/rwiki/doku.php?id=getting-started:installation:eeepc
refer to installing R on a eeepc 701 and actually
Yes, it does :)
fail for a 901 at the first step.
What exactly failed? You are actually do not need the KDE desktop.
Greeti
Thanks Gavin and sorry to all for this unnecessary question.
B.
Gavin Simpson wrote:
>
> On Tue, 2008-09-16 at 10:47 -0700, Birgitle wrote:
>> Hello R-User!
>>
>> I try to do the following:
>>
>> New<-iris[c(1:7,90:97),1:5]
>> New.rpart<-rpart(Species~., data=New, method="class")
>>
>> New.
Agustin,
which EEE PC 901 is it? The one for Windoze or Linux (Xandros)?
el
on 9/16/08 6:50 PM Agustin Lobo said the following:
> Hi!
>
> it seems that the directions in
> http://wiki.r-project.org/rwiki/doku.php?id=getting-started:installation:eeepc
>
>
> refer to installing R on a eeepc 701
If you have read the table as a data frame say, DF then the following should
work
Â
DF[DF==0]<-NA
--- On Tue, 16/9/08, Amit Patel <[EMAIL PROTECTED]> wrote:
From: Amit Patel <[EMAIL PROTECTED]>
Subject: [R] HI
To: r-help@r-project.org
Date: Tuesday, 16 September, 2008, 7:02 PM
Does anyone kno
Steven McKinney wrote:
> Not likely that anyone can explain, as
> there is not enough information in your
> email.
>
> Including the contents of the freqtest.txt file
> was a good idea, as the posting guide suggests
> (the posting guide is that clearly labeled bit
> at the bottom that looks like th
M. wrote:
Hello,
I have a matrix A with value varying from -1 to 1. I hope to use scaled
color based on its value to produce an image of this matrix.
Suppose I hope to label those data in [-1,-0.5] with blue, label those
[-0.5,0.8] with light blue (tone is proportional to its value, e.g.
rgb(1
hi, useRs-
i am using R with emacs/ESS (on xubuntu).
currently, when i send a large script to ESS process,
i have to wait for the whole script to finish.
is there a way to immediately return to the editing window
after submitting? if there is, how do i set it permanently?
thank you in advanc
Hello,
I am trying to find a permutation test that works on a general nxm table. The
data set is small enough to have cells with too small counts to make
chi2-approximation invalid. If the table was a 2x2 contingency table I would
like to use a Fsher exact test (fisher.test) but that wont work
Greetings everyone,
I'm trying to add a specific table or a specific number of rows (e.g.44) to a
table with no success.
This is my basic table
> head(dataA)
yearplot spp prop.BDCA1DCA2DCA3DCA4
1 20001 a1 0.031079 -0.0776 -0.0009 0.0259 -0.0457
2 2000
Hello,
I'm trying to plot graphs using lattice with this script :
xyplot(Y ~ X | factmod,
panel = function(x, y) {
panel.grid(h=-1, v=-1, col="gray")
panel.xyplot(x, y, type="p", pch=20)
panel.points(50,Idata, data=devdata, col="red") ->
this
Using r2dtable() you can simulate general tables nxm with given margins. Based
on these you acn calculate a Monte Carlo p-value for a conditional test for
independence.
Regards
Søren
-Oprindelig meddelelse-
Fra: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] På vegne af Magnus Petterss
In the help for scale_fill_identity, it is written
## Not run:
colour <- c("red","green","blue","yellow")
qplot(1:4, 1:4, fill=colour, geom="tile")
qplot(1:4, 1:4, fill=colour, geom="tile") + scale_fill_identity()
# To get a legend, you also need to
Greetings everyone,
I'm trying to add a specific table or a specific number of rows (e.g.44) to a
table with no success.
This is my basic table
> head(dataA)
yearplot spp prop.BDCA1DCA2DCA3DCA4
1 20001a1 0.031079 -0.0776 -0.0009 0.0259 -0.0457
2 20
Rodrigo Aluizio gmail.com> writes:
>
> Thanks for the elucidation James.
> After reading the .pdf of proj4 and your answer, I believe I won't need
> proj4 package.
> Well, I still don't know all the kind of transformations we will need here.
> I'm just getting prepared to future data that is a
Jason Gasper noaa.gov> writes:
>
> Here is the WinBUGS code
>
>
> Uwe Ligges wrote:
> > If you want that people help to translate *code*, you have to specify
> > it ...
> >
> > Uwe Ligges
Please look at the nb2WB() function in the spdep package, I think that you'll
find that it provides wh
On 17-Sep-08 05:27:40, Alexy Khrabrov wrote:
> Exactly -- also found creating horizontal vector helps:
>
> > df <- data.frame(matrix(1:5,nrow=1))
> > colnames(df) <- LETTERS[1:5]
> > df
>A B C D E
> 1 1 2 3 4 5
>
> Thanks,
> Alexy
Another (though similar) approach (e.g. if you want to bui
Hi
I have a data frame which looks like the following (it is a table of commit
logs to a source repo). The Author column
is treated as a factor, by virtue of stringsAsFactors being set to TRUE, so I
can extract the number of commits
per author by using table(logs$Author).
Is there an elegant wa
[EMAIL PROTECTED] wrote:
> Hi
>
> I have a data frame which looks like the following (it is a table of commit
> logs to a source repo). The Author column
> is treated as a factor, by virtue of stringsAsFactors being set to TRUE, so I
> can extract the number of commits
> per author by using table
Søren Højsgaard wrote:
> Using r2dtable() you can simulate general tables nxm with given margins.
> Based on these you acn calculate a Monte Carlo p-value for a conditional test
> for independence.
>
> Regards
> Søren
>
fisher.test(, simulate.p.value=TRUE) might be more direct. Also
work
[My previous message rejected, therefore I am sending same one with some
modification]
I have 3 vectors with object name : dat1, dat2, dat3
Now I want to create a loop, like :
for (i in 1:3)
{
cat(sd(dati))
}
How I can do this in R?
Regards,
Try this:
sapply(ls(patt="^dat[0-9]"), function(x)sd(get(x)))
On Wed, Sep 17, 2008 at 8:50 AM, Megh Dal <[EMAIL PROTECTED]> wrote:
> [My previous message rejected, therefore I am sending same one with some
> modification]
>
> I have 3 vectors with object name : dat1, dat2, dat3
>
> Now I want to
you need to get(), e.g., try this:
dat1 <- rnorm(5)
dat2 <- rnorm(6)
dat3 <- rnorm(7)
lis <- lapply(paste("dat", 1:3, sep = ""), get)
lis
sapply(lis, sd)
I hope it helps.
Best,
Dimitris
Megh Dal wrote:
[My previous message rejected, therefore I am sending same one with some
modification]
saggak wrote:
Dear Mr Frank,
I thank you for your prompt reply. However, I am not able to understand
(may be since for me R is a new venture) the contents of your reply. If
its a book you are referring to, I don't have access to it. How do I get
@ARTICLE{hos97com and how do I run it in R?
Hi
I'm trying to make a 3d plot showing a point cloud, the corresponding
data ellipse
and the principal axes of the ellipse as vectors.
library(rgl)
data(trees)
cov <- cov(trees)
mu <- mean(trees)
plot3d(trees, type="s", size=0.5, col="blue", cex=2)
In this step, an extra box is added. I've
Hi Michael,
You need to specify both 'box' and 'axes' to FALSE to avoid the box
lines (if you don't specify the latter one, there will still be axes
lines).
Best,
Yihui
--
Yihui Xie <[EMAIL PROTECTED]>
Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086
Mobile: +86-15810805877
Homepage: http://www.
Colleagues,
I am using R 2.7.2 on all three platforms (OS X; Windows; Linux) and I
have encountered the following problem:
I use R to issue a system call the execute a command in the OS. For
example:
system("DOSOMETHING TOSOMEFILE")
In most situations, this task completes successfully and
Hi all,
I am generally quite fond of the unix commandline keystrokes (e.g. searching
back in your history with [CTRL]-R, and cutting/pasting with [CTRL]-K/Y)
which work in the R commandline in *nix. Does anyone know if there's any
way to get similar functionality in the Windows RGUI?
I know tha
Yihui Xie wrote:
Hi Michael,
You need to specify both 'box' and 'axes' to FALSE to avoid the box
lines (if you don't specify the latter one, there will still be axes
lines).
Hi Yihui
I tried several variations, none of which would simply add the ellipse
and nothing else to the original
p
No, I mean in plot3d():
plot3d(trees, type="s", size=0.5, col="blue", cex=2, box=FALSE, axes=FALSE)
not in ellipse3d()
Yihui
--
Yihui Xie <[EMAIL PROTECTED]>
Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086
Mobile: +86-15810805877
Homepage: http://www.yihui.name
School of Statistics, Room 1037,
Why not use a script?
I feel that it is much better than using the history via [CTRL]-R in
unix, which also pulls up errorneous commands.
A script is vital for statistical analysis and research where you may
want to or be asked to repeat or reproduce the analysis months later.
Rgui (on wind
Is there a way to use the cycle variable for rowname?
v=1:6
for (a in 1:3){
for (b in 4:5) {
v=rbind(v,a.b=1)
}
}
v
This above obviously does not work, but I couldn't find out how to use a and
b to construct a rowname like 14, 15, 24, 25.
Thanks for the help.
Balazs
--
View this message in
> Is there a way to use the cycle variable for rowname?
>
> v=1:6
> for (a in 1:3){
> for (b in 4:5) {
> v=rbind(v,a.b=1)
> }
> }
> v
>
> This above obviously does not work, but I couldn't find out how to use a
and
> b to construct a rowname like 14, 15, 24, 25.
Not pretty, but this does t
On 17-Sep-08 14:22:11, Ralikwen wrote:
>
> Is there a way to use the cycle variable for rowname?
>
> v=1:6
> for (a in 1:3){
> for (b in 4:5) {
> v=rbind(v,a.b=1)
> }
> }
> v
>
> This above obviously does not work, but I couldn't find out how to use
> a and
> b to construct a rowname like 1
Hi, I'm a new user of R. My background is Electrical Engineering, so
please bear with me if this is a silly question.
I'm trying to assess whether the results of an experiment satisfy the
hypothesis of homoscedasticity (my ultimate goal is to use ANOVA).
The result of the experiment is mean delay
Emmanuel,
Please see inline for comments.
Le mer. 17 sept. à 00:59, Emmanuel Levy a écrit :
Dear list members,
I encountered this problem and the solution pointed out in a previous
thread did not work for me.
(e.g. install.packages("RCurl", repos = "http://www.omegahat.org/R";)
I work with
I wonder if anyone is familiar with the priprobit software (Sakuma M.
1998. Probit analysis of preference data. Appl Entomol Zool.
33:339–34). I would like to replicate the same analysis and output in
R. Of note, I want to determine the analytical sensitivity, defined as
95% cut-off value.
j
Hi there,
I have a dataframe D.F1
>dim (D.F1)
14351 9
This dataframe has values and for some 1000 rows it holds NULL values.I hace
found the missing values for about 500 and have those in another dataframe
D.F.sub.2
>dim(D.F.sub.2)
500 9
as dataframe is a subset of D.F1 the coulmn 1 in D.F.sub
On 9/17/2008 9:40 AM, Michael Friendly wrote:
Yihui Xie wrote:
Hi Michael,
You need to specify both 'box' and 'axes' to FALSE to avoid the box
lines (if you don't specify the latter one, there will still be axes
lines).
Hi Yihui
I tried several variations, none of which would simply add t
> I'm trying to plot graphs using lattice with this script :
>
> xyplot(Y ~ X | factmod,
> panel = function(x, y) {
> panel.grid(h=-1, v=-1, col="gray")
> panel.xyplot(x, y, type="p", pch=20)
> panel.points(50,Idata, data=devdata, col="red") ->
On 9/17/2008 9:03 AM, Dennis Fisher wrote:
Colleagues,
I am using R 2.7.2 on all three platforms (OS X; Windows; Linux) and I
have encountered the following problem:
I use R to issue a system call the execute a command in the OS. For
example:
system("DOSOMETHING TOSOMEFILE")
In most si
Is there a setting change (or other minor fix) that can be done without
an install of either the patched or development versions to address the
printing issue described below. It is my understanding that these
versions are 'source code' and not compiled for 'easy' installation.
Thanks
Tom
---
Hi Tom,
No, there are binaries for both the patched and development versions.
On the page on CRAN where you normally download the windows binary
(e.g., http://cran.wustl.edu/bin/windows/base/), there are links at
the top to the patched and devel builds.
HTH,
Jenny
At 10:31 AM 9/17/2008, To
I have a dataframe D.F1
>dim (D.F1)
14351 9
This dataframe has values and for some 1000 rows it holds NULL values.I hace
found the missing values for about 500 and have those in another dataframe
D.F.sub.2
>dim(D.F.sub.2)
500 9
as dataframe is a subset of D.F1 the coulmn 1 in D.F.sub.2 is
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Hi Emanuel.
The best thing to do is upgrade to R-2.7.2 (or any 2.7.*) and the
problem should disappear. It relates to encoding of strings.
D.
Emmanuel Levy wrote:
> Dear list members,
>
> I encountered this problem and the solution pointed ou
Others have told you about order, also note that order(order(x)) gives the
indexes to return a sorted x variable back to its original order.
Hope this helps,
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
> project.org] On Behalf Of Peng Jiang
> Sent: Tuesday,
Hi,
I would be extremely grateful if somebody can help with following
problem.
We have done an experiment using a reference design with technical and
biological replicates where the target file looks like this:
FileName Cy3 Cy5
xxx A REF
xxx
For purely educational purposes I would like to see the spectrum of a kernel. I
tried:
spectrum(kernel("daniell", 100), taper=0, log="no")
but I get:
Error in complete.cases(object) : not all arguments have the same length
I can plot the kernel but that is only in the time-domain. I would like
After you have interpolated your data, replace the values in the areas of the
ocean with NA, then filled.contour will not plot in those areas.
Hope this helps,
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
> project.org] On Behalf Of Victor Homar
> Sent: Tuesda
On Wed, Sep 17, 2008 at 5:34 AM, Pedro Barros <[EMAIL PROTECTED]> wrote:
>
> In the help for scale_fill_identity, it is written
> ## Not run:
> colour <- c("red","green","blue","yellow")
> qplot(1:4, 1:4, fill=colour, geom="tile")
> qplot(1:4, 1:4, fill=colour, geom="tile")
Hi all,
though i know this is a simple question, i really hope someone could.
I am just trying to compare percents respectively the difference in
percents by a simple test. t.test usually compares means and is
working like it is supposed to be. Now i wonder how i can use a test
in R to test
On Sep 17, 2008, at 9:26 AM, mfrumin wrote:
Hi all,
I am generally quite fond of the unix commandline keystrokes (e.g.
searching
back in your history with [CTRL]-R, and cutting/pasting with [CTRL]-
K/Y)
which work in the R commandline in *nix. Does anyone know if
there's any
way to get
Sounds like a homework problem. Consult any basic statistics text. Hint:
You don't use a t-test to test for a difference in percentages. Check
"chi-squared" and/or "contingency table."
Cheers,
Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of B
Adaikalavan, thanks.
Perhaps I was not so specific enough in what I want, for those not so
familiar with unix commandline featuers. I'm looking for the 'reverse
search' functionality where you hit CTRL-R, then start typing a bit of text
and it finds previous commands with that bit of text, which
Look at GNU Emacs with ESS, you'll get what you want out of the box
--- and much more! See
http://vgoulet.act.ulaval.ca/en/emacs
[shameless plug] for a distribution of Emacs bundled with ESS.
HTH
Vincent
Le mer. 17 sept. à 09:26, mfrumin a écrit :
Hi all,
I am generally quite f
On Sep 17, 2008, at 12:41 PM, Bunny, lautloscrew.com wrote:
Hi all,
though i know this is a simple question, i really hope someone could.
I am just trying to compare percents respectively the difference in
percents by a simple test. t.test usually compares means and is
working like it is s
Dear all,
The documentation on the function open3d in the rgl package, reads:
open3d opens a new rgl device, and sets the parameters as requested.
I want the new rgl device (window) to be sized and positioned in a specific
place on my screen. So, I try to set the "windowRect" parameter as fo
Check prop.test().
At 12:41 PM 9/17/2008, Bunny, lautloscrew.com wrote:
Hi all,
though i know this is a simple question, i really hope someone could.
I am just trying to compare percents respectively the difference in
percents by a simple test. t.test usually compares means and is
working like
Hi,
I have a dataframe with 14319rows and 9 colums. for some rows there are null
values.I want a dataframe without these null values.I wanna select only
those that have values !=NA.
kindly let me know how to do that.
Ramya
--
View this message in context:
http://www.nabble.com/selecting-dat
On Sep 17, 2008, at 1:22 PM, Rajasekaramya wrote:
Hi,
I have a dataframe with 14319rows and 9 colums. for some rows there
are null
values.I want a dataframe without these null values.I wanna select
only
those that have values !=NA.
kindly let me know how to do that.
?is.na
--
David
Dear all,
I have a column within a dataframe of values which range between 1 and 2. I
want to display graphically the distribution of these values (i.e. are they
clustered towards either exteme? Or spread evenly?). What is a good way of
doing this in R?
I've tried a few things, including usin
I have been using the nb2WB() package for the car.normal function in WinBUGS,
but it will not create the C[] matrix; it only creates adj[], num[], and
weights[]. I was planning on using this function to create the C[] matrix
(by using the num matix) required for the car.proper, but I got slipped u
Could someone please tell me how to stop the package/function name from
being included before the command prompt? This started happening today
after I made some changes to my Rprofile.site file and I don't know why.
For example: if I enter example(AIC), instead of just getting the
regular '>' befo
This is to announce that we plan to release R version 2.8.0 on Monday,
October 20, 2008.
Release procedures start Friday, September 19. The detailed schedule can
be found on http://developer.r-project.org
The source tarballs will be made available daily (barring build
troubles) and the tarballs c
On 9/17/2008 1:13 PM, Koen Stegen wrote:
Dear all,
The documentation on the function open3d in the rgl package, reads:
open3d opens a new rgl device, and sets the parameters as requested.
I want the new rgl device (window) to be sized and positioned in a specific
place on my screen. So, I
Mike,
There are several options to pursue, depending on your answer
to the question:
Why are you staying within the Rgui if you are already an
Emacs ESS user?
1. If you don't need any other feature of Rgui on Windows,
then don't use it. From emacs, enter
M-x R
and you will immediately
Farr, David wrote:
> Could someone please tell me how to stop the package/function name from
> being included before the command prompt? This started happening today
> after I made some changes to my Rprofile.site file and I don't know why.
> For example: if I enter example(AIC), instead of just g
On 9/15/08, Waichler, Scott R <[EMAIL PROTECTED]> wrote:
> When I use panel.contourplot() with filled color regions, the coloring
> follows the stair-step edge of the underlying grid instead the smooth
> contour lines themselves. How can I get the latter behavior? I would
> guess there is a mu
On 9/17/2008 1:42 PM, Farr, David wrote:
Could someone please tell me how to stop the package/function name from
being included before the command prompt? This started happening today
after I made some changes to my Rprofile.site file and I don't know why.
For example: if I enter example(AIC), i
Hi:
I've searched the archives and the Internet for hours but have yet to find a
way to do stacked area plots (like the kind in Excel) in R. I think that
polygon may work but it would require a bit of manipulation of my data. I was
hoping for an easier way.
An example of what I'm trying to do
On Wed, Sep 17, 2008 at 1:01 PM, Josip Dasovic <[EMAIL PROTECTED]> wrote:
> Hi:
>
> I've searched the archives and the Internet for hours but have yet to find a
> way to do stacked area plots (like the kind in Excel) in R. I think that
> polygon may work but it would require a bit of manipulatio
Steve Murray hotmail.com> writes:
>
>
> Dear all,
>
> I have a column within a dataframe of values which range between 1 and 2. I
want to display graphically the
> distribution of these values (i.e. are they clustered towards either exteme?
Or spread evenly?). What is
> a good way of doing thi
Hi Jathine,
>> I hope this can explain the problem a bit more clearly.
>> Why PCA gives different results on the two different platforms?
What is amazing, Jathine, is how nearly exactly identical the two sets of
results are, not that they begin to differ at the 16th decimal place. To
assuage yo
Dear Phil and Jorge,
Many thanks for your quick replies. I found that:> hist(urban.long[,3]) worked
and displayed the data as I hoped. This reveals that the data are strongly
distributed towards the value '1', and the other bars on the histogram are
difficult to distinguish from each other a
I have a problem in where i generate m independent draws from a
binomial distribution,
say
draw1 = rbinom( m , size.a, prob.a )
then I need to use each draw to generate a beta distribution. So,
like using a beta prior, binomial likelihood, and obtain beta
posterior, m many times. I have
> I have a column within a dataframe of values which range between 1 and
> 2. I want to display graphically the distribution of these values
> (i.e. are they clustered towards either exteme? Or spread evenly?).
> What is a good way of doing this in R?
>
> I've tried a few things, including using
Hi:
I am trying to run your SCHOOLS example as following and I get an error message
stated that : Error: could not find function "bugs"
Would you please help me since I would like to run some WINBUG programs from
your R by reading external data and I am new to R?
Do I missing some function? Sh
Hello R users,
I want to readout the row and column postion from a certain matrix value
into a csv file.
I have only found this syntax
"which(a == b, arr.ind = T)"
so I get
a = matrix
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,]11222112
[2,]123
Thank you for your reply.
Here are some more info, I hope this can explain the problem a bit more
clearly.
Why PCA gives different results on the two different platforms?
freqtest.txt file line text :
M1 M2 M3 M4 M5 M6 M7 M8
-1 -1 -1 -1 -1 -1 -1 -1
0 0 0 0 -1 -1 1 1
-1 -1 -1 -1 -1 -1 -1 -1
0 0
I wrote a function that reads data from a file and manipulates it to build
and return a matrix with a million rows and three columns. I also need to
return some metadata: these are about a half-dozen single numbers or short
character strings, plus one numeric vector of length 5,000. I've been
in
Thanks a lot for your suggestion but...
Can I do that using the/a land-sea information from the 'worldHires'
database (from mapdata)?
Best,
VÃctor.
Greg Snow wrote:
After you have interpolated your data, replace the values in the areas of the o
cean with NA, then filled.conto
Hi Jathine,
And then to see things more clearly still, you can do something like this on
your test results:
format(formatC(p1$var$coord, digits=15, format="f"), justify="right")
and
format(formatC(p1$var$coord, digits=16, format="f"), justify="right")
Though I do hope that the second command
Hi Chris82,
Try this:
res=c(which(a == max(a), arr.ind = T),max(a))
names(res)=c('row','col','value')
write.csv(res,"C://yourfile.csv",row.names=TRUE)
HTH,
Jorge
On Wed, Sep 17, 2008 at 1:47 PM, Chris82 <[EMAIL PROTECTED]> wrote:
>
> Hello R users,
>
> I want to readout the row and column
Gao, Aijun i3statprobe.com> writes:
>
>
> Hi:
>
> I am trying to run your SCHOOLS example as following and I get an error
message stated that : Error: could not
> find function "bugs"
> Would you please help me since I would like to run some WINBUG programs from
your R by reading external data
On 9/17/2008 2:29 PM, Victor Homar wrote:
Thanks a lot for your suggestion but...
Can I do that using the/a land-sea information from the 'worldHires'
database (from mapdata)?
Yes, you use the map.where function to identify the location of the
points. If NA, it's out in the ocean s
Hi,
you might try this:
set.seed(100)
m <- 10
size.a<- 10
prob.a<- 0.3
prior.constant = 0
draw1 = rbinom( m , size.a, prob.a )
beta.draws <- function(draw, size.a, prior.constant, n) {
rbeta(n, prior.constant + draw, prior.constant + size.a - draw)
}
bdraws <- sapply(d
Both shape parameters of rbeta can be vectors; for
x <- rbeta(n, shape1, shape2)
x[i] ~ Beta(shape1[i], shape2[i])
so
bbsim <- function(m=1000, num.post.draws=1e4, size.a=100, prob.a=.27,
prior.count=1) {
data.count <- rbinom(m, size.a, prob.a)
shape1 <- rep(prior.count + data.count, e
I submitted this to rmetrics help list as well...
I've been using R-2.4 with Rmetrics successfuly for a year or two.
I recently moved to Ubuntu Linux 8.04, loaded R via apt-get install R-base
etc. etc...
then in the R interface i ran install.packages("fMultivar"), chose the CRAN
mirror and it lo
On Wed, 17 Sep 2008, Peter Dalgaard wrote:
fisher.test(, simulate.p.value=TRUE) might be more direct. Also
works for chisq.test().
And, contrary to popular belief, fisher.test() does work for larger
than2x2 tables, although you may run into space/time limitations.
The 'exact' test comput
> I therefore wish to examine all values of urban.long[,3] which are
> greater than 1. I have tried the following, but receive error
> messages each time:
> > hist(urban.long[,3]>1)
> Error in hist.default(urban.long[, 3]> 1) : 'x' must be numeric
> > hist(urban.long[urban.long[,3]>1])
> Error i
I think this is a very useful function that I imagine has wide appeal -
thanks. Using the code below produces the plot OK but when I try and
copy/save it (as a metafile) I receive the following error and an hourglass:
Error: invalid graphics state
[using XP, 2.72, lattice 0.17.13]
Regards,
Dav
Dear R-users,
I have some very simple data where 1's represent events and zeroes non-events,
e.g.
temp <- rbind(c(0,1,0,0,1,1,1,0), c(0,0,0,1,0,0,0,0))
For each row in the matrix, I would like to replace a singelton event by a 0.
That is, any 1 surrounded by zeroes (010) should be replaced by a
On Sep 17, 2008, at 3:36 PM, Jorge Ivan Velez wrote:
Hi Chris82,
Try this:
res=c(which(a == max(a), arr.ind = T),max(a))
names(res)=c('row','col','value')
write.csv(res,"C://yourfile.csv",row.names=TRUE)
I wasn't sure that max(x) or whatever target was specified would be
unique, and in th
Try this:
do.call(rbind, strsplit(gsub("010", "000", apply(temp, 1, paste, collapse =
"")), NULL))
On Wed, Sep 17, 2008 at 4:50 PM, Stacey Burrows <[EMAIL PROTECTED]>wrote:
> Dear R-users,
>
> I have some very simple data where 1's represent events and zeroes
> non-events, e.g.
> temp <- rbind(c
On Wed, Sep 17, 2008 at 1:12 PM, David Carslaw
<[EMAIL PROTECTED]> wrote:
>
> I think this is a very useful function that I imagine has wide appeal -
> thanks. Using the code below produces the plot OK but when I try and
> copy/save it (as a metafile) I receive the following error and an hourglass
I'm trying to use the following loop to open a window multiple times
to select files, but only the last window shows up. What am I missing?
library(tcltk)
nWin <- 6
fn <- vector('list', nWin)
for (ii in nWin) {
fn[[ii]] <- tclvalue( tkgetOpenFile( filetypes =
"{{Files} {.1D}} {{All files
Javier Acuña gmail.com> writes:
>
> Hi, I'm a new user of R. My background is Electrical Engineering, so
> please bear with me if this is a silly question.
>
> I'm trying to assess whether the results of an experiment satisfy the
> hypothesis of homoscedasticity (my ultimate goal is to use ANOV
On Wed, Sep 17, 2008 at 1:25 PM, Deepayan Sarkar
<[EMAIL PROTECTED]> wrote:
> On Wed, Sep 17, 2008 at 1:12 PM, David Carslaw
> <[EMAIL PROTECTED]> wrote:
>>
>> I think this is a very useful function that I imagine has wide appeal -
>> thanks. Using the code below produces the plot OK but when I tr
Bo,
Try using RSiteSearch with the strings 'huber-white', 'sandwich' or
even 'clustered standard errors'.
You may also want to consider a mixed models approach -- see:
http://www.stat.columbia.edu/~cook/movabletype/archives/2007/11/clustered_stand.html
HTH,
Kingsford Jones
On Tue, Sep 16, 2
Hi, Sorry for bothering your guys. I am still working on boxplot.
How can I extend the whisker to the 5% and the 95% quantiles? What I should
do if I only want to show the most extreme outlier, like 0.01% and 99.99%
percentiles?
I saw something on boxplot.stat, but I even donot know how to conn
Dear Help List,
Thanks in advance for reading...I hope my questions are not too ignorant.
I have an experiment looking at evolution of wing size [centroid] in
fruitflies and the effect of 6 different experimental treatments
[treatment]. I have five replicate populations [replic] in each
treat
1 - 100 of 146 matches
Mail list logo
