On Jan 2, 2012, at 1:31 AM, 王琦 wrote:
hello:
I am trying to use R to draw a 3D picture, then color the picture
according to the value of z , how could I do this job?
this is my exalple
x<-y<-seq(-50,50,2)
m<-function(x,y) x^3+y^3-x^2*y-x*y^2
z<-outer(x,y,m)
persp(x,y,z,theta=-60,phi=30)
I
hello:
I am trying to use R to draw a 3D picture, then color the picture according to
the value of z , how could I do this job?
this is my exalple
x<-y<-seq(-50,50,2)
m<-function(x,y) x^3+y^3-x^2*y-x*y^2
z<-outer(x,y,m)
persp(x,y,z,theta=-60,phi=30)
I want to draw the Zmax in red and Zmin in b
On Sun, 1 Jan 2012 09:50:24 -0600
Michael wrote:
> Happy New Year all!
>
> I am looking for a good solution for keeping record of my experiments
> - could you please help me?
>
If you are working from scripts, which is a very good way to
standardize procedures as a work flow and analysis devel
Dear R users,
I fitted a GEE model using the function 'geese' (or 'geeglm') with user
defined
correlation matrix. I want to get the var-cov matrix of the regression
coefficients. But the output provides only limited information.
I would be very much thankful if you could kindly let me know ho
Hi folks,
Any ideas on this? This does sound like a fairly common situation - reading
in large data file into R?
Thanks.
Andy
--
View this message in context:
http://r.789695.n4.nabble.com/Reading-large-sparse-arff-files-into-R-tp4249409p4252393.html
Sent from the R help mailing list archive a
Dear Simon,
I have the same problem. I understand te(a), te(b) are nested in te(a,b)
according to your paper on tensor product. I have no enough data to perform
te(a,b,d) and only care the interactions a*b and a*d, so I did
y=te(a,b)+te(a,d). The resutl is good. I am wondering if this is the
co
I'm trying to read CODA/mcmc files (see the coda package), as
generated by jags/WinBUGS/OpenBUGS, into a big.matrix.I can't load
the whole mcmc object produced by read.coda() into memory since I'm
using a laptop for this analysis (currently I'm unfunded).
Right now I'm doing it by creat
Dear Michael,
On Sun, Jan 1, 2012 at 10:50 AM, Michael wrote:
> Happy New Year all!
[snip]
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
Nobody's yet mentioned that the latest version of Emac
Thanks. That did it!
And I get it now--in your original example, aes(x = x, y = Freq), x
refers to the column name in as.data.frame(table(x)), not the x
vector(?) you created.
Aren
On Sun, Jan 1, 2012 at 4:44 PM, Joshua Wiley wrote:
> Sorry, that was probably a really confusing example...too ma
Josh, you've solved the problem, fantastic.
Thanks for as.formula() David, that will be of great use in my work.
Next time I'll provide better examples. Thanks for your help all.
-
Isaac
Research Assistant
Quantitative Finance Faculty, UTS
--
View this message in context:
http://r.7
On Jan 1, 2012, at 5:50 PM, Rich Shepard wrote:
A data frame has one factor, one date, and one numeric column. When
I plot
these using the default pch of the open circle (first attachment),
xyplot(TDS ~ sampdate | she.s, data = sheep.cast, main = 'TDS in Sheep
Creek', ylab = 'Concentration
On Sun, 1 Jan 2012, Rich Shepard wrote:
Where do I start to look for the reason for this behavior?
Perhaps the dates are not in sequence? So I need to use zoo to order the
dates sequentially?
Rich
__
R-help@r-project.org mailing list
https://sta
A data frame has one factor, one date, and one numeric column. When I plot
these using the default pch of the open circle (first attachment),
xyplot(TDS ~ sampdate | she.s, data = sheep.cast, main = 'TDS in Sheep
Creek', ylab = 'Concentration (mg/L)', xlab = 'Time')
I see the higher concentrat
Sorry, that was probably a really confusing example...too many xs
floating around.
set.seed(10)
rawdata <- sample(0:23, 1, TRUE, prob = sin(0:23)+1)
## do theis step first for your data
tableddata <- as.data.frame(table(rawdata))
## use these names in ggplot
colnames(tableddata)
require(ggpl
Hi Ana,
most probably this is one of the more ugly solutions:
> d=dim(A)
> d
[1] 2 4
> cbind(rep(1:d[1], each=d[2]), rep(1:d[2], d[1]))
[,1] [,2]
[1,]11
[2,]12
[3,]13
[4,]14
[5,]21
[6,]22
[7,]23
[8,]24
Thanks,
wr
* Ana [2012-
Hi Randall,
This will do it. There may be more elegant ways. Formula methods are
just handy front ends (e.g., stats:::aggregate.formula), that end up
dispatching to other methods usually. It is easy to pass a character
vector to extractor functions like `[`() so with a bit more typing,
you can
This is helpful, although I can't seem to adapt it to my own data.
If I run your sample as is, I do get the nice graphs.
However, this doesn't work:
(Assume you already have a data frame "dallas" with 2057980 rows. It
has column "offense_hour", and each row has a value between 0 and 23,
inclusive
Thanks Michael.
That did the trick. Despite googling most of the day yesterday, I didn't
quite have the right search string to find that one. Almost feels like the
answer was hiding in plain sight, now that you point me to it.
I added some code to save the xyplots to a variable and then print it
Here's one way:
> data.frame(rowID=as.vector(row(A)), colID=as.vector(col(A)), A=as.vector(A))
rowID colID A
1 1 1 1
2 2 1 2
3 1 2 3
4 2 2 4
5 1 3 5
6 2 3 6
7 1 4 7
8 2 4 8
You can sort that as desired.
Sarah
On Sun, Jan 1, 2012 a
How can I extract a list of the positions in the matrix?
> A=matrix(1:8, nrow=2,ncol=4)
> A
[,1] [,2] [,3] [,4]
[1,]1357
[2,]2468
Something like this
pos.A
1 1
1 2
1 3
1 4
2 1
2 2
2 3
2 4
__
R-help@r-project.org m
Hello,
I have some code that currently works fine and I am endeavoring to
convert the major pieces of it into functions.
This involves taking "hard coded" names of variables that are used in
various places and figuring out how to
abstract them out into functions where the arguments (i.e. a list of
Look at the SWord package. It is available from rcom.univie.ac.at. It is
the integration of R and Word
similar to RExcel at the same site.
Rich
On Sun, Jan 1, 2012 at 3:31 PM, Joshua Wiley wrote:
> Hi Michael,
>
> I like Sweave and LaTeX, but I can appreciate the difficulty using it
> with col
On 12/28/2011 1:17 PM, David Winsemius wrote:
Along those lines, Thompson's "R accompaniment" to Agresti's text on
"Categorical Data Analysis" is simply wonderful:
https://home.comcast.net/~lthompson221/SplusDiscrete.PDF
It's not surprising that Agresti recommends it on the book's homepage:
htt
Hi Michael,
I like Sweave and LaTeX, but I can appreciate the difficulty using it
with collaborators. What about something similar using HTML?
Certainly integrates to any webpages nicely. There are two packages I
think do this nicely, one is the R2HTML package (on CRAN). Another
one that is not
Tal,
Plese look at this thread
https://stat.ethz.ch/pipermail/r-help/2011-October/292786.html
right justify right-axis tick values in lattice
and see if the function provided by David Winsemius will work for you.
I included it in the HH package as
HH:::panel.axis.right
Rich
On Sun, Jan 1, 2012
It should be noted that there is an Octave/Matlab clone for Android out there
(addi). It is not far fetched that someone given enough time, determination, a
rooted phone (LG phones are better for that, stay away from SE if you want a
rooted device) and Wilson :) could deploy R in the droid!
Chr
Hi Frederico,
This may not be what you're looking for, but we setup RStudio Server (
http://www.rstudio.org/ ) on a Linux box and access it from our tablets through
a web browser.
Best,
Tim
--
Timothy P. Jurka
Ph.D. Student
Department of Political Science
University of California, Davis
www.t
I don't know any language other than English, much less one written
right-to-left, but a quick search shows that the RGtk2 package includes
some support for such things.
My gut reaction is that right-to-left support would be highly device
dependent, and since devices often interact closely wit
okay thanks is there a way to programatically restart a new Console session?
On Sun, Jan 1, 2012 at 12:55 PM, Hadley Wickham wrote:
> >> My understanding is that one only clears the variables... not
> >> functions/packages, etc...
> >
> > Not exactly true. rm should remove any function you defin
If the phone is "rooted" one could hypothetically install from Debian
repositories where R is represented very comprehensively. I surmise this is
more trouble than it's worth in a major way. If stranded on a desert island and
getting R fully functional on Droid would save you, I'd advise looking
>> My understanding is that one only clears the variables... not
>> functions/packages, etc...
>
> Not exactly true. rm should remove any function you defined in your current
> session. You need to look at
>
> ?unloadNamespace
> ?detach
>
> ... in order to remove loaded packages.
And read the cav
On 12-01-01 9:05 AM, William Simpson wrote:
When using bmp() under Windows XP, I find that the saved image is a
shifted version of the correct image. Try this:
The image() function isn't designed to be able to do pixel-level
addressing, so it's not too surprising that some rounding error
some
On Jan 1, 2012, at 12:43 PM, Michael wrote:
Are there commands that can do more complete clean-up than
"rm(list=ls(all=T))"?
My understanding is that one only clears the variables... not
functions/packages, etc...
Not exactly true. rm should remove any function you defined in your
current
Are there commands that can do more complete clean-up than
"rm(list=ls(all=T))"?
My understanding is that one only clears the variables... not
functions/packages, etc...
Thanks a lot!
[[alternative HTML version deleted]]
__
R-help@r-project.or
Hi Isaac,
This is my best understanding of what you are after. A small worked
example would have been helpful (even if you had to manually write the
code, we would at least have a 'gold standard' for what we were trying
to automate to take lists of arbitrary size). This is not necessarily
an eff
I am not at my desk but you might search the CRAN web site for "reproducible
research".
Respectfully,
Frank Lawrence
-Original Message-
From: Michael
Sender: r-help-bounces@r-project.orgDate: Sun, 1 Jan 2012 09:50:24
To: r-help
Subject: [R] R report generator (for Word)?
Happy New Ye
Have you seen r2wd?
http://www.r-bloggers.com/exporting-r-output-to-ms-word-with-r2wd-an-example-session/
On Jan 1, 2012, at 10:50 AM, Michael wrote:
> Happy New Year all!
>
> I am looking for a good solution for keeping record of my experiments -
> could you please help me?
>
> My work is
On Jan 1, 2012, at 9:13 AM, iliketurtles wrote:
HI all,
I'm new to R.
Say I have a multi-layered list called newlist.
Which you have not provided in a form that lends itself well to
constructing worked examples.
str(newlist)
List of 2
$ :List of 5
..$ : num [1:8088] NA 4
Hi Aren,
I was busy thinking about how to make what you wanted, and I missed
that you were working with hours from a day. That being the case, you
may think about a circular graph. The attached plots show two
different ways of working with the same data.
Cheers,
Josh
set.seed(10)
x <- sample(
On Jan 1, 2012, at 11:07 AM, David Winsemius wrote:
On Jan 1, 2012, at 10:05 AM, Federico J. Villatoro wrote:
Hello,I am new to this list.
So new that you have not yet read the posting guide it appears.
I have ben wondering if there is an app to run
R on mobile devices? Specially android
On Jan 1, 2012, at 10:05 AM, Federico J. Villatoro wrote:
Hello,I am new to this list.
So new that you have not yet read the posting guide it appears.
I have ben wondering if there is an app to run
R on mobile devices? Specially android...Thanks in advance
Quick answer: Not yet (as such).
On Sun, Jan 1, 2012 at 5:29 AM, peter dalgaard wrote:
> Exactly. If what you want is a barplot, make a barplot; histograms are for
> continuous data. Just remember that you may need to set the levels
> explicitly in case of empty groups: barplot(table(factor(x,levels=0:23))).
> (This is irrel
HI all,
I'm new to R.
Say I have a multi-layered list called newlist.
> str(newlist)
List of 2
$ :List of 5
..$ : num [1:8088] NA 464 482 535 557 ...
..$ : num [1:8088, 1:2] NA 464 482 535 557 ...
..$ : num [1:8088, 1:3] NA 464 482 535 557 ...
..$ : num [1:8088, 1:4] NA 464
Thanks for passing that on, Uwe.
There are two separate things here:
- Uwe's problem is a bug: there is an experimental source-code-debugging
feature in 'debug' that doesn't work properly yet (pending a change in the
'parser' package) and for some stupid reason I'd set the default to "use
exp
Hello,I am new to this list. I have ben wondering if there is an app to run
R on mobile devices? Specially android...Thanks in advance
Federico
--
Federico
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat
When using bmp() under Windows XP, I find that the saved image is a
shifted version of the correct image. Try this:
n<-5
fn<-"01.bmp"
x<-matrix(runif(n*n),nrow=n)
image(x,col=gray(0:255/255),axes=F,frame.plot=F)
bmp(filename = fn,width = n, height = n, units = "px")
par(mar=c(0,0,0,0),pty="s")
ima
dear all,
i am looking for R software resources persons from India. Please reply.
Thanks & Regards
Sridhar
_
Sridhar Gutam PhD, ARS, PG Dip Patent Laws (NALSAR), IP & Biotechnology
(WIPO)
Senior Scientist (Plant Physiology) &
Happy New Year all!
I am looking for a good solution for keeping record of my experiments -
could you please help me?
My work is about analysing data... My current work-flow:
1. Everyday my bosses give me some small steps/tasks for analysing data -
which are parts of one bigger/whole project.
2.
On Jan 1, 2012, at 07:40 , Joshua Wiley wrote:
> If you just want a plot of the frequencies at each hour why not just call
> barplot on the output of table? Histograms create bins and count in those,
> which doesn't sound like what you're after.
>
Exactly. If what you want is a barplot, make
Hi Duncan,
Thank you for your reply.
I am also using Win 7.
And I would be surprised if this would be different in any OS.
I guess the answer is that there is no way for making text in a graph in R
be "right-to-left".
Thanks again, and happy new year,
Tal
Contact
Details:---
I'm guessing R FAQ 7.22: http://cran.r-project.org/doc/FAQ/R-FAQ.html
The subtlety is that in an interactive session print is automatically
called at the final evaluation of most everything, but you have to
prompt it in interactive use (and depending on details, in some
function calls)
Michael We
51 matches
Mail list logo