similar to: Can not restore domain from a shared state file

Hello, I've been working with "nnet" and now I'd like to use the weigths, from the fitted model, to iterpret some of variables impornatce. I used the following command: mts <- nnet(y=Y,x=X,size =4, rang = 0.1, decay = 5e-4, maxit = 5000,linout=TRUE) X is (m x n) Y is (m x 1) And then I get the coeficients by: Wts<-coef(mts) b->h1 i1->h1

I am trying to multiply a 3D array of 4x4x4 by the 4 3D arrays of a 4D array with dimensions 4x4x4x4 (the last dimension being the one that I want to split by). (4x4x4 array) > hiaAry , , a1 i1 i2 i3 i4 h1 9.5936098 6.001040 0.08772 0.3138600 h2 1.2003500 1.454570 2.79248 0.0000000 h3 0.1346500 0.201220 0.39256 0.5464000 h4 0.0109000 0.012270 0.16417 0.2766900 ,

Hello, I'm trying to calculate a certain distance estimator for my thesis. I have a program in C that works fine when I call it with .C() in R, but since I'm dealing with big matrices like 30000x20000 it was getting a stack overflow. Now I have the same program but more efficeintly coded using malloc, and it works perfectlry in C, compiles well with R CMD SHLIB but when I call it with

Hi, I'm having trouble within my function CalcPos to get it to call CalcHorz with values from each row. I *think* it's calling CalcHorz with the final values of the inputs and not the values from each row. How can I do this properly in R? The values aa,bb,cc,dd are inputs. CalcPos first calculates V1 and V2 vertically, and then I attempt to call CalcHorz to handle H1, H2 & H3

Hi, Am working on neural network. Below is the coding and the output > library (nnet) > uplift.nn<-nnet (PVU~ConsumerValue+Duration+PromoVolShare,y,size=3) # weights: 16 initial value 4068.052704 final value 3434.194253 converged > summary (uplift.nn) a 3-3-1 network with 16 weights options were - b->h1 i1->h1 i2->h1 i3->h1 16.64 6.62 149.93

Hi R-devel. Is the following behavior in g1() and h1() expected? It seems to make "..." arguments work slightly differently from named arguments. #missing() has the property that it looks "up the chain" #for example, "z" can be missing in f3 even if #that argument did have a name ("y") in f2 f1 <- function(x, ...) { cat("In f1, missing(x) is

Hi all, Here's my story. I have 3 hosts, 1 remote and 2 local. We'll say for the sake of clarity that h1 is remote and h2 is the local rsync pull host and h3 is also a local host on the net. Using rsync, I want to pull a data set first from h1 to both h2 and h3. The first pull copies all the data from h1 to h2 and h3. OK, that gives me a baseline of the data on both h2 and h3.

Full_Name: Gerald Guglielmo Version: 2.8.1 (2008-12-22) OS: OSX Leopard Submission from: (NULL) (131.225.103.35) When I attempt to use the hist() function to bin integers the behavior seems very odd as the bin boundary seems inconsistent across the various bins. For some bins the upper boundary includes the next integer value, while in others it does not. If I add 0.1 to every value, then the

Hi folks, here is a short RFC regarding hashing for the implementation of switch statements and my preliminary patch. I posted this patch on 2014-01-16 21:58 at llvm-commits at cs.uiuc.edu. You can find a copy e.g. on http://lists.cs.uiuc.edu/pipermail/llvm-commits/Week-of-Mon-20140113/201782.html. Best regards Jasper === Preliminary: Special identifiers n number of given switch labels

Hi Dear all, I have some gene expression data samples from different tissue types ----------------------------------------------- - 120 samples from blood (B) - 20 samples from Liver (L) - 15 samples from Kidney (K) - 6 samples from heart (H) ----------------------------------------------- All the samples are from different individuals, so there are in total 161 individuals from which the DNA was

Hi there. Thanks for your time in advance. I am using R 2.2.0 and OS: Windows XP. My final goal is to calculate 1/2*integral of (f1(x)^1/2-f2(x)^(1/2))^2dx (Latex codes: $\frac{1}{2}\int^{{\infty}}_{\infty} (\sqrt{f_1(x)}-\sqrt{f_2(x)})^2dx $.) where f1(x) and f2(x) are two marginal densities. My problem: I have the following R codes using "adapt" package. Although "adapt"

hi everyone! i've got a problem in optimizing the following function fun <- function(theta){ theta <- rbind(beta,lambda) fun <- sum(exp(h(beta)%*%lambda)) fun } where h(beta) is h <- function(beta,...){ g1 <- matrix(0,b,M) g2 <- matrix(0,b,M) h1 <- matrix(0,b,1) h2 <- matrix(0,b,1) for(f in 1:M){ for(F in 1:b){ g1[F,f] <-

Hi all, I would like create a new column in a data.frame (a1) to store 0, 1 data converted from a factor as below. a1$h2<-NULL for (i in 1:dim(a1)[1]) { if (a1$h1[i]=="H") a1$h2[i]<-1 else a1$h2[i]<-0 } My question: is it possible to remove the loop from above code to achieve the desired result? Thanks in advance, Jin Geoscience Australia Disclaimer: This e-mail

Hi, I am computing the Nelson-Aalen (NA) estimate of baseline cumulative hazard in two different ways using the "survival" package. I am expecting that they should be identical. However, they are not. Their difference is a monotonically increasing with time. This difference is probably not large to make any impact in the application, but is annoyingly non-trivial for me to just

Hi R users I am trying to run a non linear parameter optimization using the function optim() and I have problems regarding memory allocation. My data are in a dataframe with 9 columns. There are 656100 rows. >head(org_results) comb.id p H1 H2 Range Rep no.steps dist aver.hab.amount 1 1 0.1 0 0 1 100 0 0.2528321

#I have a data set that looks like this. A bit more complicated actually with # three factor levels but these calculations need to be done on one factor at a #I then have a set of different rates that are applied #to it. #dataset cata <- c( 1,1,6,1,1,2) catb <- c( 1,2,3,4,5,6) doga <- c(3,5,3,6,4, 0) data1 <- data.frame(cata, catb, doga) rm(cata,catb,doga) data1 # start rates #

Hi I managed to compile (and use) C code from R, but now I have to use C++ code. After reading "Writing R extensions", I came up with the following code, which compiles fine ------------------- // emdL1_R.cc: #include <time.h> #include <stdio.h> #include "emdL1.h" extern "C" { void emd_L1( double *H1, int *n1, double *H2, int *n2,

Hello! I''m everywhere searching the answer for my question. How to get parameters to controller from view using <input tag. But I want to do this without db. I have the code: class TranslateController < ApplicationController def start @time = Time.now end def result @rez = params[:name] end end <html> <head> <title> Translate numbers from 10 to

Hello. I'm trying to find a way in R to calculate between and within subject variances and confidence intervals for some analytical method development data. I've found a reference to a method in Burdick, R. K. & Graybill, F. A. 1992, Confidence Intervals on variance components, CRC Press. This example is for Balanced Data confidence interval calculation from Pg 62. The data are

Hello, I have just subscribed to this list. I will introduce myself: For some time, I have kept a markdown implementation in awk for personal use, different from other implementations. Now, I'm in the process of rewriting it and I'm trying to do it as compatible as possible. There are many questions I have, I know some test suites and am trying to pass those tests. When I don't know