similar to: The best solver for non-smooth functions?

Dear R helpers exposure <- data.frame(id = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20), ead = c(9483.686,50000,6843.4968,10509.37125,21297.8905,50000,706152.8354, 62670.5625, 687.801995,50641.4875,59227.125,43818.5778,52887.72534,601788.7937, 56813.14859,4012356.056,1419501.179,210853.4743,749961,6599.0862), pd =

Dear, I am using the solnp command (package Rsolnp) for a problem with equality and inequality constraints. I am not getting convergence for my problem but apart that I get 1 warning per iteration saying: ?In cbind(temp, funv) : number of rows of result is not a multiple of vector length (arg 1)?. I checked for equality and inequality functions and they seem fine to me. Where this

I have a dataframe called cvar, with two variables (among many others) called MSA and ACTUP. Both are numeric. This was working fine. Then I found out that for two MSAs, ACTUP should be 1, not 0. so I tried cvar$ACTUP[cvar$MSA == 6840] <- 1 cvar$ACTUP[cvar$MSA == 5360] <- 1 but when I try table(cvar$MSA, cvar$ACTUP) the level of ACTUP for those two MSAs has not changed, and is still

Hi List, Merlin Magix hardware v02 I'm trying to get asterisk to act as a voicemail server for a lucent merlin magix PBX that we purchased used. We have 4 FXO channels between the two PBXs on a Sangoma A200 card. The 770 dialgroup is working properly, in that calls to 770 are answered by Asterisk. The magix is sending mode codes in the format #XX#XXX#, where the 2nd block of digits

Dear R helpers, Thanks a lot for your earlier guidance esp. Mr Davind Winsemius Sir. However, there seems to be mis-communication from my end corresponding to my requirement. As I had mentioned in my earlier mail, I am dealing with a very large database of borrowers and I had given a part of it in my earlier mail as given below. For a given rating say "A", I needed to have the bad-wise

Hi I have a large credit portfolio (exceeding 50000 borrowers). For particular process I need to add up the exposures based on the bands. I am giving a small test data below. rating <- c("A", "AAA", "A", "BBB","AA","A","BB", "BBB", "AA", "AA", "AA", "A", "A",

-------------- next part -------------- An embedded message was scrubbed... From: Gildas Mazo <gildas.mazo at curie.fr> Subject: Re: [R] optimization subject to constraints Date: Tue, 10 Aug 2010 15:49:19 +0200 Size: 4924 URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100810/78862894/attachment.eml>

R-help, I'm trying to plot the following: year lgd 1 1986 136.97479 2 1987 69.10377 3 1988 67.66744 4 1989 71.60316 5 1990 62.06897 6 1992 6.25000 7 1993 27.72021 8 1995 23.83648 9 1996 10.29412 10 1997 95.67487 11 1998 82.09367 12 1999 56.60401 13 2000 29.80864 14 2001 23.77535 15 2002 48.30378 16 2003 83.47571 17 2004 74.58711 There are 2 missing years

When I try citypro <- sas.get('c:/ndri/cvar/data', member = 'citypro2') I get the following errors and warnings: 'sas' is not recognized as an internal or external command, operable program or batch file. 'less' is not recognized as an internal or external command, operable program or batch file. Error in sas.get("c:/ndri/cvar/data", member =

Hi, I am not sure how to fix the following error. LGD <- c(11.6, 12.3, 15.8, 33.1, 43.5, 51.3, 67.3, 84.9) cor (x=(file [1:47231,3:10]), y= rep (LGD, 47231), method = "pearson") Error in cor(x = (file[1:47231, 3:10]), y = rep(LGD, 47231), method = "pearson") : incompatible dimensions > sessionInfo() R version 2.13.0 (2011-04-13)

I have a dataset with 4 years of students, and normally I want to estimate things using each individual year, so I have a for loop as follows for (i in 1:4){} However, the only way I know how to calculate estimates using all four years of data is to put the estimations outside of the loop. Is there anyway to make a for loop that uses all four years at once, then uses each individual year?

Dear R helpers, I am using Beta distribution to generate the random no.s (recovery rates in my example). However, each time I need to save these random no.s in a csv format. To distinguish different csv files, one way I thought was use of Sys.time in the file name. My code is as follows - # My code rr = rbeta(25, 6.14, 8.12) lgd = 1 - mean(rr) write.csv(data.frame(recovery_rates = rr), file =

Dear R forum I have following variables - EAD = 10000 LGD = 0.45 PD = 0.47 M = 3 # Equation 1 R = 0.12*(1-exp(-50*PD))/(1-exp(-50)) + 0.24*(1-(1-exp(-50*PD))/(1-exp(-50))) b = (0.11852 - 0.05478 * log(PD))^2 K = (LGD * pnorm((1 - R)^(-0.5) * qnorm(PD) + (R / (1 - R))^0.5 * qnorm(0.999)) - PD * LGD) * (1 - 1.5 * b)^(-1) * (1 + (M - 2.5) * b) RWA = K * 12.5 * EAD > RWA [1] 22845.07 #

Hi, I'm currently beginning to use R for financial analysis (mainly Basel II benchmarks) and I would like to know if any R-User can give me some initial directions on packages and tutorials which I can use to calculate capital requirements, default probabilities, and related stuff. Thanks in advance, Rick [[alternative HTML version deleted]]

Hello, Is there any way we convert SAS dataset into R dataset? Kindest Regards, Peppy Adi-Purnomo ------ Peppy Adi-Purnomo Energy Market Analyst Energy Link Ltd Dunedin - New Zealand Ph.: +64 3 479 2475 Fax: +64 3 477 8424 Email: s.adi.purnomo at energylink.co.nz www.EnergyLink.co.nz -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello again I have a file from SPSS in .sav format. when I run library(foreign) cvar<-as.data.frame(read.spss("c:\\NDRI\\cvar\\data\\cvar2rev3.sav")) I get a warning Warning message: c:\NDRI\cvar\data\cvar2rev3.sav: Compression bias (0) is not the usual value of 100. The data appear to be OK, but I am concerned. (I tried searching the archives and the documenation for data

I used garchFit function to fit 1600 observations of EURO/USD 2-day returns in GARCH(1,1) model. As part of the summary I got warning message: NaNs produced in: sqrt(diag(fit$cvar)) And didn't get any estimates for 3 params' std.error, t value or probability: Error Analysis: Estimate Std. Error t value Pr(>|t|) mu -0.004827 0.020141 -0.240 0.811 ar1 0.010311

Hello, I just configured a computer on my local network to work as a PDC with Samba, although I'm not able to join the Domain from the windows computer I want to add. I though it was that the server was unreachable, but I can ping the server from the computer and viceversa. If someone please can help me, I'll appreciate it. The smb.conf file is here: ------------------- [global]

hey, i have this website (with ajax also) with lot of records of employee, and as an extra feature i want that the customer can import his employeelist (perhaps 100+ records), lot of work to do this manually. i would create a template like this Firstname Lastname Phone xxxx xxxxx xxxxx yyyy yyyyy yyyyy then they upload the excel file or save it as csv file

I recently had a problem with the configure script from perfparse running properly so that I could 'make install' it to work with Nagios 2.10 on latest CentOS 4.x + RPMForge (see http://readlist.com/lists/centos.org/centos/5/28966.html). Thinking it was a problem that was a result of "yum remove gd", I formatted the box and reinstalled everything. Now, oddly, I still have the