similar to: repeat a function

Hi Petr, Many thanks for your response. Basically I want to create a probability matrix to be used in a trinomial tree going forward. This is the reason why I thought to build the matrix around 0 would be much more efficient. I need to loop through because the probabilities will depend on my node and is not always the same per row (e.g. if N> jmax, jmax being defined in another function) I

Hi Well, I am not an expert in this field so I cannot comment your approach. I wanted only to point out that building matrix your way is like scratching your left ear with right hand, especially in R. What if you want increase size of your matrix? E.g. you use function ProbUP once for row "0" and than for rows different from jmax (if I correctly understand your code). Use of any

Hi Eric I did not see any answer and frankly speaking I cannot provide you with canned help. AFAIK if a function is defined within another function (which is your case) it cannot be called directly so it is necessary to define it in global environment. > fff <- function(x) { + myf <- function(a) a+2 + myf(x)^2} > > fff(5) [1] 49 > myf(5) Error in myf(5) : could not find

I had some fun this afternoon coding up a 'maze generator' in R. I thought I'd pass along the fruits of my labor for everyone's amusement. As written, every point is connected to every other point, so feel free to 'start' and 'finish' anywhere you like. Have fun! --Robert PS. Feel free to pass along suggestions or comments.

Hi I am trying to Run AsteriskNow 1.5 which uses Centos 5.3 Distribution with Readonly Root Filesytem. Right now if I shutdown the system abrupty or there is power failure linux kernel doesn't not boot "Kernel panic" and other short of issues.prb due to corrupted files. Thats why i want to make root filesystem as read-only and placing some directories /var and /etc/asterisk etc as

I have a matrix of data. I need to scan the matrix and find every sequence from maxima to maxima across a row. I can write a loop to do this easily. Problem is, I can't figure out how to store the results. Each result is a vector of widely varying lengths. Ideally I'd like a vector of these, i.e. a vector of vectors, so I can quickly iterate through them and compute correlation

Hi all. This is more of a stats question, I suppose. Let's say I have two separate simple regressions of weight on year from two different datasets. I want to combine the regressions so that I can come up with a single equation for the total weight regressed on year. In reality, there is missing data, so I can't just sum the data across datasets and come up with a regression on the

I am trying to get a GSM gateway (Alpha Tech GSM Gateway Blue Gate Dual Band Analoog FXO) working with Asterisk. I had a working FXO configuration to a analog port of a small home 1/4 ISDN pbx. I used this same configuration to connect a GSM Gateway that is supposed to be connected to the external(FXO) analog port of a pbx. I can get my configuration to dial the mobile number via the gateway, but

I'd like to expand the following data to perform a daily time series. It should cover from '2012-07-01' to '2012-10-06' with the values I actually have being the mean from one point measurement to another. Does anyone has a clue to perform this task. structure(list(Date.beg = structure(c(15635, 15617, 15615, 15610, 15609, 15605, 15604, 15601, 15593, 15593, 15586, 15581,

hi all, my appl works well but i get the following error suddenly for few pages. to resolve this, i had to restart the server everytime. then it works fine for some time, but after some time it appears again. this is happening contineously. i searched for solution on this, but didn''t get much. can anybody help in this plaese? any help will be greatly appreciated... i''ve pasted

Dear all, I realised a correspondence analysis with function cca() of vegan library. Just like in Okansen (2010) in the example of R help: library(vegan) data(varespec) data(varechem) vare.cca<-cca(varespec~ Al + P + K, varechem) With plot.cca() function I represented the species matrix in the next way: plot(vare.cca,display="species") Being similar to: plot((c(-2,2)),(c(-2,2)),

hi, how can we assign an xml element value to some variable in controller? e.g. if we have <order> <id>1</id> <name>abc</name> <location>xyz</location> </order> and if we want value of name i.e. ''abc'' to be assigned in a variable declared in a controller, say order_controller, like this- temp = "value of xml element

I have constructed the function mml2 (below) based on the likelihood function described in the minimal latex I have pasted below for anyone who wants to look at it. This function finds parameter estimates for a basic Rasch (IRT) model. Using the function without the gradient, using either nlminb or optim returns the correct parameter estimates and, in the case of optim, the correct standard

Hi, Our firm is looking to see if it is possible to wrap some of more advanced lib in R like modreg.dll and specifically Loess function in there for both regression and prediction into MSVC++. We will probably take it to either VB or Python for implementation issues. Given that locfit has many utitlities built in for data transformation, this would be a good thing. If that's not

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hello vorbis, uhmm... i've just installed the RC3 packages (for debian, but i'm pretty sure it's not going to matter what linux (or even *nix) you use) and the man pages for "oggenc" still says "release candidate 2" up the top of the man page. This should possibly be checked next time <grins>. The rest of

Is there any R package that implements a bifurcating autoregression, aka the BAR(n) model? I've been reading the Huggins and Staudte paper, "Variance Components Models for Dependent Cell Populations", from the Journal of the American Statistical Association, 1994. Shawn Garbett <shawn.p.garbett at vanderbilt.edu> Vanderbilt Cancer Biology 220 Pierce Ave, PRB 715AA

Hello there ? There were folks in the recent time running DC++ on wine. I'm also one of them, and I have a small problem, which turns out to be very annoying. The action buttons at the top of the main DC++ window as well as tabs at the bottom tend to show "tooltips" containing information about the last files to be hashed. However, there is an error somewhere that these tooltip

Hello! I have problems with using of "readBin" to read files, which was written in C with "fwrite". In the C-File there is the following Code: fwrite(MyitINI,sizeof(itINItype),1,outfile); where MyitINI is a structure of the following form typedef struct{ int KernelFileSave; /* Determined, if Systemmatrix saved or not.*/ char KernelFileName[200]; /* A-Matrix name

Dear all, I?m trying to get a output table for age and the summary of a and b, stratified by epo as follows using summary.formula h<-data.frame(a=sample(c("A",NA),100,replace=T),b=sample(c("B",NA),100,replace=T),age=rnorm(100,50,25),epo=sample(c("Y","N"),100,T)) library(rms)

glm(A~B+C+D+E+F,family = binomial(link = "logit"),data=tre,na.action=na.omit) Error in `contrasts<-`(`*tmp*`, value = "contr.treatment") : contrasts can be applied only to factors with 2 or more levels however, glm(A~B+C+D+E,family = binomial(link = "logit"),data=tre,na.action=na.omit) runs fine glm(A~B+C+D+F,family = binomial(link =