similar to: Bayesian estimates for the 1st-order Spatial Autoregressive model

We are hiring an R programmer to make biologists better at crop improvement: http://jobs.monsanto.com/missouri/research-and-development/jobid6130734-r-programmer-jobs Please apply if you are passionate about building R culture and infrastructure. . . .. . ... . ... . .. . . ... . ... . . . .. .. .... .. . ... . Barrett Foat, PhD Genome Data Analytics Team Lead

dear R experts: ?I am writing my own little newey-west standard error function, with heteroskedasticity and arbitrary x period autocorrelation corrections. ?including my function in this post here may help others searching for something similar. it is working quite well, except on occasion, it complains that Error in solve.default(crossprod(x.na.omitted, x.na.omitted)) : system is

thank you, achim. I will try chol2inv. sandwich is a very nice package, but let me make some short suggestions. I am not a good econometrician, so I do not know what prewhitening is, and the vignette did not explain it. "?coeftest" did not work after I loaded the library. automatic bandwidth selection can be a good thing, but is not always. as to my own little function, I like the

Hello, for my diploma thesis I need to program a solver for Merton?s respectively Black?s and Scholes? Option pricing formula, which should be achieved for several dates. What I want to do is to estimate the value of a firm?s assets "vA" (x[2] denotes vA) and the option-implied volatility of firm?s assets "sigA" (x[1] denotes sigA) by solving it simultaneous using the Black

predict.warn() -- a function to display factor levels in new data for linear model prediction that do not exist in the estimating data. Date: 2005-8-24 From: John C. Nash (with thanks to Uwe Ligges for suggestions) nashjc at uottawa.ca Motivation: In computing predictions from a linear model using factors, it is possible to introduce new factor levels. This was encountered on a practical

Hi, I don't know about the lazyeval package or what you are trying to do but to answer the main question "How to create columns dynamically using a function?" I would do something like that: # dataset dem <- structure(list(id = c("L1", "L2", "L3", "M1", "M2", "M3"), TEST_SET_NAME = c("A", "A",

Your code contains the lines stopifnot(!(any(data1$norm_sd >= 1))) if (!(any(data1$norm_sd >= 1))) { df1 <- dat1 return(df1) } stop() "throws an error", causing the current function and all functions in the call stack to abort and return nothing. It does not mean to stop now and return a result. Does the function give

Hi, I am writing a function to dynamically create column names and fill those columns with some basic calculations. My function "demo_fn" takes argument "blup_datacut" and I like to use the contents of those arguments to dynamically create new columns in my dataset. Please note that I have another function called "calc_gg" within the function "demo_fn".

Eric: I will try and see if I can figure out the issue by debugging as you suggested. I don?t know why my code after stopifnot is not getting executed where I like the code to run the funlp2 function when the if statement is TRUE but when it is false, I like it to keep running until the stopifnot condition is met. When the stopifnot condition is met, I like to get the output from if statement

Hi! All, I was thinking to use R to handle my large scale data of 3 Mi records On Unix. Could you let me know if there is limits of R for its matrices and vectors and what is the maximum allowed size if there is ? Thank you advance Yuefu St Louis --------------------------------------------------------------------------------------------------------- This e-mail message may contain

Hi, I have coding question on mixed model in R. I am using R2.11.0 in windows. I have an experiment with 2 fixed effect factors - A and B. The levels of B are within the levels of A factor. The model is very similar to a split plot design except the nesting relationship between the 2 fixed effect factors. For example: there are 2 levels for A - GM and ZM. There are 7 levels of B in total

I've a problem with trusting domain authentication , I've searched the newsgroup but found no answer . The configuration is very simple . A windows 2000 domain (DOMA) , which contains the samba server , and a Windows 2003 domain (DOMB). The two are trusted by a bidirectional trust . The problem is that users from the 2003 domain can't get authenticated by the samba server and are

Eric: Thanks for taking time to look into my problem. Despite of making the change you suggested, I am still getting the same error. I am wondering if the logic I am using in the stopifnot and if functions is a problem. I like the recursive function to stop whenever the norm_sd column has zero values that are above or equal to 1. Below is the calclp function after the changes you suggested.

The message is coming from your stopifnot() condition being met. On Thu, Dec 14, 2017 at 5:31 PM, DIGHE, NILESH [AG/2362] < nilesh.dighe at monsanto.com> wrote: > Hi, I accidently left out few lines of code from the calclp function. > Updated function is pasted below. > > I am still getting the same error ?Error: !(any(data1$norm_sd >= 1)) is > not TRUE? > > >

Hi, I accidently left out few lines of code from the calclp function. Updated function is pasted below. I am still getting the same error ?Error: !(any(data1$norm_sd >= 1)) is not TRUE? I would appreciate any help. Nilesh dput(calclp) function (dataset) { dat1 <- funlp1(dataset) recursive_funlp <- function(dataset = dat1, func = funlp2) { dat2 <- dataset %>%

Hi, The help file on princomp says that 'The signs of the columns of the loadings and scores are arbitrary, and so may differ between differnt programs for PCA, and even between different builds of R'. Does anyone know 1. Whether this applys to all functions like svd, eigen ... etc? 2. During a leave-one-out cross step, I notice the loadings also change sige when a different

If you are trying to understand why the "stopifnot" condition is met you can replace it by something like: if ( any(dat2$norm_sd >= 1) ) browser() This will put you in a debugging session where you can examine your variables, e.g. > dat$norm_sd HTH, Eric On Thu, Dec 14, 2017 at 5:33 PM, Eric Berger <ericjberger at gmail.com> wrote: > The message is coming from

I received the following e-mail in response to an e-mail I had sent to SSH communications questioning the wisdom of their requesting OpenSSH to change it's name. Contained in the message is that statement that SSH Communications did not exert their trademark rights earlier becuase it's only recently that OpenSSH has become more visible. In the United States, this would invalidate the

samba samba ? how are you? After a serious accident, I am ready to go back to work. At monsanto they let the solutions dry out in these silly racks. I don' work ther any more. Then I was working at the Danforth plant sience center. tThe whole building was down toone autoclave so I was there early But I dropped a 4l erlenmeyer and slipped and fell on the broken glass. They had to take

ARE YOU LOOKING FOR A WORK-AT-HOME JOB? -- EITHER PART-TIME OR FULL-TIME?? DO WANT ANOTHER PRODUCT TO COMPLEMENT YOUR EXISTING MLM PROGRAM?? OR... ARE YOU LOOKING TO REPLACE AN MLM THAT YOU ARE NOT TOTALLY PLEASED WITH?? If the answers to all of the above questions are no, then you can delete this message now. If the answer to any of the above questions is YES, then PLEASE READ ON!