Hello R users,
I have a montly time series over a several year period.
It's easy to compute a monthly climatology (12 values), Now, I would like to
calculate the corresponding standard deviation, ie the 12 values
calculated from the january values, february values, etc.
What's the best way for suc
Paul
you might find the view() function in the 'elliptic' package useful.
This function implements various methods to visualize functions
over the complex plane.
HTH
rksh
Paul Smith wrote:
Dear All,
The function curve() draws the graph of functions from R to R. Is
there some homologous f
Hello,
Thanks for the response.
However, I have this problem with view.
> View(fit1)
Error: could not find function "View"
I have R 2.4.1 and even tried loading 'utils'.
As for saving residuals in a spreadsheet I used something like this:
res.fit1<-fit1$residuals
res.fit2<-fit2$residuals
Dear Peter,
thank you so much for the information and quick help, this was my problem
to get the "configure" run.
please see the results:
R is now configured for hppa2.0w-hp-hpux11.11
Source directory: .
Installation directory:/opt/R
C compiler:gcc -std=gnu99
DumpsterBear wrote:
I am getting "Error: cannot allocate vector of size 197 MB".
I know that similar problems were discussed a lot already, but I
didn't find any satisfactory answers so far!
Details:
*** I have XP (32bit) with 4GB ram. At the time when the problem
appeared I had 1.5GB of avail
[EMAIL PROTECTED] wrote:
> Make: Don't know how to make #. Stop.
> *** Error exit code 1
>
Sounds like you are using a "make" which doesn' t understand comments
(??!). You might need to install GNU make and possibly set MAKE=gmake or
some such to make sure that you are using the right one. At l
Dear all,
I am currently learning to run logistic regression models with R. Would
someone tell me what this exactly means:
Estimated scale (compare to 1 ) 1.746724
If the value is higher or lower than 1, what should I do?
The complete results of the model were as follows:
Generalized linear
"7" is an integer, but it's also a real.
In R '?is' and '?is.integer' are clear that you're testing the class(es) of
objects, not their values.
I can't comment on the relationship with "S Programming"
hth
Keith J
<[EMAIL PROTECTED]> wrote in message
news:[EMAIL PROTECTED]
> This is really bo
Dear UseRs,
please find below a simple example which exemplifies my question.
In words, I'd like to write a function which prints formula and outcomes only
when results are not assigned or recalling the already assigned object, as e.g.
glm.
I'm sure a solution is rather simple and straightforwa
Em Ter, 2008-09-23 às 21:42 -0400, DumpsterBear escreveu:
> I am getting "Error: cannot allocate vector of size 197 MB".
> I know that similar problems were discussed a lot already, but I
> didn't find any satisfactory answers so far!
>
> Details:
> *** I have XP (32bit) with 4GB ram. At the time
Ajay ohri wrote:
Hi List,
Graphical output to PDF's ,RTF ,CSV is known through R.
Can it be modified for outputting to Google Docs (which is basically
uploaded files ,published to become html pages)
Is there any package on this ?
Hi Ajay,
Do you mean transforming R output into HTML? The R2H
Do I mean R output to Google docs, like in Google spreadsheets which are
different from a spreadsheet published as an html
or Presentations in Google Docs, which are much different tahn ordinary
html presentations.
See examples here of Google Docs-
https://spreadsheets.google.com/ccc?key=pS9vSxWu
Werner Wernersen wrote:
Hi,
I am trying to prepare a report with R2HTML using
HTMLStart(outdir="./html", filename="report", echo=T, HTMLframe=F)
then, for instance, I want to get the output of a loop:
for (i in 1:20) print(summary(rnorm(1000)))
but only the first of summaries really ends up
Hi,
I am working with a phylog tree and I would like to extract a subset of the
tree based on the species names (conserving the evolutionary distance and
relationships between the pairs of species I am interested in). I see there
is an option to select a subset of the tree using node names
(phylog
Quoting Peter Dalgaard <[EMAIL PROTECTED]>:
[EMAIL PROTECTED] wrote:
Hi,
I have been using 'read.table' regularly to read tab-delimited text
files with data. No problem, until now.
Now I have a file that appeared to have read fine, and the data inside
looks correct (structure etc), except I
The result was 11, 24001 times, as I expected originally. hmmm...
JOse
Quoting Gabor Grothendieck <[EMAIL PROTECTED]>:
Try looking at the result of count.fields to diagnose it.
On Tue, Sep 23, 2008 at 5:19 AM, <[EMAIL PROTECTED]> wrote:
Hi,
I have been using 'read.table' regularly to r
Dear list,
with the code below I produce the right graph, but the colours of the
legend are different from the colours of the graph. The colours of the
graph are the desired colours.
Thanks for any help.
Patrick
library(lattice)
pal1 <- rgb(196, 255, 255, max = 255)
pal2 <- rgb( 0, 35,
Shengqiao Li wrote:
On Tue, 23 Sep 2008, Duncan Murdoch wrote:
On 23/09/2008 4:00 PM, Shengqiao Li wrote:
How to use sub, gsub, etc. to replace "\" in a string to "/"?
For example,convert "C:\foo\bar" to "C:/foo/bar".
If those are R strings, there are no backslashes in the firs
Hi,
you can use par.settings to get a consistent color scheme,
pal1 <- rgb(196, 255, 255, max = 255)
pal2 <- rgb( 0, 35, 196, max = 255)
df <- data.frame( Gruppe = c("A", "B", "A", "B"),
Kat = c("x1", "x1", "w1", "w1"),
value= c(1,2, 4, 5))
barchart(value ~ Kat, group=
Is it possible to convert an hexbin object like
hexbin(rep(1:100,each=100),rep(1:100,100),xbins=10)
and convert it to a vectorial (shp) format?
Paulo
-
This electronic mail
Hello list,
I've been working on this problem for a while and I haven't been able
to come up with a solution.
I have a couple of functions that plot a bunch of data, then a single
point on top of it. What I want is to be able to change the plot of
the point without replotting all the data. Cons
Check out JohnsonFit in SuppDists.
On Tue, Sep 23, 2008 at 10:59 AM, zhijie zhang <[EMAIL PROTECTED]> wrote:
> Dear R users,
> I hope to explain the concepts of skewness and kurtosis by generating
> series of distributions with same skewness and different kurtosis or with
> same kurtosis and dif
> "KJ" == Keith Jewell <[EMAIL PROTECTED]>
> on Wed, 24 Sep 2008 09:46:08 +0100 writes:
KJ> "7" is an integer, but it's also a real.
KJ> In R '?is' and '?is.integer' are clear that you're testing the
class(es) of
KJ> objects, not their values.
KJ> I can't comment on
On Wed, 24 Sep 2008, Duncan Murdoch wrote:
Shengqiao Li wrote:
On Tue, 23 Sep 2008, Duncan Murdoch wrote:
On 23/09/2008 4:00 PM, Shengqiao Li wrote:
How to use sub, gsub, etc. to replace "\" in a string to "/"?
For example,convert "C:\foo\bar" to "C:/foo/bar".
If those are R strings, th
On Wed, Sep 24, 2008 at 4:05 AM, Ben Bolker <[EMAIL PROTECTED]> wrote:
>> The function curve() draws the graph of functions from R to R. Is
>> there some homologous function to curve() to draw functions from R^2
>> to R?
>
> There is a curve3d function in the emdbook package on CRAN.
Thanks, Ben
It probably depends on how you are prompting the user for input. Here
is what happens with 'readline' and you will see that the string after
input does have the backslashes escaped.
> x <- readline("Input File: ")
Input File: c:\dir1\dir2\file.name
> x
[1] "c:\\dir1\\dir2\\file.name"
> gsub("
Last week I asked about data ellipses with rgl:::ellipse3d() with lines
showing the principal axes.
(The goal is a visual demonstration of PCA as a rotation of variable
space to component space.)
I was trying, unsuccessfully, to use princomp() to generate the PCA axes
and plot them using
segmen
Shengqiao Li wrote:
On Wed, 24 Sep 2008, Duncan Murdoch wrote:
Shengqiao Li wrote:
On Tue, 23 Sep 2008, Duncan Murdoch wrote:
On 23/09/2008 4:00 PM, Shengqiao Li wrote:
How to use sub, gsub, etc. to replace "\" in a string to "/"?
For example,convert "C:\foo\bar"
Dear All,
I am not sure whether the following can be considered a bug:
> x <- seq(-1,1,length=20)
> y <- x
> z <- matrix(1,20,20)
> persp(x,y,z)
Error in persp.default(x, y, z) : invalid 'z' limits
It works with:
persp(x,y,z,zlim=c(0,1.5))
Should not persp set the z limits by default when the
Don't think it's in kernSmooth, but locfit has such facility. See
Loader's book for details.
Best,
Andy
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Chinthaka Kuruwita
> Sent: Tuesday, September 23, 2008 12:18 PM
> To: r-help@r-project.org
> S
Dear list,
I would like to apply my own colours to a stacked area plot, done with
qplot, but I have not succeeded...
What do I have so far (I am dealing with the development of cover of
specific groups of plants):
library(ggplot2)
library(RODBC)
channel <- odbcConnect("myusername", ca
On Wed, 24 Sep 2008, jim holtman wrote:
It probably depends on how you are prompting the user for input. Here
is what happens with 'readline' and you will see that the string after
input does have the backslashes escaped.
x <- readline("Input File: ")
Input File: c:\dir1\dir2\file.name
x
Hello list,
I want to compare time needed for a simulation between R and Ra (time
optimized loops).
Does anybody know, how I can count time needed for a simulation in R?
Bests,
Ralph
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/li
On Wed, Sep 24, 2008 at 8:39 AM, Albin Blaschka
<[EMAIL PROTECTED]> wrote:
>
> Dear list,
>
> I would like to apply my own colours to a stacked area plot, done with
> qplot, but I have not succeeded...
>
> What do I have so far (I am dealing with the development of cover of
> specific groups of pla
?system.time
On Wed, Sep 24, 2008 at 9:46 AM, Ralph Scherer
<[EMAIL PROTECTED]> wrote:
> Hello list,
>
> I want to compare time needed for a simulation between R and Ra (time
> optimized loops).
> Does anybody know, how I can count time needed for a simulation in R?
>
> Bests,
> Ralph
>
>
Dear Yogesh,
Send us a minimal reproducive code, with a sample of your data structure.
It will help us to help you :-)
Best wishes,
miltinho astronauta
brazil
On Wed, Sep 24, 2008 at 2:00 AM, Yogesh Tiwari <[EMAIL PROTECTED]>wrote:
> Hello R users,
>
> I have a montly time series over a severa
Michael Friendly wrote:
Last week I asked about data ellipses with rgl:::ellipse3d() with lines
showing the principal axes.
(The goal is a visual demonstration of PCA as a rotation of variable
space to component space.)
I was trying, unsuccessfully, to use princomp() to generate the PCA axes
an
Say I have the following data:
testDat <- data.frame(A = c(1,NA,3), B = c(NA, NA, 3))
> testDat
A B
1 1 NA
2 NA NA
3 3 3
rowsums() with na.rm=TRUE generates the following, which is not desired:
> rowSums(testDat[, c('A', 'B')], na.rm=T)
[1] 1 0 6
rowsums() with na.rm=F generates the fol
Your users will have to learn that if they are inputting quoted
strings into R, then, by convention, a backslash is used to 'escape'
certain character sequences and if you want a backslash, you have to
escape it ('\\'). You can also have your users upgrade to a system
that does not use backslashes
I have an MS Access database with one table and one column holding rep(1:10)
I use:
library(RODBC)
channel <- odbcConnect("test")
sqlQuery(channel, paste("SELECT col FROM tblTest"), believeNRows=FALSE)
and get:
49
50
51
52
53
54
55
56
57
12337
What? The above should, of course, read
Hi,
How can I generate lower case letters for my axis-titles?
Thanks,
Jörg
__
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 comm
>> I am getting "Error: cannot allocate vector of size 197 MB".
>> I know that similar problems were discussed a lot already, but I
>> didn't find any satisfactory answers so far!
>>
>> Details:
>> *** I have XP (32bit) with 4GB ram. At the time when the problem
>> appeared I had 1.5GB of available
R Help gmail.com> writes:
>
> Hello list,
>
> I've been working on this problem for a while and I haven't been able
> to come up with a solution.
>
> I have a couple of functions that plot a bunch of data, then a single
> point on top of it. What I want is to be able to change the plot of
> t
Duncan Murdoch wrote:
The normals component contains the surface normals. It is used to
help in rendering the surface, but isn't much use for your purposes.
Unfortunately, I'm not familiar enough with the internals of princomp
to tell you how to get the axes you want.
Duncan Murdoch
OK, t
hadley wickham schrieb:
On Wed, Sep 24, 2008 at 8:39 AM, Albin Blaschka
<[EMAIL PROTECTED]> wrote:
[snip]
As I have those categories (cover grasses, cover herbs...) I would like to
assign special colours to each (for example "cover grasses" => dark green,
"cover open soil" => brown ...), but I
check:
?tolower()
Best,
Dimitris
Jörg Groß wrote:
Hi,
How can I generate lower case letters for my axis-titles?
Thanks,
Jörg
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
htt
I guess this would be the fastest way would be:
rs <- rowSums( testDat, na.rm=T)
rs[ which( rowMeans(is.na(testDat)) == 1 ) ] <- NA
since both rowSums and rowMeans are internally coded in C.
Regards, Adai
Doran, Harold wrote:
Say I have the following data:
testDat <- data.frame(A = c(1,N
One way is to keep a copy of the original and then return to it when you
need it.
x <- rnorm(100,1,0.5)
y <- rnorm(100,1,0.5)
plot(x,y,pch=16)
original <- recordPlot()
for( i in 1:10 ){
points( x[i], y[i], pch=19, col="yellow", cex=3)
points( x[i], y[i], pch=16)
Sys.sleep(1)
On 9/24/2008 10:06 AM, Doran, Harold wrote:
> Say I have the following data:
>
> testDat <- data.frame(A = c(1,NA,3), B = c(NA, NA, 3))
>
>> testDat
>A B
> 1 1 NA
> 2 NA NA
> 3 3 3
>
> rowsums() with na.rm=TRUE generates the following, which is not desired:
>
>> rowSums(testDat[, c('A',
Oh, ok ...
perhaps it was the wrong description;
I want lowered cases...
In LaTeX I get them with x_1 for example.
Or x^2 for upper cases.
But how can I make this in R?
Am 24.09.2008 um 16:17 schrieb Dimitris Rizopoulos:
check:
?tolower()
Best,
Dimitris
Jörg Groß wrote:
Hi,
How can I
On Wed, Sep 24, 2008 at 9:59 AM, Shengqiao Li <[EMAIL PROTECTED]> wrote:
>
> On Wed, 24 Sep 2008, jim holtman wrote:
>
>> It probably depends on how you are prompting the user for input. Here
>> is what happens with 'readline' and you will see that the string after
>> input does have the backslash
so you mean subscripts and superscripts, then check:
?plotmath
Best,
Dimitris
Jörg Groß wrote:
Oh, ok ...
perhaps it was the wrong description;
I want lowered cases...
In LaTeX I get them with x_1 for example.
Or x^2 for upper cases.
But how can I make this in R?
Am 24.09.2008 um 16:17
try the following:
testDat <- data.frame(A = c(1,NA,3), B = c(NA, NA, 3))
ind <- rowSums(is.na(testDat)) == length(testDat)
out <- rowSums(testDat, na.rm = TRUE)
out[ind] <- NA
out
I hope it helps.
Best,
Dimitris
Doran, Harold wrote:
Say I have the following data:
testDat <- data.frame(A =
An example would help.
You generally control the titles using arguments like main, xlab, ylab,
sub in the plotting functions or afterwards using title() function. You
can get the upper/lower case using toupper()/tolower() functions. See
help(par), help(title), help(tolower). Here is an example
on 09/24/2008 09:06 AM Doran, Harold wrote:
> Say I have the following data:
>
> testDat <- data.frame(A = c(1,NA,3), B = c(NA, NA, 3))
>
>> testDat
>A B
> 1 1 NA
> 2 NA NA
> 3 3 3
>
> rowsums() with na.rm=TRUE generates the following, which is not desired:
>
>> rowSums(testDat[, c('A',
On 9/24/2008 10:38 AM, Marc Schwartz wrote:
> on 09/24/2008 09:06 AM Doran, Harold wrote:
>> Say I have the following data:
>>
>> testDat <- data.frame(A = c(1,NA,3), B = c(NA, NA, 3))
>>
>>> testDat
>>A B
>> 1 1 NA
>> 2 NA NA
>> 3 3 3
>>
>> rowsums() with na.rm=TRUE generates the following
Great thanks and sorry for my bad english!
Am 24.09.2008 um 16:32 schrieb Dimitris Rizopoulos:
so you mean subscripts and superscripts, then check:
?plotmath
Best,
Dimitris
Jörg Groß wrote:
Oh, ok ...
perhaps it was the wrong description;
I want lowered cases...
In LaTeX I get them with
> I agree that that some sort of facility would be convenient. Creating
> latex output is
> another situation where not having to escape backslashes would be convenient.
And regular expressions, of course.
Hadley
--
http://had.co.nz/
__
R-help@r-pro
Thanks Charles, ftable() works perfectly.
-Original Message-
From: Charles C. Berry [mailto:[EMAIL PROTECTED]
Sent: Tuesday, September 23, 2008 5:06 PM
To: Hutchinson,David [PYR]
Cc: r-help@r-project.org
Subject: Re: [R] Counting character occurrences in data frame
See
?ftabl
Is there an elegant way in R to change a number reported as a less-than
number in text format, "<1" for example, to the numeric equivalent 1? I have
been trying to use as.numeric, but have not come up with anything clever
yet.
Tom
--
View this message in context:
http://www.nabble.com/How-do-
On 24/09/2008 10:12 AM, Michael Friendly wrote:
Duncan Murdoch wrote:
The normals component contains the surface normals. It is used to
help in rendering the surface, but isn't much use for your purposes.
Unfortunately, I'm not familiar enough with the internals of princomp
to tell you how to
You could always try
F(x,y) = f(x) + 0*y
That is "zero out" the degenerate dimensions. Of course you'll be
plotting what is essentially a two dimensional object as if it were
three dimensional. The degeneracy in y means a 2-D curve will be
"extruded" along the Y dimension.
Robert Farley
Met
Hello everyone,
I have the following problem:
My analysis includes many predictor variables (>50) in the form of
raster maps (asc), but I am trying to avoid having to type all their
names over and over again in the analysis (e.g. for vectorisation, for
deletion of NA's, etc.)
So ideally I wou
On 24/09/2008 10:53 AM, hadley wickham wrote:
I agree that that some sort of facility would be convenient. Creating
latex output is
another situation where not having to escape backslashes would be convenient.
And regular expressions, of course.
There were proposals to do this last year (I t
Hi, R-users,
If I have a data frame like this:
>x<-data.frame(g=c("g1","g2","g1","g1","g2"),v=c(1,7,3,2,8))
g v
1 g1 1
2 g2 7
3 g1 3
4 g1 2
5 g2 8
It contains two groups, g1 and g2. Now for each group I want the max v:
> aggregate(x$v,list(g=x$g),max)
g x
1 g1 3
2 g2 8
Beautiful. But wh
On 24/09/2008 11:11 AM, Duncan Murdoch wrote:
On 24/09/2008 10:12 AM, Michael Friendly wrote:
Duncan Murdoch wrote:
The normals component contains the surface normals. It is used to
help in rendering the surface, but isn't much use for your purposes.
Unfortunately, I'm not familiar enough with
Thanks for directing my attention to the survey package and svyolr().
Best,
Greg
.
Gregory Wawro [EMAIL PROTECTED]
Associate Professor phone: 212-854-8540
Dept. of Political Science fax:212-222-0598
741 Inte
Hey list!
It looks simple, though it's not possible for me to plot the following
properly:
(some made-up data)
GrSeClone1 Clone2 Clone3 Clone4 Clone...
G1999 2 3 6 5
G2000 2 5 7 4
G2001 5 3 7 3
G2002 4 5 8 3
Rather than aggregate, use order and duplicated as in this post:
https://stat.ethz.ch/pipermail/r-help/2008-September/173139.html
On Wed, Sep 24, 2008 at 11:21 AM, zhihuali <[EMAIL PROTECTED]> wrote:
>
> Hi, R-users,
>
> If I have a data frame like this:
>>x<-data.frame(g=c("g1","g2","g1","g1","g
Dear readers,
I have a basic question about how to use lme for my design. I haven't been able
to find an example in r-help that made it clear to me how to tackle this
problem and unfortunately I also cannot get hold of Pinheiro & Bates 2000. I
hope someone can help.
Data for the response va
If I have a set of data comprising a list of numbers of eggs on mosquito guts
that range from 1 to 157. How can I get R to draw a barchart of the
distribution of the data (i.e. x axis= number of eggs on a gut, y axis=number
of mosquitoes found with that number of eggs)? Data is listed below.
Sorry for re-posting the question, I did not get any reply.
Kindly reply please if any one can.
###
Hello R users,
I have a montly time series over a several year period.
It's easy to compute a monthly climatology (12 values), Now, I would like to
calculate the corresponding standard deviation, i
Dear all,
My question concerns using repetitions and simulations (loops?) in R. I am
very new R user, so any help that can be offered would be greatly
appreciated!
I am using fitdistr() to determine the distribution of empirical univariate
datasets, and ks.test to assess the goodness of fit.
I have a set of data that comprises genome numbers in single eggs from three
different parasite clones - 3D7, HB3, and MIX. I can draw a boxplot of the
genome numbers for each clonefed but how do I carry out a t test or ANOVA to
compare if the means are signifcantly different? (Data is listed b
Dear Tom,
I don't know whether it qualifies as elegant or clever, but you should
be able to use sub(), as in
> as.numeric(sub("<", "", c("1", "< -2", "<1")))
[1] 1 -2 1
I hope this helps,
John
On Wed, 24 Sep 2008 08:08:21 -0700 (PDT)
Tom La Bone <[EMAIL PROTECTED]> wrote:
>
> Is there an e
Thank you for all of you. Intuitively, 7 is an integer for people who
live in this planet. It is just very difficult for me to believe that
R does not think 7 is an integer but 7L is.
is.integer(7) # R 2.7.2
[1] FALSE
Thus, based on Martin's comments, I try it again on the S-PLUS 8.0 and
i
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Yogesh Tiwari
> Sent: Wednesday, September 24, 2008 5:43 AM
> To: [EMAIL PROTECTED]
> Subject: [R] climatological standard deviation- (question re-posted)
>
> Sorry for re-posting the question, I did n
On Wed, Sep 24, 2008 at 8:43 AM, Yogesh Tiwari
<[EMAIL PROTECTED]> wrote:
> Sorry for re-posting the question, I did not get any reply.
> Kindly reply please if any one can.
Most likely that's because you didn't clearly explain the problem. Please read
the posting guide (see link at bottom of eac
Since x$GrSe isn't a number, what do you expect R to do?
On Wed, Sep 24, 2008 at 11:28 AM, Schreiber, Stefan
<[EMAIL PROTECTED]> wrote:
> Hey list!
>
>
> It looks simple, though it's not possible for me to plot the following
> properly:
>
> (some made-up data)
>
> GrSeClone1 Clone2 Clone3 C
package zoo rollapply if I had to guess at what you want to do, but
heed the advice given above.
On Wed, Sep 24, 2008 at 12:11 PM, Sarah Goslee <[EMAIL PROTECTED]> wrote:
> On Wed, Sep 24, 2008 at 8:43 AM, Yogesh Tiwari
> <[EMAIL PROTECTED]> wrote:
>> Sorry for re-posting the question, I did not g
look at
histogram and densityplot in the lattice package
On Wed, Sep 24, 2008 at 9:26 AM, Georgina Sarah Humphreys
<[EMAIL PROTECTED]> wrote:
> If I have a set of data comprising a list of numbers of eggs on mosquito guts
> that range from 1 to 157. How can I get R to draw a barchart of the
> d
Hi all, I
am trying to run a linear mixed effect models in lmer() from the lme4
package using the weights option.
I am using the
R version 2.7.2 (2008-08-25) and lmer version in lme4_0.999375-26, which I
think it is the latest version!
I am getting and error message when I add the
optio
Thanks Duncan (& others)
Here is a function that does what I want in this case, and tries to do
it to work generally
with ellipse3d. (Note that I reverse the order of centre and scale
'cause I was bitten
by trying ellipse3d.axes(cov, mu))
# draw axes in the data ellipse computed by ellipse3d
Have you tried is.integer(7.0) in S-Plus? (I have)
Do you think 7.0 is integer?
IMHO in R there is nothing to be fixed (in this regard) except your
understanding.
This is a computer language, not English; intuition isn't reliable, so we
have help pages.
is.integer(x) is not intended to indicate
I'm new to R and to validation of logistic regression models via bootstrapping.
I see that there are various approaches to bootstrapping, and I am wondering
(1) which approach would be best for internal validation of predictive
discrimination, (2) are there any good resources that could be reco
Van Dyke, Alison wrote:
I'm new to R and to validation of logistic regression models via bootstrapping.
I see that there are various approaches to bootstrapping, and I am wondering
(1) which approach would be best for internal validation of predictive
discrimination, (2) are there any good re
On 24/09/2008 12:32 PM, Michael Friendly wrote:
Thanks Duncan (& others)
Here is a function that does what I want in this case, and tries to do
it to work generally
with ellipse3d. (Note that I reverse the order of centre and scale
'cause I was bitten
by trying ellipse3d.axes(cov, mu))
# dr
On Wed, Sep 24, 2008 at 4:14 PM, Farley, Robert <[EMAIL PROTECTED]> wrote:
> You could always try
> F(x,y) = f(x) + 0*y
>
> That is "zero out" the degenerate dimensions. Of course you'll be
> plotting what is essentially a two dimensional object as if it were
> three dimensional. The degeneracy
IMHO Keith Jewell has a couple of "jewels" here (sorry Keith -- couldn't
resist): nice possible fortunes package candidates:
"IMHO in R there is nothing to be fixed (in this regard) except your
understanding."
" This is a computer language, not English; intuition isn't reliable, so we
have hel
Hi Keith,
No doubt, 7.0 is integer in math. But if people can write 7 why people
need to write 7.0 (I do not see any reason to do this). My point is
that R maybe can do something like S-plus. No point to argue. don't
you think so?
Thanks
Chunhao
Quoting Keith Jewell <[EMAIL PROTECTED]>:
I don't remember which version View was added in, but it was more recent than
2.4.1. I would suggest upgrading (2.4.1 is ancient in R terms).
You can get something similar with:
> invisible(edit(mydata))
But you will have to close the window to continue with R (View lets you look at
the data
I would like to use the data below where the plots are close to what I
want. Instead of color I would like to use different symbols, and
have the symbols in the legend match the graphs. I am also going to
add a regression line to these I know about the type="r" (which is
fine for these particular
Dear R People:
I finally (Yay!) got R installed in a classroom!
Anyhow, I have a respectful request, please: could anyone recommend
some nice undergrad projects in R, please?
This is in a statistical computation class; first time being run.
Thanks,
Erin
--
Erin Hodgess
Associate Professor
D
Hello,
I have been using sink to create text file outputs.
>sink("summ_model1and2.txt")
> summary(model1)
> summary(model2)
> sink()
Q1: Is there a way I could add a line of the text above the summary to act
like a title?
Also, I have been using the following to save plots from the lm function:
And here is one way of using JohnsonFit to display how a distribution changes
with changes in the moments:
library(TeachingDemos)
library(SuppDists)
tmpfun <- function(m1=0,m2=1,m3=-0.2,m4=3.7, xmin=-3, xmax=3){
parms <- JohnsonFit( c(m1,m2,m3,m4), 'use')
plot(function(x) dJohnso
Thank you Peter. That is incredibly helpfyul, and much much smaller!
Best,
Collin.
__
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
The command line app is fine, although there are some potential security
issues or annoyances with a smooth work flow.
But it would be more flexible and powerful to do it directly in R.
And we don't need to know or program in C. The RCurl (https, getForm()
and postForm()) and XML packages provi
Wow.
Dont you get to name it on me or something (like the RPackageGoogleAjayDoc
?) Just kidding
This is fabulous and the reason this community is so innovative
Now if only someone can create a package for me to use R from Amazon EC2 and
be free of the damn RAM problem forever ...
Am 23.09.2008 um 23:57 schrieb Peter Dalgaard:
For this kind of problem I'd go directly for the binomial
distribution. If the actual probability is 0, this is essentially
deterministic and you can look at
> binom.test(0,99,p=.03, alt="less")
This means that you don't sample from the p=.03
Remember that computers are not as smart as you. Some smart people have
written instructions for the computer on what to do in certain cases, but they
can't anticipate everything, so when you tell the computer to do something that
was different from what is anticipated, it either gives an error
1 - 100 of 157 matches
Mail list logo