Displaying 20 results from an estimated 1000 matches similar to: "recursive function help"
2012 Jul 27
3
bivariate normal
Dear list members
I need a function that calculates the bivariate normal distribution for each observation. It is part of a likelihood function and I have 1000's of cases. As I understand it I cannot use packages like "mvtnorm" because it requres a covariance matrix of the same dimension as the number of observations. Basically what I need is a function that takes as arguments a
2007 Jul 28
8
generating symmetric matrices
Greetings,
I have a seemingly simple task which I have not been able to solve today. I want to construct a symmetric matrix of arbtriray size w/o using loops. The following I thought would do it:
p <- 6
Rmat <- diag(p)
dat.cor <- rnorm(p*(p-1)/2)
Rmat[outer(1:p, 1:p, "<")] <- Rmat[outer(1:p, 1:p, ">")] <- dat.cor
However, the problem is that the matrix
2008 Nov 30
4
Sweex 1000VA UPS (Lakeview Research)
Hi all,
i've read a thread about this UPS from Peter van Valderen. He tried to
develop a
driver for this specific type of UPS. I've downloaded nut-2.2.2 and tried to
apply
the patches (lakeview.patch & lakeview.v2.patch) but both resulted in
rejected
chunks. I've tried to ./configure with type lakeview and then a make, but
the make
command fails.
Does anybody has any
2007 Jul 30
3
Constructing correlation matrices
Greetings,
I have a seemingly simple task which I have not been able to solve today and I checked all of the help archives on this and have been unable to find anything useful. I want to construct a symmetric matrix of arbtriray size w/o using loops. The following I thought would do it:
p <- 6
Rmat <- diag(p)
dat.cor <- rnorm(p*(p-1)/2)
Rmat[outer(1:p, 1:p, "<")] <-
2010 Apr 29
1
image function with date-time on X axis
I am trying to plot a image where the x axis has the units of time.
When I issue the
image(x,y,z) command with x as a POSIXct object, it fails to put a
time stamp on the
x axis.
Instead I get a warning "Incompatible methods" warning and no dates on
my x axis.
This example shows my problem:
Rmat=t(matrix(data=rnorm(1:500),ncol=10,nrow=50))
tax=seq(ISOdate(2010,4,14,12,0,0),
2002 Feb 27
1
Bug in glm.fit? (PR#1331)
G'day all,
I had a look at the GLM code of R (1.4.1) and I believe that there are
problems with the function "glm.fit" that may bite in rare
circumstances. Note, I have no data set with which I ran into
trouble. This report is solely based on having a look at the code.
Below I append a listing of the glm.fit function as produced by my
system. I have added line numbers so that I
2017 Dec 29
1
winbuilder warning message wrt function pointers
And remove the cast on the return value of R_GETCCallable. And check
that your function is found before using it.
#include <R.h>
#include <Rinternals.h>
#include <R_ext/Rdynload.h>
void bdsmatrix_prod4(int nrow, int nblock, int *bsize,
double *bmat, double *rmat,
int nfrail, double *y) {
DL_FUNC fun = NULL;
if (fun==NULL)
2010 Sep 15
1
optim with BFGS--what may lead to this, a strange thing happened
Dear R Users
on a self-written function for calculating maximum likelihood probability (plz
check function code at the bottom of this message), one value, wden, suddenly
jump to zero. detail info as following:
w[11]=2.14
lnw =2.37 2.90 3.76 ...
regw =1.96 1.77 1.82 ....
wden=0.182 0.178 0.179...
w[11]=2.14
lnw=2.37 2.90 3.76 ...
regw =1.96 1.77 1.82 ....
wden=0.182
2017 Dec 29
3
winbuilder warning message wrt function pointers
I've recently updated the coxme package, which calls internal routines from the bdsmatrix
package.? (It is in fact mentioned as an example of this in the Extensions manual.)
The call connections are a blocks like this, one for each of the 9 called C routines.
void bdsmatrix_prod4(int nrow,??? int nblock,?? int *bsize,
??????????????????? double *bmat, double *rmat,
??????????????????? int
2009 Nov 02
2
a prolem with constrOptim
Hi,
I apologize for the long message but the problem I encountered can't be stated in a few lines.
I am having some problems with the function constrOptim. My goal is to maximize the likelihood of product of K multinomials, each with four catagories under linear constraints on the parameter values. I have found that the function does not work for many data configurations.
#The likelihood
2009 Feb 10
3
summary of a list
Hello,
I'm using the following for loop to find regression curves using a list of functions (formList), a list of starting
values (startList), uppervalues (upperList) and lower values (lowerList).
A sample of the list of function I use in the loop is the following:
FormList <- list(PTG.P ~ fz1(Portata, a, b), PTG.P ~ fz2(Portata, a, b), PTG.P ~ fz3(Portata,a, b, d, e),
PTG.P ~
2009 Jun 20
1
png() resolution problem {was "Silhouette ..."}
Hallo Sebastian,
>>>>> "SP" == Sebastian P?lsterl <sebp at k-d-w.org>
>>>>> on Sun, 14 Jun 2009 14:04:52 +0200 writes:
SP> Hello Martin,
SP> I plotting the silhouette of a clustering and storing it as png. When I
SP> try to store the image as png the bars are missing. The bars are plotted
SP> when I use x11 or
2011 Oct 12
2
Nonlinear regression aborting due to error
Colleagues,
I am fitting an Emax model using nls. The code is:
START <- list(EMAX=INITEMAX, EFFECT=INITEFFECT, C50=INITC50)
CONTROL <- list(maxiter=1000, warnOnly=T)
#FORMULA <- as.formula(YVAR ~ EMAX - EFFECT * XVAR^GAMMA / (XVAR^GAMMA + C50^GAMMA)) ## alternate version of formula
FORMULA <- as.formula(YVAR ~ EMAX - EFFECT / (1 + (C50/XVAR)^GAMMA))
FIT <-
2004 Sep 04
5
R question
Hi,
Would you help me solve the following question? Thanks.
Question: If I try to set the probability=0.05 and find the approximate x. (The answer should be somewhere between 2.1782 and 2.1783.)
I write about this R program as follows but I don¡¦t know how to get the value of x which is between 2.1782 and 2.1783.
library(mvtnorm)
value<-array(1000)
a<-array(1000)
2010 Aug 13
2
Unable to retrieve residual sum of squares from nls output
Colleagues,
I am using "nls" successfully (2.11.1, OS X) but I am having difficulties retrieving part of the output - residual sum of squares. I have assigned the output to FIT:
> > FIT
> Nonlinear regression model
> model: NEWY ~ PMESOR + PAMPLITUDE * cos(2 * pi * (NEWX - POFFSET)/PERIOD)
> data: parent.frame()
> PMESOR PAMPLITUDE POFFSET
>
2001 Oct 06
0
calculating DNA mismatch distributions for large populations
Hi all,
I am interested in calculating and displaying the distributions of
pairwise comparisons between DNA sequences in populations. The
comparisons are the number of nucleotide sites that differ between the
two sequences (mismatches). My sequences are stored in a vector of
strings. There is an additional vector of the same length that
provides the indices to the DNA sequences. Finally, I
2018 Mar 22
1
Cannot install broom package
Hello,
I've problems installing several packages in my R on Fedora 27 64 bit. I
found out that it has to do something with a missing compiler
(libgfortran.so.3, see below).
It works if I downgrade the current version of libgfortran to the
specified version by downloading libgfortran-6.2.1-2.fc25.x86_64.rpm and
manually installing it.
However, I don't want to mess up my system,
2017 Oct 02
0
Issues with 'Miwa' algorithm in mvtnorm package
Good point. Now this returns 0.04062184. Hmmm.....
On Mon, Oct 2, 2017 at 6:30 PM, Hollie Johnson (PGR) <
h.a.johnson at newcastle.ac.uk> wrote:
> Hi Eric,
>
>
> Thanks for having a look into this. I think you have a small typo...
>
> B <- matrix(x, nrow=3, byrow = TRUE) should read B <- matrix(y, nrow=3,
> byrow = TRUE)
>
>
> Regards, Hollie
>
2017 Oct 02
5
Issues with 'Miwa' algorithm in mvtnorm package
Currently doing some work on local maxima on a random field and have encountered an issue with the Miwa algorithm used with the pmvnorm function in the mvtnorm R package.
Based on recommendations by Mi et al., we ran the mvtnorm package using the Miwa algorithm, since we have a maximum of 4 dimensions with non-singular matrices. However, running the estimation procedure in this way, we obtained
2006 Jan 23
1
mutlivariate normal and t distributions
Dear R-help list members,
I have created a package 'mnormt' with facilities for the multivariate
normal and t distributions. The core part is simply an interface to
Fortran routines by Alan Genz for computing the integral of two
densities over rectangular regions, using an adaptive integration
method. Other R functions compute densities and generate random
numbers.
The starting