On Tue, 6 Nov 2007, Charilaos Skiadas wrote:
> Hello all,
>
> I ran into the following, to me unexpected, behavior. I have (for
> reasons that don't necessarily pertain to the question at hand, hence
> I won't go into them) the need/desire to use an empty string for the
> name of a vector entry.
Hello all,
I ran into the following, to me unexpected, behavior. I have (for
reasons that don't necessarily pertain to the question at hand, hence
I won't go into them) the need/desire to use an empty string for the
name of a vector entry. Perhaps I did not read ?"[" very carefully,
but it
Does this do any better, or perhaps offer other ideas? I could not
understand what the "Hct" and step()'s were doing from the description
of your desired result, so they were omitted from geq(). I also could
not understand why you were rounding off the integrals or why there wee
two integral te
Maura,
The workaround from https://bugzilla.novell.com/show_bug.cgi?id=330991
seems to fix this for most people (but not me): just add
export LIBXCB_ALLOW_SLOPPY_LOCK=1
to your ~/.profile file (or to /etc/profile to make it available for all users).
Ted.
Maura E Monville wrote on 11/03/200
Hi,
In the documentarion of sort() it reads (Arguments section):
"...: arguments to be passed to or from methods or (for the default
methods and objects without a class) to 'sim.int'."
should it be 'sort.int' instead of 'sim.int'?
Thank's
--
Henrique Dallazuanna
Curitiba-Paraná-Brasi
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Rohan7 wrote:
> Hello everybody,
>I'm using SSOAP and RCurl to call webservice. I got success in calling
> webservice with string as
> an input. Now i want to pass array as an input(can be single or
> multidimensional). When i do that i get
> err
A paper that may help you:
"Methods for Studying Coincidences", Persi Diaconis; Frederick
Mosteller. Journal of the American Statistical Association, vol 84, no.
408 (Dec., 1989), 853-861.
And remember that the birthday problem assumes independence, but if you
have 2 students that studied togeth
On Wed, 2007-11-07 at 09:28 +1300, David Scott wrote:
> In regard to validation of R, there are test directories with test
> code. See for example:
> https://svn.r-project.org/R/trunk/tests/
> and
> https://svn.r-project.org/R/trunk/src/library/stats/tests/
>
> David Scott
Those are a part of th
Thanks to all who helped out with this question. I have a version
that's working reasonably well for me right now, and will post code if I
manage a more generalizable version.
Best,
Ari
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailma
Alexandre,
Try rereading FAX 7.10, it explains why as.numeric() won't do it:
"In any case, do not call as.numeric() or their likes directly for the
task at hand as as.numeric() or unclass() give the internal codes"
I.e., the INTERNAL CODE of the factor is what as.numeric() is working
on rather t
On 7/11/2007, at 9:12 AM, Prof Brian Ripley wrote:
> 1) Did you merge the resources or restart X? You need to in order
> to get new resources to be recognized.
>
> xrdb -merge ~/.Xresources would be needed to merge resources into
> the current X session.
No I didn't. Didn't know or
OK, well, the way I suggested will allow you to modify the C code too --
but reading your question I realize you probably don't want to do this.
If you want to use the definitions of the C functions from stats, you can
skip 2. below and add another modification to the calls to .C in your
def. of bw
I tested
as.numeric(as.character(Ratio))
and it works perfectly!
I still don't get why as.numeric(Ratio) was not enough, but at least
now I know how to deal with it.
Thanks for the tip, and sorry for missing the R-FAQ issue 7.10.
Cheers,
Alexandre Santos
2007/11/6, John Kane <[EMAIL PROTECTE
Simple solution provide by Prof. Ripley is to modify the relevant code as
follows (he posted in the R-development group):
SDh <- function(x, h, n, d) .C(stats:::R_band_phi4_bin,
as.integer(n),
as.integer(length(x)), as.double(d), x, as.double(h),
u = double(1))$u
TDh
In regard to validation of R, there are test directories with test
code. See for example:
https://svn.r-project.org/R/trunk/tests/
and
https://svn.r-project.org/R/trunk/src/library/stats/tests/
David Scott
_
David Scott Depart
bw.SJ calls C code from the file
/usr/local/bin/R-2.6.0/src/library/stats/src/bandwidths.c
You have to make this code available.
You can do the following:
1. copy bandwidths.c and the definition of bw.SJ (the latter renamed) to a
directory
2. compile bandwidths.c into a shared library with the com
1) Did you merge the resources or restart X? You need to in order to get
new resources to be recognized.
xrdb -merge ~/.Xresources would be needed to merge resources into the
current X session.
I don't know if MacOS X11 is standard enough to support this, but it is a
very long-standing basic
Have a look at the R-FAQ issue 7.10. It's a standard
problem
For more information about your variable try
str(variable).
--- Alexandre Santos <[EMAIL PROTECTED]>
wrote:
> Dear R list,
>
> I would like to perform an ANOVA in a set of
> measurements, but I have
> problems formatting the data.
On 7/11/2007, at 6:24 AM, Dieter Menne wrote:
> you have little chance to get a C program to work
> without some experience.
A fortune???
cheers,
Rolf Turner
##
Attention:\ This
On 7/11/2007, at 6:39 AM, Prof Brian Ripley wrote:
> Better to use resources. But you can use (in .Rprofile)
>
> options(device=function() x11(xpos=-1))
>
> The resources version would be to put in ~/.Xresources something like
>
> R_x11*geometry: 700x700-0+0
Neither of these options app
I am a true R novice aonly using it for this function ;)
I am trying to use color2D.matplot to form a image of my data using the
following conditions
color2D.matplot(fi1, c(dr), c(dg), c(db), nslices=7, ylab='Species',
xlab="gene", show.legend=TRUE) where fi1 is my matrix.
I have a matrix with
On Tue, 06-Nov-2007 at 08:17AM -0200, Alberto Monteiro wrote:
|>
|> Prof Brian Ripley wrote:
|> >
|> >> I would like to know the created time and date of specific object.
|> >> Is there any function for it?
|> >
|> > There isn't even the concept. Most objects in R are a collection of
|> > SEXP
I'm looking at algorithms for determining coincidences. In educational
testing, it is interesting to look at cheating via the birthday problem
where I can assess the probability of n students having the same test
score in a class of size k.
I was writing my own code for the b-day problem until I r
I tried to modify the R function bw.SJ (from the stats package) for my own
use. But if I just get the source code and run directly (without any
modification), it can no longer find some key functions called within it. I
understand this has something to do with searching path which I do not
un
"jim holtman" <[EMAIL PROTECTED]> wrote in message
news:[EMAIL PROTECTED]
> --max-mem-size=N
> (Windows only) Specify a limit for the amount of memory to be used
> both for R objects and working areas. This is set by default to the
> smaller of 1.5Gb24 and the amount of physical RAM in the machine
8rino-Luca Pantani wrote:
> I see.
> Thanks for your highly appreciated help, Prof. Ripley.
>
>
> I've been reading the x11 help file, only after people from this list
> told me that it is the default in Linux
>
> Please consider that I recently (and definitely) moved to linux from
> Windows, and
Charles C. Berry tajo.ucsd.edu> writes:
> As the posting guide suggests, you could perform a search using
>
> > RSiteSearch("repeated measures", restric="functions")
>
> say.
>
> That generates an inventory of functions that pertain to repeated
> measures designs.
The problem I noted in prac
On Tue, 6 Nov 2007, 8rino-Luca Pantani wrote:
> I see.
> Thanks for your highly appreciated help, Prof. Ripley.
>
>
> I've been reading the x11 help file, only after people from this list told me
> that it is the default in Linux
>
> Please consider that I recently (and definitely) moved to linux
For the survival package, you might want to reference appendix E of Therneau
and Grambsh, Modeling Survival Data, Springer-Verlag, 2000. This contains a
set
of small example data sets for which results have been worked out in closed
form. S-plus and R pass all the tests: they are a portion o
Pedro Mardones gmail.com> writes:
>
> I would like to try to implement a big series of nested loops in a C
> code and then call it from R; however I'm not familiar with C
> programming. Does anyone know about some sort of reference I can use
> that help me to translate my code to C? (I'm thinkin
I see.
Thanks for your highly appreciated help, Prof. Ripley.
I've been reading the x11 help file, only after people from this list
told me that it is the default in Linux
Please consider that I recently (and definitely) moved to linux from
Windows, and I still need to re-learn a lot of things
Thanks for the help.
On Tue, 2007-11-06 at 14:48 -0200, Alberto Monteiro wrote:
> Deepankar Basu wrote:
> >
> > I am trying to do an ML estimation in R. My likelihood function has
> > several nested loops and so it takes a lot of time (days when I use the
> > genetic algorithm for optimization) f
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Of course, Brian is correct that the problem requires a great
deal more specificity. However, one could imagine that it would
be useful to tag objects created by top-level expressions
with the time they were created in much the same way that the
file
Dear R list,
I would like to perform an ANOVA in a set of measurements, but I have
problems formatting the data.
The data is a two dimensional array containing two columns:
- "Stim" : the type of stimulation (string)
- "Ratio" : a ratio of two numeric values
Now, because some values are missing
Deepankar Basu wrote:
>
> I am trying to do an ML estimation in R. My likelihood function has
> several nested loops and so it takes a lot of time (days when I use the
> genetic algorithm for optimization) for the optimization to finish.
> Unable to avoid loops, I am thinking of writing the likeli
You didn't mention your platform but if its windows then try this:
plot(1:10)
for(p in c("ps", "pdf", "jpg")) savePlot("myplot", type=p)
and if not look up ?dev.copy and RSiteSearch("dev.copy")
On Nov 6, 2007 9:18 AM, Cuvelier Etienne <[EMAIL PROTECTED]> wrote:
> Is it possible to produce the s
On Tue, 6 Nov 2007, Deepankar Basu wrote:
> Hi All,
>
> I am trying to do an ML estimation in R. My likelihood function has
> several nested loops and so it takes a lot of time (days when I use the
> genetic algorithm for optimization) for the optimization to finish.
> Unable to avoid loops, I am
Yes, this can be done and the call is the same, however the C++ needs to
be indicated as , see details in the writing R extensions
manual.
hth, Ingmar
On 6 Nov 2007, at 17:05, Deepankar Basu wrote:
> Hi All,
>
> I am trying to do an ML estimation in R. My likelihood function has
> several neste
On 06-Nov-07 15:53:53, Oarabile Molaodi wrote:
> I am trying to determine whether two samples are identical or not.
> I'm aware that somebody can use the Kolmogorov-Smirnoff test to
> compare empirical distributions, but since my samples have ties
> I'm not sure if I'm getting the right p-values fo
On Linux the default graphics device is X11(), and not so on Windows.
?X11 tells you no less than *two* ways to do this: I do wonder why you
don't just read the help file?
?options (see 'device') tells you how to set up a default device which you
could use to change default arguments, but using
I would like to try to implement a big series of nested loops in a C
code and then call it from R; however I'm not familiar with C
programming. Does anyone know about some sort of reference I can use
that help me to translate my code to C? (I'm thinking on something
like to Octave to R reference po
-- Forwarded message --
From: Mark Wardle <[EMAIL PROTECTED]>
Date: 6 Nov 2007 16:01
Subject: Re: [R] Produce a multiple formats graphic
To: Cuvelier Etienne <[EMAIL PROTECTED]>
You have several options:
1. Write a function that creates your plot, and call that function
repeatedl
I am trying to determine whether two samples are identical or not. I'm
aware that somebody can use the Kolmogorov-Smirnoff test to compare
empirical distributions, but since my samples have ties I'm not sure if
I'm getting the right p-values for the comparison. Can the
Kolmogorov-Smirnoff test
Hello All,
I trying to find some way to fill in missing values in a seasonal time
series. All the function that I find until now, don't have any reference to
seasonal data and the output is very different of what I looking for.
I also searched the forum but this problem don't have many informatio
Hi All,
I am trying to do an ML estimation in R. My likelihood function has
several nested loops and so it takes a lot of time (days when I use the
genetic algorithm for optimization) for the optimization to finish.
Unable to avoid loops, I am thinking of writing the likelihood function
in C++ and
Cuvelier Etienne wrote:
> Is it possible to produce the same graphic in different formats
> in one shot?
> Now if I want to produce the same graphic in pdf, ps and jpeg
> I run the producing code, but I change the graphic device, and
> I choose successively :
> pdf(...)
> ps(...)
> jpeg(...)
>
>
Hi Shubha,
Have a look at the Hmisc package and try ?wtd.mean as
the help command. This may be what you need. Also you
might want to look at ?summerize in the Hmisc package.
--- Shubha Vishwanath Karanth
<[EMAIL PROTECTED]> wrote:
> Hi R,
>
>
>
> I have the data,
>
>
>
> Value
>
>
Is it possible to produce the same graphic in different formats
in one shot?
Now if I want to produce the same graphic in pdf, ps and jpeg
I run the producing code, but I change the graphic device, and
I choose successively :
pdf(...)
ps(...)
jpeg(...)
So in this case I run the same code 3 times,
Does anyone know if there exists a package for non-parametric tolerance
limits in R? Or any other way to do it?
--
View this message in context:
http://www.nabble.com/Non-parametric-tolerance-limits-tf4750263.html#a13583088
Sent from the R help mailing list archive at Nabble.com.
__
The lines are:
## Initial position of the graphics window
## (pixels, <0 values from opposite edge)
xgraphics = 200
ygraphics = 10
On 06/11/2007, 8rino-Luca Pantani <[EMAIL PROTECTED]> wrote:
>
> There's a Renviron file under /etc/R, but there are no lines with
> "position (x, y)".
>
> Henrique
Hi all,
Usually when I do a command in shell and have the output read from R, I
like to do this:
pp <- pipe(cmd)
sol <- scan(pp, what="character", sep="\n")
I prefer this one over read.table since read.table somewhere has a gc
call and it's slow.
However, sometim
Hi R,
I have the data,
Value
0
1
2
3
4
5
6
Frequency:
258
125
106
94
85
12
5
Say:
Value=c(0,1,2,3,4,5,6)
Frequency=c(258,125,106,94,85,12,5)
In R, how can I find the mean, variance, skewness and kurtosis for this
grouped data? Is there a direct way of handling th
On Tue, 6 Nov 2007, Diogo Alagador wrote:
> Hi all,
>
> Basically I have 2 data frames with equal dimension and I want to apply
> the wilcox.test to compare columns in the same position (i.e. 1st of
> df.A with 1st of df.B, 2nd of df.A with 2nd of df.B,...). Anyone give me
> an hint on this, as
There's a Renviron file under /etc/R, but there are no lines with
"position (x, y)".
Henrique Dallazuanna ha scritto:
> In the Windows the file Rconsole contains the position(x,y) of the graphics.
> I don't know if file have the same name in Ubuntu.
>
I don't think that the loop-part puts a great penalty in this case;
however, it'd be better to convert the data.frames to matrices, since
matrix indexing is faster, e.g.,
# I presume that you only want to keep the p-values
mat.A <- data.matrix(df.A)
mat.B <- data.matrix(df.B)
pvalues <- numeric(
On Tue, 2007-11-06 at 08:03 -0600, Dirk Eddelbuettel wrote:
> On Tue, Nov 06, 2007 at 01:08:40PM +0100, Sven Garbade wrote:
> > Dear list,
> >
> > I cannot load the tcltk library:
> >
> > > library(tcltk)
> > Loading Tcl/Tk interface ... Error in fun(...) :
> > Can't find a usable tk.tcl in t
On Tue, Nov 06, 2007 at 01:08:40PM +0100, Sven Garbade wrote:
> Dear list,
>
> I cannot load the tcltk library:
>
> > library(tcltk)
> Loading Tcl/Tk interface ... Error in fun(...) :
> Can't find a usable tk.tcl in the following directories:
> /usr/share/tcltk/tcl8.4/tk8.4 ./lib/tk8.4 /u
Hi all,
Basically I have 2 data frames with equal dimension and I want to apply the
wilcox.test to compare columns in the same position (i.e. 1st of df.A with 1st
of df.B, 2nd of df.A with 2nd of df.B,...).
Anyone give me an hint on this, as I think it is nicer to avoid loops,
specially for hu
Hello everybody!
I have problems with integrating my function.
My primary function is a survival function of the following type:
surviv <- exp(-k*x)/(1+exp(alpha*(x-tau)))
I would like to integrate this function over a defined range and obtain a
vector with all the values from integrate(survi
Hi,
you may want to look at Macdonald P (2006) Mixdist package for R
www.math.mcmaster.ca/peter/mix/mix.html
It also has an anova test for your model so you would know if your "possible
sub-populations" are statistically significant of not. I don't know why
actually this package is not pos
Thanks to Greg Snow and John Kane I now have a working function that
does what I wanted, that is compares two scores in a plot.
Here is the function:
## compare.ratings: plots two lists corresponding to two different
## ratings. For each element, a line connects the position of that
## element in
On 11/6/07, Van Campenhout Bjorn <[EMAIL PROTECTED]> wrote:
> Hi all,
>
> I made a dotplot() with lattice, which comes out nice on the graphics
> device. I can save this as a eps using postscript() and include this in
> a word document. This prints nice, but does not look good on screen.
> If I p
In the Windows the file Rconsole contains the position(x,y) of the graphics.
I don't know if file have the same name in Ubuntu.
On 06/11/2007, 8rino-Luca Pantani <[EMAIL PROTECTED]> wrote:
>
> Dear R users,
> is there a way to tell R where to draw the windows where the plots are
> drawn ?
>
> I a
Hi all,
I made a dotplot() with lattice, which comes out nice on the graphics
device. I can save this as a eps using postscript() and include this in
a word document. This prints nice, but does not look good on screen.
If I produce a pdf, it is nice on screen, but not on paper. How can I
save a
Thank you very much.
Worked perfectly !!
El dl 05 de 11 del 2007 a les 22:17 +, en/na Prof Brian Ripley va
escriure:
> It would help if we knew what that command was supposed to do. My guess
> is that you want the Tr$size with the smallest value of Tr$dev, in which
> case
>
> Tr$size[which
Dear R users,
is there a way to tell R where to draw the windows where the plots are
drawn ?
I always get that window on the right of my screen, where it interferes
with my Xemacs + ESS window where I'm writing lines.
Of course I can move it manually, but I would like to have it
permanently on
Dear list,
I cannot load the tcltk library:
> library(tcltk)
Loading Tcl/Tk interface ... Error in fun(...) :
Can't find a usable tk.tcl in the following directories:
/usr/share/tcltk/tcl8.4/tk8.4 ./lib/tk8.4 /usr/local/lib/tcltk/tk8.4
/usr/local/share/tcltk/tk8.4 /usr/lib/tcltk/tk8.4 /u
Quin Wills wrote:
> Hello all
>
> (1) My departmental IT guru has installed the latest Wine (0.9.48) on a
> well maintained Fedora 7 server, so that I can run WinBUGS from R via
> R2WinBUGS (this seems like the only real option for running BUGS from R
> on a Linux platform?)
> (2) WinBUGS see
Hello all
(1) My departmental IT guru has installed the latest Wine (0.9.48) on a
well maintained Fedora 7 server, so that I can run WinBUGS from R via
R2WinBUGS (this seems like the only real option for running BUGS from R
on a Linux platform?)
(2) WinBUGS seems to run okay, except for the fo
On Tue, 6 Nov 2007, Alberto Monteiro wrote:
>
> Prof Brian Ripley wrote:
>>
>>> I would like to know the created time and date of specific object.
>>> Is there any function for it?
>>
>> There isn't even the concept. Most objects in R are a collection of
>> SEXPRECs created at different times. S
The attachment(s) that you sent with the following mail
had Viruses in it!
=
The Mail came from: [EMAIL PROTECTED]
The Mail recipient: [EMAIL PROTECTED]
Subject of the Mail : Delivery reports about your e-mail
Message-ID
The attachment(s) that you sent with the following mail
had Viruses in it!
=
The Mail came from: [EMAIL PROTECTED]
The Mail recipient: [EMAIL PROTECTED]
Subject of the Mail : Delivery reports about your e-mail
Message-ID
Dear,
I am having problem importing data from other packages or spreadsheet
to R. The data file contains characters and numerics variables.
[EMAIL PROTECTED] wrote:
Send R-help mailing list submissions to
r-help@r-project.org
To subscribe or unsubscribe via the World Wide Web, visit
ht
Prof Brian Ripley wrote:
>
>> I would like to know the created time and date of specific object.
>> Is there any function for it?
>
> There isn't even the concept. Most objects in R are a collection of
> SEXPRECs created at different times. Suppose you create a data
> frame out of existing co
Hi again,
I found a solution (but maybe it is not the best one) : I converted the
poynom into a function.
pderiv<-as.function((deriv(polynomial(pfit$coefficients
uniroot(pderiv, c(0.8,0.2))[[1]]
[1] 0.4778413
Now I just have to find the corresponding "b" to this "a" but I think it
must be q
A point of clarification: the measurements are continuous; the road surface
is discrete with four levels, and the road capacity is discrete with five
levels. Each measured point has both a surface and capacity associated with
it.
Thanks,
Ari
I have a series of measurements
>
[[alterna
Hello,
I have 3 columns : a, b and a*b
I would like to find the pair (a,b) so that a*b is the minimum but not from
the points I measured but from the fit of the curve (I have more points that
the ones given below but I fit only on this part because I
know that the minimum a*b is in this interval
77 matches
Mail list logo