similar to: Problem with integrate function

Hi, I have R 2.6.0 with updated lme4 and Matrix packages, and I am trying to fit a nonlinear multilevel model. I get the following error message: Error in nlmer(f ~ grModel(x, w, Tmin, Tmax, Topt, kopt, m) ~ kopt | flat, : gradient attribute of evaluated model must be a numeric matrix and I wonder what this may indicate. The nonlinear model I try to fit is as follows: > grModel

hi! i send again my question because there was a problem earlier that someone did not see my attached file. If you really can't download it, this is the attached file. Please help me how to convert this S plus file to R. Is there a quick method to do it? I don't have an S plus installer here. --------------------------------------------------- # Computes a possible choice for

This is significant because it's a plugin which was from the last stable release containing undefined symbols (in the internal plugins). $ ldd -r tests/old-plugins/x86_64/Linux/v1.18.2/nbdkit-file-plugin.so linux-vdso.so.1 (0x00007ffe86cdd000) libgcc_s.so.1 => /lib64/libgcc_s.so.1 (0x00007f551b3e5000) libc.so.6 => /lib64/libc.so.6 (0x00007f551b21b000) /lib64/ld-linux-x86-64.so.2

Hello, I wuold like to propose my modifications of the original cor.test to you : I tried to calcolate the correct p-value for Spearman and Kendall's test with ties. Let me know what you think. Thanks you for your time. Antonietta di Salvatore test <- function(x, ...) UseMethod("test") test.default <- function(x, y, alternative = c("two.sided",

let''s try this with a shorter subject, which was originally this, which probably tripped the list-admin-approval triggers.... Xen 3.2 Debian 2.6.18 on Dual Xeon with Remote Serial access - grub.conf config advice needed the subject says it all - anyone who knows the answer can stop reading and hit "send" with a few lines of grub.conf - however, for those people who are bored,

Hello, I have written a function in order to analyse gaze paths. It works with the test data but when I try to apply the function to a data frame that stores "the real data" in columns I receive the error message that the " In if (pp > 1) { : condition has length > 1 only the first element will be used " I interpret this error message as saying that only the first

hi there ! the current version of tinc crashes with a segfault. it seems, as if has to to with fact, that i don`t have configured any subnets.... if the client tries to start up subnet_update() ist called for myself. the for iteration does not get an entry (because the subnet avl-tree is empty). the later call of net2str causes a segfault, because it gets an empty subnet and tries to access

Dear All, The function I wrote can run well with the small data, but with the large data, the function runs very very slowly. How can I correct it? Thank you very much. My function as below: a<-c(1:240) b<-c(1:240) l=function(a,b){ v=0 u=0 uv=0 v[1]=0 u[1]=0 uv[1]=0 for (i in 1:(length(s)-1)){ v[i]<-((gx[[i]][b,(gx[[i]][a,1]+1)])-(gx[[i]][a,gx[[i]][a,1]+1]))/(gx[[i]][a,gx[[i]][a,1]+1])

Hello, I am struggling to overcome following problem: I have matrix Vf and matrix V, which both have 3 columns and I want to create a spline function basing on coordinates from this matrices (more precisely coordinates of function nr 1: x are taken from Vf[,1] and y are taken from V[,1] respectively), because function apply() base on 1 argument I made it in this way:

于 2013/3/24 0:07, Joerg Sonnenberger 写道: > On Sat, Mar 23, 2013 at 04:59:49PM +0100, Tobias Grosser wrote: >> With Chris's reply you also got the general OK to upstream the >> patch, _AFTER_ the actual patches have been reviewed. > One thing I was asking on IRC for is whether it makes sense for new > backends to be committable in incremental steps. Especially for a >

Dear List: The square of the noncentral t-statistic with noncentrality parameter \delta is a noncentral F with noncentrality parameter \lambda=\delta^2. So, t^2_{\nu,\delta} = F_{1,\nu,\lambda=\delta^2}. Consequently, it should follow that t^2_{1-\alpha/2,\nu,\delta} = f_{1-alpha,1,\vu,\lambda=\delta^2}. However, this is not what is happening with the following code. The central

Hi. guys, I am trying to write my own Stochastic Gradient Ascent for logistic regression in R. But it seems that I am having convergence problem. Am I doing anything wrong, or just the data is off? Here is my code in R - lbw <- read.table("http://www.biostat.jhsph.edu/~ririzarr/Teaching/754/lbw.dat" , header=TRUE) attach(lbw) lbw[1:2,] low age lwt race smoke ptl ht ui ftv

On Sat, Mar 23, 2013 at 12:34 PM, Jiong Wang <jiwang at tilera.com> wrote: > 于 2013/3/24 0:07, Joerg Sonnenberger 写道: > > On Sat, Mar 23, 2013 at 04:59:49PM +0100, Tobias Grosser wrote: >> >>> With Chris's reply you also got the general OK to upstream the >>> patch, _AFTER_ the actual patches have been reviewed. >>> >> One thing I was

Hi, as Conrad suggested, I've made a complete list of all points in the specification, which I beleive are errors, or where the explanation is unclear, contains unneccessary steps and so on. I hope someone has time to look through the points and if and when accepting or rejecting them be so very kind and inform me about it. I will also once again try to work through the residue

Add an option to specify the number of MSI-X vectors for PCI NIC cards. This can also be used to disable MSI-X, for compatibility with old qemu. This option currently only affects virtio cards. Signed-off-by: Michael S. Tsirkin <mst at redhat.com> --- hw/virtio-net.c | 5 ++++- net.c | 18 +++++++++++++++++- net.h | 4 ++++ qemu-options.hx | 14 ++++++++------

Add an option to specify the number of MSI-X vectors for PCI NIC cards. This can also be used to disable MSI-X, for compatibility with old qemu. This option currently only affects virtio cards. Signed-off-by: Michael S. Tsirkin <mst at redhat.com> --- hw/virtio-net.c | 5 ++++- net.c | 18 +++++++++++++++++- net.h | 4 ++++ qemu-options.hx | 14 ++++++++------

Dear all, I have a big data frame: str(data1) 'data.frame': 18272 obs. of 11 variables: $ tag : int 100001 100002 100003 100005 100007 100008 100009 100011 100012 100014 ... $ sp : Factor w/ 18 levels "acassp","acocar",..: 13 5 7 14 14 18 3 11 13 10 ... $ gx : num 20 10 35 68 88 63 123 115 137 136 ... $ gy : num 30 25 24 1 10 40 45 25 23 45 ...

Hi, My apologies for asking this question, but could not find the answer elsewhere. I understand the glmnet package uses Fortran code. For example, the lognet.R file includes the lines of code shown below. But how can I see the Fortran code that is being referenced in the code below? Is that provided somewhere in the package source code? .Fortran("lognet",

Dear friends, Now, when i use the argument return(x=x,y=y,prob=prob) , R displays the waring message: Warning message: The return value for multiple variables wasn't used in: return(x = x, y = gy, prob = prob) I used the methods of "help.search("return")" and "?return" to get some help, but didn't find info on it. Anybody knows how it should be used

Hi all, Updated the patches for TILE-Gx backend: 1. added initial regression tests for tilegx codegen. 2. added initial regression tests for MC Layer. 3. fixed those commenting style issues. please review, thanks. I have tried to understand the new backend requirement for LLVM from the mailiing list archive, it's sure TILE-Gx backend will be actively maintained & improved, it's