similar to: [R-lme] Extract estimated variances from output of lme?

Does anybody know of a function that implements the derivative (gradient) of the multivariate normal density with respect to the *parameters*? It?s easy enough to implement myself, but I?d like to avoid reinventing the wheel (with some bugs) if possible. Here?s a simple example of the result I?d like, using numerical differentiation: library(mvtnorm) library(numDeriv) f=function(pars, xx, yy)

dear list, as a part my problem. I have to estimate some parameters using ML estimation. The form of the likelihood function is not straight forward and I had to use a for loop to define the function. I used "optim" to maximise the result but was not sure of the programme. To validate my results, I tried to write a function to obtain the MLE of a bivariate normal in the same manner. On

Dear R-help, I'm using the 'optim' functions to minimise functions, and have read the documentation, but I'm still not sure how to determine sensible values to use for the 'fnscale' and 'parscale' options. If I have understood everything correctly, 'fnscale' should be used to scale the objective function, so that for example if the default is

R-users E-mail: r-help@r-project.org Hi! R-users. I am just wondering what the definition of "dffits" in R language is. Let me show you an simple example. function() { library(MASS) xx <- c(1,2,3,4,5) yy <- c(1,3,4,2,4) data1 <- data.frame(x=xx, y=yy) lm.out <- lm(y~., data=data1, x=T) lev1 <- lm.influence(lm.out)$hat sig1 <-

Hi, I??m trying to use the Win2BUGS package from R and I have a similar problem that reurns with the message: Error in FUN(X[[1]], ...) : symbol print-name too long But, there is no stray ` character in the file ( Sugestions given by: Duncan Temple Lang <duncan> Date: Mon, 26 Sep 2005 07:31:08 -0700 ) The progam in R is: library(R2WinBUGS) library(rbugs) dat <-

Dear All- I am trying to merge two data files - they have different date formats and different times zones. I need to match up the date/time of the datasets and then invoke a conditional statement, such as: if dataC$mph is >= 12 then keep dataM$co23 for the corresponding time/date stamp. snippets of data files: *dataC.txt* LST in mph Deg DegF DegF2 % volts Deg

Dear knowing people, Dennis Murphy helped me a lot with my first loop last week. Thanks again - I could have made more than 10 "Thank-You cakes" in the time it saved me! But now I want to complicate the thing. My ideas didn't work. Let's see if anyone is smarter ;-) The following packages are needed: library(adehabitatHR) library(rgdal) library(plyr) # My dataframe looks

Dear all, I'm trying to get a Brownian bridge kernel (kernelbb) for each combination of two consecutive animal locations (see commands below) and put them, with a loop, inside a list. It works well at the beginning but after 42 runs, it appears the following warning : >Error in seq.default(yli[1], yli[2], by = diff(xg[1:2])) :  >  invalid (to - from)/by in seq(.) I looked at the

On Tue, Feb 09, 2016 at 05:53:55PM +0300, Roman Kagan wrote: > +and s_nic_model = Source_rtl8139 | Source_e1000 | Source_virtio_net [...] > +and source_video = Source_Cirrus | Source_QXL As you know there could be a lot of other input devices. At the moment the code basically ignores these (but it prints a warning "unknown [video|network] adapter model %s ignored"). I think we

libswfdec/swfdec_bits.c | 27 +- libswfdec/swfdec_swf_decoder.c | 46 ---- test/Makefile.am | 19 + test/swfedit.c | 81 +++++++ test/swfedit_file.c | 209 ++++++++++++++++++ test/swfedit_file.h | 59 +++++ test/swfedit_tag.c | 61 +++++ test/swfedit_tag.h | 56 ++++ test/swfedit_token.c |

Hello everybody, I need an help because I don´t know if the command for the ANOVA analysis I am performing in R is correct. Indeed using the function aov I get the following error:"In aov (......) Error() model is singular" The structure of my table is the following: subject, stimulus, condition, sex, response Example: subject stimulus condition sex response

Hi @all, I've created a Distributed-Replicated Volume consisting of 4 bricks on 2 servers. # gluster volume create glusterfs replica 2 transport tcp \ gluster0{0..1}:/srv/gluster/exp0 gluster0{0..1}:/srv/gluster/exp1 Now I have the following very nice replication schema: +-------------+ +-------------+ | gluster00 | | gluster01 | +-------------+ +-------------+ | exp0 | exp1 |

Package: xen-utils-4.0 Version: 4.0.1-1 Severity: wishlist [resending from another address since my ISP seems to be queueing emails for up to 25 hours and counting...] I did some experiments with the xen remus HA system. With overlord3:~$ debdiff xen_4.0.1-1.dsc xen_4.0.1-1lindi1.dsc diff -Nru xen-4.0.1/debian/rules.real xen-4.0.1/debian/rules.real --- xen-4.0.1/debian/rules.real 2010-08-02

Package: xen-utils-4.0 Version: 4.0.1-1 Severity: wishlist I did some experiments with the xen remus HA system. With overlord3:~$ debdiff xen_4.0.1-1.dsc xen_4.0.1-1lindi1.dsc diff -Nru xen-4.0.1/debian/rules.real xen-4.0.1/debian/rules.real --- xen-4.0.1/debian/rules.real 2010-08-02 16:10:13.000000000 +0300 +++ xen-4.0.1/debian/rules.real 2011-01-05 00:32:41.000000000 +0200 @@ -159,6 +159,7 @@

On 02/05/2014 11:26, Adam Strzelecki wrote: >> It should be fixed. Can you confirm? > > Nope. Sorry. Still doesn't work here: > > Get:3 http://llvm.org/apt/trusty/ llvm-toolchain-trusty/main libllvm3.5 amd64 1:3.5~svn207822-1~exp1 [7,300 kB] > Get:4 http://llvm.org/apt/trusty/ llvm-toolchain-trusty/main llvm-3.5-runtime amd64 1:3.5~svn207822-1~exp1 [52.0 kB] > Get:5

This was working just a few hours ago, and now my CI server gives: The following packages have unmet dependencies: clang-6.0 : Depends: libllvm6.0 (>= 1:6.0~svn298832-1~) but it is not going to be installed Depends: libstdc++6 (>= 4.9) but 4.8.4-2ubuntu1~14.04.4 is to be installed Depends: libstdc++-4.9-dev but it is not installable Depends:

below is a code to compute hessian matrix , which i need to generate 29 number of different matrices for example first element in x1 and x2 is use to generate let say matrix (M1) and second element in x1 and x2 give matrix (M2) upto  matrix (M29) corresponding to the total number of observations and b1 and b2 are constant.  can some one guide me or help to implement this please. I did not

Hi, I want to create a data frame similar to the following, but greatly scaled up: df <- data.frame(aaa= c("a","b","c"), integer(3), integer(3)) names(df)[2:3] <- paste("var",1:2,sep="") which yields aaa var1 var2 1 a 0 0 2 b 0 0 3 c 0 0 I would not relish having to paste 'integer(3)' 5000 times :(

Commit 47b5f245bec908f803f0a89c3b1e3166cfe33aad originally introduced the filtering of files by extension. For the QXL driver at least, there is a qxldd.dll file which is part of the driver, so that must be copied to the guest as well. This patch will cause 'WdfCoInstaller*.dll' files to be copied too. While I'm not sure if this is correct, it seems as if it will be safe as although

Hi, Using indicator variables I have been able to fit and run the code for fitting a bivariate mixed model using unstructured covariance matrix The code is lme.fit1<- lme(one.var~-1+indic1+indic2+I(indic1*d.time)+I(indic2*d.time), random =~ -1+indic1+indic2|m.unit, weights = varIdent(~1|indic1) ,data = new.data) My variables are one.var :- the two response variables stacked one after