Hello,
I want to do a plot of two series, but adding a secondary axis in the rigth
side to take into account the different scales of the two series. Anyone can
tell me haw can I do that in R?
Thank you very much
Maria
--
View this message in context:
http://www.nabble.com/secondary-axis-in-a-plo
I have a hclust object that looks like this:
> test77
Call:
hclust(d = input)
Cluster method : complete
Number of objects: 11
> test77$height
[1] 0.000 0.000 0.000 0.000 0.000 0.000 0.900
[8] 0.9473684 1.7894737 8.5771948
> test77$merge
[,1] [,2]
[1,] -
Dear all -
With the following code I get Boxplots and Stripcharts in one Plot:
with(InsectSprays, boxplot(count ~ spray, boxwex = 0.3))
with(InsectSprays, stripchart(count ~ spray, col = "red", vertical =
TRUE, add = TRUE))
But the dots from the stripchart are plotted over the Boxes.
Is there any
Dear all -
With the following code I get Boxplots and Stripcharts in one Plot:
with(InsectSprays, boxplot(count ~ spray, boxwex = 0.3))
with(InsectSprays, stripchart(count ~ spray, col = "red", vertical =
TRUE, add = TRUE))
But the dots from the stripchart are plotted over the Boxes.
Is there any
MarySR:
> I want to do a plot of two series, but adding a secondary axis in the
> rigth side to take into account the different scales of the two series.
> Anyone can tell me haw can I do that in R?
I would *strongly* advice you to read this article on dual-scaled axes
before proceeding:
http://
MarySR wrote:
Hello,
I want to do a plot of two series, but adding a secondary axis in the rigth
side to take into account the different scales of the two series. Anyone can
tell me haw can I do that in R?
Hi Maria,
Have a look at twoord.plot in the plotrix package.
Jim
_
If you need to do it:
plot(PLOT1, bty='c')
par(new=TRUE)
plot(PLOT2, axes=FALSE, ylab='', xlab='', bty='c')
axis(4)
mtext("secondary y-axis label",4)
On Thu, May 22, 2008 at 7:16 AM, MarySR <[EMAIL PROTECTED]> wrote:
>
> Hello,
> I want to do a plot of two series, but adding a secondary axis in
Not exactly clear what you are asking for. Your data.frame.table does not
seem related to the original 'orig'. What exactly are you expecting as
output?
On Wed, May 21, 2008 at 10:16 PM, maiya <[EMAIL PROTECTED]> wrote:
>
> i appologise for the trivialness of this post - but i've been searching
Dieter Wirz:
> with(InsectSprays, boxplot(count ~ spray, boxwex = 0.3))
> with(InsectSprays, stripchart(count ~ spray, col = "red", vertical =
> TRUE, add = TRUE))
>
> But the dots from the stripchart are plotted over the Boxes.
> Is there any possibility to have Stripchart and Boxplots side-by-s
sorry, my mistake!
the data frame should read:
orig<-as.data.frame.table(orig)
orig
Var1 Var2 Freq
1AA 40
2BA5
3AB 30
4BB 25
but basicaly i would simply like a sample of the original matrix ( which is
a frequency table/contingency table/crosstabulation)
On May 22, 2008, at 8:56 AM, maiya wrote:
sorry, my mistake!
the data frame should read:
orig<-as.data.frame.table(orig)
orig
Var1 Var2 Freq
1AA 40
2BA5
3AB 30
4BB 25
but basicaly i would simply like a sample of the original matrix
( which is
a fre
Is this what you want?
> xtabs(Freq ~ Var1 + Var2, data = orig)
Var2
Var1 A B
A 40 30
B 5 25
See ?xtabs
Or is this what you want?
expand.dft <- function(x, na.strings = "NA", as.is = FALSE, dec = ".")
{
DF <- sapply(1:nrow(x), function(i) x[rep(i, each = x$Freq[i]), ],
I do this all the time. See
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/105099.html
(Vista is a bit harder than XP, but similar.)
As far as determining when to proceed with the script,
you could have your db updater create a semaphore file your batch or R
script could wait for.
HTH,
David L.
Hi, I am fairly new to R so this may be easy for some.
I have one column that is the Type="A, B' T"
and I have another column that lists Price='894, 895, 896'
I wanted to plot three different series based upon the Type...If B than
xts..If T yts...
How exactly do I do this?
Thanks in advance
--
The weights given should correspond to the ordering of levels(y) where y
contains the class labels. If in doubt, you can also give the classwt
as a named vector (e.g., classwt=c(B=3, A=2, C=1)).
Search in the R-help archive to see other options and why you probably
shouldn't use classwt.
Andy
F
Hello,
Does anyone know what the functions p_length and s_length do in SPLUS and
what is their equivalent in R?
Thanks!
--
View this message in context:
http://www.nabble.com/length-functions-tp17406498p17406498.html
Sent from the R help mailing list archive at Nabble.com.
___
On Tue 20 May 2008 at 11:52:32 AM -0400, George Georgalis wrote:
>
>But this makes me realize, since our OS is generally
>backwards compatible, I may be able to build on the old
>systems and use on the newer ones as well. Will give
>that a shot.
>
that doesn't work, I think because gcc changed fro
Hi,
I am doing a very informal presentation for my office about R capabilities to
deal with and analyze spatial data, display data and maps, and connections with
GIS. I've used in my presentation info from the CRAN, the spatial Task view,
and the more striking graphics examples from
http://ad
On Thu, May 22, 2008 at 5:00 PM, Monica Pisica <[EMAIL PROTECTED]> wrote:
> - It seems that it is hard to install correctly under Linux.
Actually, it is quite easy to install R under Linux, at least in some
distributions. For instance, on Fedora:
yum install R R-devel
and that is it.
Paul
Monica Pisica wrote:
>
> Cons:
>
> - R has a very steep learning curve.
I don't think the learning curve is any steeper than SAS programming,
it is just a different kind of curve.
--
Kevin E. Thorpe
Biostatistician/Trialist, Knowledge Translation Program
Assistant Professor, Department of Pub
agree,
i think R is more like a standard program language than SAS. however,
SAS programmer might not feel intuitive to pick up R.
On Thu, May 22, 2008 at 12:14 PM, Kevin E. Thorpe
<[EMAIL PROTECTED]> wrote:
> Monica Pisica wrote:
>>
>> Cons:
>>
>> - R has a very steep learning curve.
>
> I don't
Wensui Liu wrote:
agree,
i think R is more like a standard program language than SAS. however,
SAS programmer might not feel intuitive to pick up R.
That says more about SAS programmers than it does about R.
__
R-help@r-project.org mailing list
htt
Hello,
I enjoy using R, and find the excellence of this newsgroup removes most of the
altitude of the "learning curve". Any learning curve, IMHO, is more than offset
by the benefits and capabilities of the language and its awesome plotting tools.
A grateful advocate,
John
Kevin E. Thorpe wro
What do you mean, "there are no up-grades?" There are 1,401 ancillary
packages - all free. That sounds like upgrades to me.
Of course there is no *perfect* beginner's book, but Peter Dalgaard's
Introductory Statistics with R (Paperback), Springer, 3d printing edition
(January 9, 2004) is pretty
Hi ngupta:
i didn't get what you mean with "If B than
xts..If T yts..."
i suppose that you want to plot the prices of type A, then the prices of type
B and so on...
# well, to extract prices of tipe A, you can use
A <- which(table[,1] == "A")
# and then
A.prices <- table[,2][A]
hope that helps,
On the 64-bit part, I tested 2.7.0 on a Dual Core Lenovo ThinkPad and was able
to allocate memory beyond 2G. Have not done much else otherwise but it seems to
work just fine.
H
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Monica Pisica
Sent: Thursday
Marc, it's the second "expansion" type transformation I was after, although
your expand.dft looks quite complicated? here's what I finaly came up with -
the bold lines correspond to what expand.dft does?
> orig<-matrix(c(40,5,30,25), c(2,2))
> orig
[,1] [,2]
[1,] 40 30
[2,]5 25
>
On Thu, 22 May 2008, francogrex wrote:
f >
f > Hi, this is probably quite simple but I can't seem to do it correctly. I
have
f > a data frame of counts of infections in different ages; something like:
f > count=c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 7, 8, 8, 9, 9,
f > 10, 11, 15, 17,
Hello,
I'm a new user to R, and I was trying to install "packfor" in my R console
and I always run into an error.
I tried to search, in in the forum, someone that already might have made the
same question, but all I could find was a general information how to install
libraries in R.
Pass through
This is about the 8761st time that this old warhorse has been discussed on
this list. I suggest you search the archives (e.g. via CRAN, Gmane,
Google,...) for these prior discussions.
-- Bert Gunter
Genentech
P.S. Incidentally, the answer to your question is "It depends..." -- on what
your appli
Dear list,
A new mailing list *R-sig-genetics* has just been created (thanks to
Martin Maechler).
This list is dedicated to discussions about the *analysis of genetic
markers* using the R software. The covered topics, along with various
information about the list, are available from the R-si
Hello,
I am trying to run some AMOVA analyses with the amova function in the
package ade4.
When running the example dataset provided in ade4, I noticed a
difference between the published results from the same data (Excoffier
et al. 1992) and what ade4 calculates.
Below are the data for "wi
Charles,
I don't consider packages as up-grades, i consider them as enhancements .
When a new R version is in place you cannot up-grade your old R version, you
have to do a new installation and re-load all the packages you used to have and
delete / un-install the old version of course
Thanks for the heads up - i will add that to the "pros".
Monica> From: [EMAIL PROTECTED]> To: [EMAIL PROTECTED]; r-help@r-project.org>
Date: Thu, 22 May 2008 09:32:06 -0700> Subject: RE: [R] Pros and Cons of R> >
On the 64-bit part, I tested 2.7.0 on a Dual Core Lenovo ThinkPad and was able
t
Hi,
Monica Pisica wrote:
> - There is no perfect “beginner” book.
How about
- Crawley, Michael (2007). The R book, Wiley & Sons.
- Maindonald, John & John Braun (2007): Data Analysis and Graphics Using
R (2nd edition), Cambridge University Press.
As a political scientist (with programming exp
Kim Milferstedt uiuc.edu> writes:
>
> Hello,
>
> I am trying to run some AMOVA analyses with the amova function in the
> package ade4.
>
> When running the example dataset provided in ade4, I noticed a
> difference between the published results from the same data (Excoffier
> et al. 1992) a
> Hi,
>
> I am doing a very informal presentation for my office about R capabilities
> to deal with and analyze spatial data, display data and maps, and
> connections with GIS. I've used in my presentation info from the CRAN, the
> spatial Task view, and the more striking graphics examples from
> h
I think it is a fairly safe bet that at any given time there
are numerous 64-bit R sessions using more than 4GB
of memory -- I periodically have one.
In my opinion 64-bit implementation is definitely in the
"pro" column. It is as transparent going from 32-bit to
64-bit as moving between 32-bit pl
Hi,
IMHO, "- Everything you want to do is a command line, minimal GUI", this should
be considered a "Pro". The commandline is much more powerful than a GUI.
Cheers,
John
Monica Pisica wrote:
>
> Charles,
>
> I don't consider packages as up-grades, i consider them as enhancements .
>
Hi,
I have a table with three columns of data, let's call these a, b and
c. I have produced a scatter plot of a against b, but I would like to
colour the points according to their corresponding c values. c takes
the value of either 0,1,2,3 or 4 and so I would like to colour all of
the poi
I have an equation describing the best-fit model for a set of points (just 2
axes) that is in the form:
y=b+mx+px^2
Where b is the intercept, m is the slope describing a linear term, and p is
a slope of the quadratic term.
I would like to plot this equation on a curve (I know the equation is
y=(
Colleagues,
I am executing a large batch file in R (2.7.0) in OS X (Leopard) or
Linux (Fedora).
Usually, R completes normally and executes wrap-up commands. However,
occasionally, R is "killed" externally, e.g., "kill " from the
OS. When that happens, the wrap-up is not executed.
Is
There is a new book coming out this year on spatial data analysis in R.
http://www.springer.com/statistics/statistical+theory+and+methods/book/978-0-387-78170-9
I am not sure there are any other book length treatments of spatial data
analysis that utilize any of the other major statistical pa
On Thu, 22 May 2008, Monica Pisica wrote:
I don't consider packages as up-grades, i consider them as
enhancements . When a new R version is in place you
cannot up-grade your old R version, you have to do a new
installation and re-load all the packages you used to have
and delete / un-inst
Hi,
I have an object from the prediction function from the ROCR package and I would
like to extract one of the slots from the object, for example the cutoffs slot.
However the usual techniques ($, [["name"]]) of subsetting don't work. How can
I assess the lists in the slots?
Here is an example
Hi Ryan,
See the example in ?predict.lm . I think it would be useful in this case.
HTH,
Jorge
On Wed, Dec 31, 1969 at 8:16 PM, Ryan Harrigan <[EMAIL PROTECTED]> wrote:
> I have an equation describing the best-fit model for a set of points (just
> 2
> axes) that is in the form:
>
> y=b+mx+px^2
On Thu, May 22, 2008 at 12:00 PM, Monica Pisica <[EMAIL PROTECTED]> wrote:
>
> Hi,
>
> I am doing a very informal presentation for my office about R capabilities to
> deal with and analyze spatial data, display data and maps, and connections
> with GIS. I've used in my presentation info from the
Hi Stacey,
Try this:
library(ROCR)
data(ROCR.simple)
pred <- prediction(ROCR.simple$predictions,ROCR.simple$labels)
perf <- performance(pred,"sens","spec")
cuts <- unlist(slot(perf,"alpha.values"))
cuts
HTH,
Jorge
On Thu, May 22, 2008 at 2:08 PM, Stacey Burrows <[EMAIL PROTECTED]>
wrote:
> Hi
Hi,
I have no idea about SAS but i know that when i've started using R i've
struggled quite a bit until i started to get some insight and results. I
cannot say it is easier now, but certainly i like it much better, i would not
be ale to do the work i am doing without it and i would hate
On 5/22/08, Sarah Kinnings <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I have a table with three columns of data, let's call these a, b and c. I
> have produced a scatter plot of a against b, but I would like to colour the
> points according to their corresponding c values. c takes the value of
> either
On Thu, May 22, 2008 at 02:07:01PM -0400, R P Herrold wrote:
> On Thu, 22 May 2008, Monica Pisica wrote:
[...]
> >When a new R version is in place you
> >cannot up-grade your old R version, you have to do a new
> >installation and re-load all the packages you used to have
> >and delete / un-in
Hi Stacey,
ROCR uses S4 classes. The elements are accessed using "@" instead of "$".
You can find an example on slide 12 of the following slide deck:
http://rocr.bioinf.mpi-sb.mpg.de/ROCR_Talk_Tobias_Sing.ppt
Also have a look at the R code that appears when you type "demo(ROCR)" in R
which contai
Never mind! In fact I just discovered that the underscore and the "<-" are
used equivalently in SPLUS.
Applejus wrote:
>
> Hello,
>
> Does anyone know what the functions p_length and s_length do in SPLUS and
> what is their equivalent in R?
>
> Thanks!
>
--
View this message in context:
Has anyone come across any alternatives to rJava and JRI? Are they any
good? Better perhaps?
Please give your reasons.
Thanks,
Danish
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- - - -
This message is intended only for the personal and confidential use of t
I do not know if this is "better" as we are just starting to look
into it (and "better" is in the eye of the beholder) but look at:
http://tolstoy.newcastle.edu.au/R/e3/help/07/11/4253.html
in the archives. The link to the software is http://www.ebi.ac.uk/
microarray-srv/frontendapp/
HTH,
There's SJava, http://www.omegahat.org/RSJava.
On Thu, May 22, 2008 at 12:50 PM, Munir, Danish <[EMAIL PROTECTED]>
wrote:
>
> Has anyone come across any alternatives to rJava and JRI? Are they any
> good? Better perhaps?
>
> Please give your reasons.
>
> Thanks,
> Danish
> - - - - - - - - - - - -
Although this is an update of the same code base, RServe is nice and has
newer objects that one may use. For example, the exception handling has been
improved and now there is a external server ability. Before you give up on
JRI look into the newer code, it might be more useful.
http://www.rforge.
Dear All,
I would be most grateful for help on this problem.
I have a binary response variable. It is classified by various factors.
I would like to have a table with;
column 1
frequency response is 1 / (frequency the response is 1 + frequency the
response is 0) #given the classifying factors
Hello,
maybe someone can help me. I am looking for a possibility to plot a 3D area
diagram like in Excel:
http://www.microsoft.com/germany/mac/excel/images/chartbefore.jpg Watch
this!
Would be nice if someone had any idea about that.
Thank you
Alex
--
View this message in context:
http://www
Hi all, I tried to follow an online tutorial to run openBUgs but the package
BRugs has been removed from R repository.
Could someone provide a link for where to download BRugs?
Thanks a lot!
--
View this message in context:
http://www.nabble.com/Where-to-download-BRugs-tp17411410p17411410.htm
Hi,
While executing the make after successful configuration (./configure
--with-readline=no --without-iconv), the following error occurs:
mbf2n11s (Regatta) /saswork/R/R-2.7.0$ make
Target "R" is up to date.
Target "R" is up to date.
Target "R" is up to date.
Target "R" is up to date.
Target "Ma
Hi,
Ben is correct. The alternative hypothesis is not good for this test
(within pop). I will modify it and submit a new version of the package
on CRAN (probably next week).
BTW, for question about ade4 please send your email to adelist
(http://listes.univ-lyon1.fr/wws/info/adelist)
Cheers
On 22/05/2008 3:04 PM, AlGates wrote:
Hello,
maybe someone can help me. I am looking for a possibility to plot a 3D area
diagram like in Excel:
http://www.microsoft.com/germany/mac/excel/images/chartbefore.jpg Watch
this!
Would be nice if someone had any idea about that.
I'd say the world w
On May 22, 2008, at 3:04 PM, AlGates wrote:
Hello,
maybe someone can help me. I am looking for a possibility to plot a
3D area
diagram like in Excel:
http://www.microsoft.com/germany/mac/excel/images/chartbefore.jpg
Watch
this!
Would be nice if someone had any idea about that.
I don't
AlGates algates.de> writes:
>
>
> Hello,
>
> maybe someone can help me. I am looking for a possibility to plot a 3D area
> diagram like in Excel:
> http://www.microsoft.com/germany/mac/excel/images/chartbefore.jpg Watch
> this!
>
This would be theoretically possible, although not that eas
Hello,
I want to plot 21 scatter plots in a 7*3 matrix on a single A4 page
(using mfcol). Below is an example, which -by and large- looks very
close to the desired result. But I do not get it to fill out the whole
page. If I specify the a different page size, the plots still do not
fill out the wh
Could you mean RBugs?
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of qqw
Sent: Thursday, May 22, 2008 3:10 PM
To: r-help@r-project.org
The BRugs package is maintained by Uwe Ligges. So it is possible he
forgot to place it in the new version of R/ repositories.
Aldi
Charles Annis, P.E. wrote:
Could you mean RBugs?
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.
Thank you all for your help. I will check the packages you suggested asap and
maybe I find a good way to get a such a plot.
The picture I posted was only an example taken from the microsoft-website,
just to make clear what kind of diagram I am looking for. I know that one
part is hidden by the ot
Hi,
To continue the test of the differences between mixed model SAS (9.1.3)
and lme (R 7.0), I removed the intercept of the random effect in the R
case. In that case lme is showing that beta coefficients are almost all
of them identical (-0.08424262) and somehow all of them are spread along
the
Hi,
I'm having difficulties with Unicode characters with R 2.7.0 on Windows
XP SP2. I'd like to use a delta character, which I'd also like to mix
with other formatted output using paste() or sprintf(). On my system, I
can successfully do this with a static string:
print("\u0394Q = 2.2 L/s")
Hi,
Below I have a function mlogl_k,
later it's called with "nlm" .
__BEGIN__
vsamples<- c(14.7, 18.8, 14, 15.9, 9.7, 12.8)
mlogl_k <- function( k_func, x_func, theta_func, samp) {
tot_mll <- 0
for (comp in 1:k_func) {
curr_mll <- (- sum(dgamma(samp, shape = x_func,
scale=theta_
Is this what you want:
> x <- by(frame, list(frame$p.a, frame$p.b), function(.row){
+ c(f=sum(.row$r)/ nrow(.row), p.a=.row$p.a[1],
p.b=.row$p.b[1],freq=nrow(.row))
+ })
> do.call(rbind, x)
f p.a p.b freq
[1,] 0.50 10 22
[2,] 0.50 20 32
[3,] 0.25 50 54
On Thu, M
Dear all,
I used metaMDS (vegan) to generate two NMDS axis and now I would like
to know the amount of variance explained for each axis. Is there a way to do
that?
kind regards,
miltinho
[[alternative HTML version deleted]]
__
R-help@r-project
On May 22, 2008, at 9:14 PM, Edward Wijaya wrote:
Hi,
Below I have a function mlogl_k,
later it's called with "nlm" .
__BEGIN__
vsamples<- c(14.7, 18.8, 14, 15.9, 9.7, 12.8)
mlogl_k <- function( k_func, x_func, theta_func, samp) {
tot_mll <- 0
for (comp in 1:k_func) {
curr_ml
Thanks Aldi, do you know anywhere we could download previous package of
BRugs? Do you have a copy which you could send to me? (I only found one old
version on Splus site which is not usable under R2.6 :()
Thanks!
--
View this message in context:
http://www.nabble.com/Where-to-download-BRugs-tp1
Dear Mark and Haris,
Thanks so much to you both. I really appreciate it.
- Edward
On Fri, May 23, 2008 at 10:54 AM, <[EMAIL PROTECTED]> wrote:
> to add to what Haris said: I think you are trying to find the mle's of the
> gamma, given some data, and either what you're doing is not the right
On 22/05/2008 8:35 PM, qqw wrote:
Thanks Aldi, do you know anywhere we could download previous package of
BRugs? Do you have a copy which you could send to me? (I only found one old
version on Splus site which is not usable under R2.6 :()
Thanks!
You can get older versions from CRAN, in
http:/
Hi,
I just installed it from
http://cran.wustl.edu/bin/windows/contrib/2.7/BRugs_0.4-1.zip
Change 2.7 to 2.6 and you get the older version.
HTH,
Aldi
>
> Thanks Aldi, do you know anywhere we could download previous package of
> BRugs? Do you have a copy which you could send to me? (I only foun
Hi all,
What's the construct to import a built-in package
(e.g. ARIMA) into the R script.
I can't seem to find them in the documentation.
Regards,
Edward
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
R has very good GUI packages that I have used -
R Commander RCmdr and Rattle rattle . Since I work on multiple packages with
constraints of time, I almost always use the GUI rather go through the
intricacies of command line .
The log of these GUIs shows the relevant R command that was used, so you
Completely agreed!
BR, Shubha
Shubha Karanth | Amba Research
Ph +91 80 3980 8031 | Mob +91 94 4886 4510
Bangalore * Colombo * London * New York * San José * Singapore *
www.ambaresearch.com
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ruben Roa Ureta
Look like you need to use a Fortran compiler set to accept Fortran 77
(which the R sources are).
We don't know what 'f95' is, and this is not really an R issue so please
ask your sysadmins for help.
On Thu, 22 May 2008, Senthilkumaran wrote:
Hi,
While executing the make after successful c
install.packages("BRugs") still works on Windows -- the binaries have been
moved to the CRAN extras collections. If you want a direct URL
http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.7/BRugs_0.4-1.zip
On Thu, 22 May 2008, qqw wrote:
Hi all, I tried to follow an online tutorial
Dear all,
I am trying to model number of samples from
a given series. The series are modelled according
Gamma function.
In order to estimate the # samples, I use BIC/AIC
with MLE (computed from dgamma function).
Here is the code I have.
__BEGIN__
mlogl <- function( x_func, theta_func, samp) {
Hi!!!
I am trying to run a bootstrap but I have troubles. My code is as follows
theta <- function(x){mean(x,)}
boot1<-boot(x,theta,R=1000)
Error in statistic(data, original, ...) : unused argument(s) (1:487)
Any help is welcome
Thanks
[[alternative HTML version deleted]]
Dear List,
Does there exist a function that calculates a cumulative average?
Neither running() from library(gregmisc) nor running.mean() from
library(igraph) seems to be able to give a cumulative average.
Any help or pointers would be greatly appreciated.
Jacques
___
i appologise for the trivialness of this post - but i've been searching the
forum wothout luck - probably simply because it's late and my brain is
starting to go..
i have a frequency table as a matrix:
orig<-matrix(c(40,5,30,25), c(2,2))
orig
[,1] [,2]
[1,] 40 30
[2,]5 25
i basic
Hi, this is probably quite simple but I can't seem to do it correctly. I have
a data frame of counts of infections in different ages; something like:
count=c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 7, 8, 8, 9, 9,
10, 11, 15, 17, 17, 17, 17, 19, 19, 19, 19, 20, 20, 20, 21, 21,
21, 22, 22, 22
On 5/21/2008 4:52 PM, Sharai Gomez-Raygoza wrote:
Hi!!!
I am trying to run a bootstrap but I have troubles. My code is as follows
theta <- function(x){mean(x,)}
boot1<-boot(x,theta,R=1000)
Error in statistic(data, original, ...) : unused argument(s) (1:487)
Any help is welcome
Thanks
When s
does e.g.
x <- c(8,5,4,2)
cumsum(x) / 1:length(x)
do what you need?
Albert
Am Mittwoch, den 21.05.2008, 20:36 -0500 schrieb Jacques Wagnor:
> Dear List,
>
> Does there exist a function that calculates a cumulative average?
> Neither running() from library(gregmisc) nor running.mean() from
> libra
Hi,
how to estimate a the following model in R:
y(t)=beta0+beta1*x1(t)+beta2*x2(t)+...+beta5*x5(t)+beta6*y(t-1)+beta7*y(t-2)+beta8*y(t-3)
1) using "lm" :
dates <- as.Date(data.df[,1])
selection<-which(dates>=as.Date("1986-1-1") &
dates<=as.Date("2007-12-31"))
dep <- ts(data.df[selection,c("dep
Is this what you want:
> count=c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 7, 8, 8, 9, 9,
+ 10, 11, 15, 17, 17, 17, 17, 19, 19, 19, 19, 20, 20, 20, 21, 21,
+ 21, 22, 22, 22, 22, 23, 23, 23, 23, 23, 24, 24, 25, 27, 31, 33
+ )
> age=c(3, 8, 9, 7, 56, 58, 10, 13, 53, 55, 11, 12, 14, 51, 54,
+ 15,
Please do not send HTML (as the posting guide asks) -- it has made your
code unreadable.
Also, without seeing an example we don't know if the 'completely
different' is in the output or corresponds to very different fitted
models.
There is no 'armaFit' in R: it seems you may be using package
If you just want to bin the current ages:
> count=c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 7, 8, 8, 9, 9,
+ 10, 11, 15, 17, 17, 17, 17, 19, 19, 19, 19, 20, 20, 20, 21, 21,
+ 21, 22, 22, 22, 22, 23, 23, 23, 23, 23, 24, 24, 25, 27, 31, 33
+ )
> age=c(3, 8, 9, 7, 56, 58, 10, 13, 53, 55, 11, 12,
Hi,
I was wondering if there is a way to figure out why in SAS random beta
coefficients are 0 vs. in R the beta-s are non zero.
The variables of the data are nidl, time, and sub (for subject). Time and
nidl are continuous variables. I am applying random coefficients model.
Any input is greatly app
96 matches
Mail list logo
