Displaying 20 results from an estimated 800 matches similar to: "r: intergrate behaviour"
2009 Jan 02
1
R: numerical integration problems
hello all
happy new year and hope you r having a good holiday.
i would like to calculate the expectation of a particular random variable and would like to approximate it using a number of the functions contained in R. decided to do some experimentation on a trivial example.
example
========
suppose x(i)~N(0,s2) where s2 = the variance
the prior for s2 = p(s2)~IG(a,b)
so the posterior is
2010 Feb 11
1
Fwd: Recall function: "evaluation nested too deeply: infinite recursion / options(expressions=)?"
hello all
i searched the archives and couldn't get a solution to the following question.
i have the following function:
F=function(z,v)
{
if (v==-.5) return(1) else return(2*v/z + 1/Recall(z,v-1))
}
and while testing whether it works (ie F(z,v) is approx = besselK(z,1+v)/besselK(z,V). the recursion formula allows one to calculate besselK(z,1+v)/besselK(z,V) for large values of z )
2008 Jan 22
1
R: determinants and inverses
hello all
sorry for the following "none" R related question.
does anyone know of a reference to calculate the following identity:
|I + ABC|
where I is an identity matrix and A, B,C may not have to be square matrices?
you help will be greatly appreciated.
H. V. Henderson; S. R. Searle
SIAM Review, Vol. 23, No. 1. (Jan., 1981), pp. 53-60.
provides a result to
2007 May 02
1
missing values
hello,
I need your help for this example
> for(k in LR) {
+ donGeno[[k]] <- as.numeric(levels(factor(subset(don2, Id_Essai == 1006961 & Id_Cara == LC[1] & Id_Rep == k, select = Id_Geno)[,1])))
+ print(donGeno[[k]])}
[1] 65125 65126 65127 65128 65129 65130 65131 65132 65133 65134 65135 65136 65137 65138 65139 65140 65141 65142 65143 65144 65171
[1] 65126 65127 65128 65129 65130
2010 Feb 11
0
Recall function: "evaluation nested too deeply: infinite recursion / options(expressions=)?"
hello all
i searched the archives and couldn't get a solution to the following question.
i have the following function:
F=function(z,v)
{
if (v==-.5) return(1) else return(2*v/z + 1/Recall(z,v-1))
}
and while testing whether it works (ie F(z,v) is approx = besselK(z,1+v)/besselK(z,V). the recursion formula allows one to calculate besselK(z,1+v)/besselK(z,V) for large values of z )
2010 Jan 18
0
linking r and C
hello
i have a few quick questions.
i am using a windows platform.
how does one include functions from a particular r package in the C code. e.g. rgamma can be included if we include #include <R.h>
#include <Rmath.h> , and #include <math.h> in the C code. what should be included if i want to use functions that are inside the MASS package (for example)?
When performing mcmc
2009 May 05
1
Intergrate Jasper report to Ruby on Rails Application
Hi,
Im very new in Ruby on Rails. I have to do intergrate jasper report
with Rails applications. I found the link
http://wiki.rubyonrails.org/rails/pages/HowtoIntegrateJasperReports.
But Im not clear how to do. Plz help me this how to do. If anybody
have sample application it will be very better. I look forward your
reply...
Thanks
Velu
2006 Apr 26
1
Intergrate Asterisk IP PBX with Legacy PBX, continuing existing funtionality of legacy pbx
Hi All,
I would like to explain the layout that i am trying to achive. I am so
helpless on this regard.
So here is the story ........
" This is with regard to the setup which you can find at the
"Asterisk The Future of Telephony" , chapter 11, page # 196-197, I am
attaching the picture for your information.
Now I am taking a challenging step to of integrate IP PBX with our
2010 Nov 25
1
coxph strange result
The following fit does not make sense to me, please, correct me if I have a logical error.
> moddowsn
Call:
coxph(formula = Surv(start, stop, resp) ~ sn + matfac2, data = coxsn1,
method = "efron")
coef exp(coef) se(coef) z p
sn2 0.0497 1.051 0.02030 2.450 1.4e-02
sn3 -0.0532 0.948 0.02038 -2.610 9.0e-03
sn4 -0.0410 0.960
2007 Sep 20
2
referencing packages?
Hi
I know how to referenc R in a scientific paper - but is there a
standardised way to reference packages?
Thanks
Rainer
--
NEW EMAIL ADDRESS AND ADDRESS:
Rainer.Krug at uct.ac.za
RKrug at sun.ac.za WILL BE DISCONTINUED END OF MARCH
Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation
Biology (UCT)
Plant Conservation Unit
Department of Botany
University of Cape Town
Rondebosch 7701
2007 Oct 01
4
Concatenating one character vector into one string
Hi
I am sure this is simple - but how can I convert one charecter vector
into one string?
example:
x <- c("This ", "is ", "one ", "sentence.")
should become
"This is one entence"
Thanks
Rainer
--
NEW EMAIL ADDRESS AND ADDRESS:
Rainer.Krug at uct.ac.za
RKrug at sun.ac.za WILL BE DISCONTINUED END OF MARCH
Rainer M. Krug, Dipl. Phys.
2004 Oct 16
3
Cox PH Warning Message
Hi,
Can anybody tell me what the message below means and how to overcome it.
Thanks,
Neil
Warning message:
X matrix deemed to be singular; variable 2 in: coxph(Surv(age_at_death,
death) ~ project$pluralgp + project$yrborn + .........
>
2008 Mar 14
2
Selecting elements in vector
Hi
Consider the following code
> x <- rep(1:13, 13)
> y <- 1:3
I want to select all elements in x which are equal to 1, 2 or 3.
I know that I could use
> sel <- x==y[1] | x==y[2] | x==y[3]
> x[sel]
to obtain the values, but in my analysis, the y-vector is thousands of
elements long.
Is there any way, that I can do that easily?
Thanks
Rainer
--
Rainer M. Krug, Dipl.
2009 Feb 16
0
odd GARCH(1,1) results
Hi everybody,
I'm trying to fit a Garch(1,1) process to the DAX returns. My data
consists of about 2300 10day-logreturns in chronologically descending
order (see attachment). But if I use the garch function I get a very
high alpha_1 and a quite low beta, which doesn't make that much sense. I
think I am missing something, but have no idea what it might be. I'd
appreciate it a lot
2004 Mar 10
0
Newbie: Samba Intergration into W2K Domain/ADS
Hello everyone. I humbly admit I need some direction. I have search newsgroups and the Samba help files, but I must be missing something or just not getting it.
I am comfortable witrh Samba in a small workgroup environment, I feel like a fish out of water when trying to intergrate Samba into a W2K Domain/Active Directory environment. While I consider myself technical, I've never had to manage
2007 Aug 16
6
several plots on several pages
Hi
> version
_
platform i686-pc-linux-gnu
arch i686
os linux-gnu
system i686, linux-gnu
status
major 2
minor 5.1
year 2007
month 06
day 27
svn rev 42083
language R
version.string R version 2.5.1 (2007-06-27)
I want to create a pdf withe three graphs on a page and with two pages:
2007 May 28
1
Where to find "nprq"?
Hi
I am trying to install the package "pheno", but it needs the package
"nprq" by Roger Koenker et al. which I can I find this package? It does
not seem to be on CRAN and googling also doesn't give me an URL - is it
still somewhere available?
Thanks,
Rainer
--
NEW EMAIL ADDRESS AND ADDRESS:
Rainer.Krug at uct.ac.za
RKrug at sun.ac.za WILL BE DISCONTINUED END OF
2008 Mar 07
1
How to navigate in layout() created graph?
Hi
I created a complex layout with using layout() and it looks exactly as
I need it. But I don't want to print in the order in which the
subfigure are numbered, but in a different order.
How can I navigate in the layout so that I can specify the subfigure
in which to plot?
At the moment I am using a function which is converting the number to
mfg parameter for par, but it does not seem to
2008 Jul 09
1
version problems of rkward in ubuntu hardy repository
Hi
I tried to install rkward under ubuntu hardy heron, but it tried to
use the one from the cran repository which was newer, but it did not
install. To be able to install rkward, I had to disable the cran
repository, install rkward, lock it's version and enable the
repository again.
I have the feeling, that the dependencies are not correct (it is
looking for libqt4GUI, but I have a
2008 Sep 08
1
Vorticity and Divergence
Hi all,
I have some wind data (U and V components) and I would like to compute
Vorticity and Divergence of these fields. Is there any R function that
can easily do that?
Thanks in advance for any help
Igor Oliveira
CSAG, Dept. Environmental & Geographical Science,
University of Cape Town,
Private Bag X3,
Rondebosch 7701. Tel.: +27 (0)21 650 5774
South Africa Fax: +27 (0)21