search for: prm

Dear group, Here is my df : df <- structure(list(DESCRIPTION = c("PRM HGH GD ALUMINIUM USD 09/07/10 ", "PRM HGH GD ALUMINIUM USD 09/07/10 ", "PRIMARY NICKEL USD 04/06/10 " ), CREATED.DATE = structure(c(18361, 18361, 18325), class = "Date"), QUANITY = c(-1L, 1L, 1L), CLOSING.PRICE = c("2,415.90", "2,415.90&quo...

...parameters from the .C call are not being passed correctly to the function. As an example, I have attached a GDB run of the code. I set a breakpoint on entry to the function I am calling from R. What is bothering me (and probably causing the segmentation faults I am seeing) is that the parameter prm=as.double(c(3.2,1.1)) is not 3.2,1.1 IMMEDIATELY after the call to .C. I am sure I am missing something very basic. Thanks, Sean > sessionInfo() R version 2.2.0, 2005-08-11, powerpc-apple-darwin7.9.0 attached base packages: [1] "methods" "stats" "graphics"...

...o a discrete distribution (concretely frequencies from 0 to 200). However, the program crashes in the attempt. The problem seems to be that these vectors only have positive integer numbers (+ zero). When I add a random very small positive decimal to the non-decimal part everything works fine (files prm1 & prpmr1). I attach two of these vectors so you can run the following code. I've also thought to divide both vectors by a real constant such as pi. Do you think these two approaches are acceptable? ? setwd("") require(CvM2SL2Test) prm = scan('prm.txt') prpmr = scan('p...

Dear group, Here is my data.frame : df <- structure(list(DESCRIPTION = c("PRM HGH GD ALU", "PRM HGH GD ALU", "PRIMARY NICKEL", "PRIMARY NICKEL", "PRIMARY NICKEL", "PRIMARY NICKEL", "STANDARD LEAD ", "STANDARD LEAD ", "STANDARD LEAD ", "STANDARD LEAD ", "STANDARD LEAD "...

Hi, I am running centos x686 and upgraded from 4 to 5 when it came out now I want to try the same with centos 6 where can I get the centos i686 release rpm I used does not have the rpm I used any help wood be greatfull -------------- next part -------------- An HTML attachment was scrubbed... URL:

...it does the job. > a=as.data.frame(substr(lme$DESCRIPTION,1,14)) > colnames(a)=c("DESCRIPTION") > lme=as.data.frame(c(a,lme[,2:3])) > lme DESCRIPTION CLOSING.PRICE POSITION 1 PRIMARY NICKEL 25,755.7100 0 2 PRIMARY NICKEL 25,760.8600 0 3 PRM HGH GD ALU 2,415.9000 0 4 SPCL HIGH GRAD 2,420.1000 -1 5 SPCL HIGH GRAD 2,420.4100 -1 6 SPCL HIGH GRAD 2,420.7300 1 7 SPCL HIGH GRAD 2,421.0500 1 8 SPCL HIGH GRAD 2,388.4300 0 9 SPCL HIGH GRAD 2,389.0000 0 10 SPCL HIGH GRAD...

Is there a repository of Facter plugins that people have made? I''m working on a few that i wouldn''t mind sharing when complete.

...del dat as a mixture of Gammas then plot. ## modelGammas <- function(dat, which='BIC') { set.seed(939458) fmo <- stepFlexmix(dat ~ 1, k=1:3, model=FLXMRglm(family='Gamma')) mdl <- getModel(fmo, which=which) print(smry <- summary(mdl)) print(prm <- parameters(mdl)) plotGammaMixture(dat, prm['shape', ], prm['shape', ]*prm['coef.(Intercept)', ], smry at comptab[, 'prior']) } ## ## Works well for a single Gamma distribution. ## set.seed(78483) dat1 <- rgamma(6000,...

...input on this strange issue much appreciated.. Hope I haven't missed something obvious. Thanks a lot, Aditi (P.s.- I've tried converting columns one by one to factors this time, and that works. >P1L55<-factor(P1L55) >levels(P1L55) >[1] "0" "1" Code: prm<-read.table("P:\\..... .csv", header=T, ...sep=",", ...) prmdf<-data.frame(prm) prmdf[2:13]<-lapply(prmdf[2:13], factor) ## action performed, no error message ##I tried to pick random columns and check >levels(P1L55) >NULL >is.factor(P1L96) >FALSE -...

...f : Pose1 <- structure(list(DESCRIPTION = structure(c(1L, 2L, 3L, 4L, 5L, 8L), .Label = c(" SUGAR NO.11 May/10 ", "COTTON NO.2 May/10 ", "PLATINUM Jul/10 ", "ROBUSTA COFFEE (10) May/10 ", "WHEAT May/10 ", "PRIMARY NICKEL USD", "PRM HGH GD ALUMINIUM USD", "SPCL HIGH GRADE ZINC USD", "STANDARD LEAD USD"), class = "factor"), POSITION = c(5, 3, -1, 15, 4, 2), SETTLEMENT = structure(c(3L, 5L, 2L, 1L, 4L, 8L), .Label = c("1,353.0000", "1,739.4000", "16.5400", &quot...

...<-read.table("...) modeldf<-data.frame(model) modeldf[2:13<-lapply(modeldf[2:13],factor) colms<-(modeldf)[4:13] ## 10 markers only in this file se<-c(1:1000) for(f in colms) { print("Marker") { for( i in 1:1000) { print("perm no.") print(se[i]) { peg.no.prm<-sample(peg.no, length(peg.no)) try(fit5<-lmer(data=modeldf, peg.no.prm~1 + (1|family/f))) print(summary(fit5)) capture.output(fit5, file="testperm5.txt", append=T) }}} } The data files are at: >> >> >> <http://www.4shared.com/file/131980362/460bdafe/Testvcomp...

Dear List, I am trying to run a mixed model which, on the R console, prints output as follows: [1] "Marker" [1] "perm no." [1] NA Linear mixed model fit by REML Formula: peg.no.prm ~ 1 + (1 | family/f) Data: modeldf AIC BIC logLik deviance REMLdev 3119 3134 -1555 3112 3111 Random effects: Groups Name Variance Std.Dev. f:family (Intercept) 0.0 0.000 family (Intercept) 0.0 0.000 Residual 178.5 13.360 Number of obs: 388, g...

Hallo, I own a NETUPS 450 (Model on the sticker: PRM 450i; upsstats reports "UPS 450 VA FW -0024") which perfectly works with the upscode2 driver. Maybe you want to add it to your list, even if the model is >10 years old. I suppose this is true for the whole series. According manual the come in the 450VA, 700VA, 1000VA and 1500VA stren...

> Dear All > > I have data as follows. > > D T M L > 0.20 1 03 141 > 0.32 1 07 62 > 0.50 1 05 49 > 0.80 1 04 46 > 0.20 2 14 130 > 0.32 2 17 52 > 0.50 2 13 41 > 0.80 2 14 36 > 0.20 3 24 120 > 0.32

...ple right now, from what I can tell. I''m starting with the introduction, since, of course, that''s where everyone else starts, and it has not kept pace with how I talk about Puppet in presentations and such. In the course of rewriting it, I had to create some wiki pages for PRM and modules: http://reductivelabs.com/trac/puppet/wiki/PuppetRecipeManager http://reductivelabs.com/trac/puppet/wiki/PuppetModules The Modules documentation is pretty sparse, so please add anything you think is important. Also, in general, I just want to make this clear: I am always looking...

Dear R users, The package 'mfp' that fits fractional polynomial terms to predictors. Example: data(GBSG) f <- mfp(Surv(rfst, cens) ~ fp(age, df = 4, select = 0.05) + fp(prm, df = 4, select = 0.05), family = cox, data = GBSG) print(f) To describe the association between the original predictor, eg. age and risk for different values of age I can plot it the polynomials and fitted coefficients as: plot(0.407*I((age/100)^-2) + -4.96*I((age/100)^-0.5) ~ age, GBSG) But I...

...or embedding that metadata in the manifest ? At some point, there''s also going to be a need to include dependency information in the module metadata (module webserver needs module funky_yumrepos), though that would probably not affect puppet itself, just the tool downloading modules, like prm. * My previous writeup at [2] should probably die; I''ll delete the page in a few days David [1] http://reductivelabs.com/cgi-bin/puppet.cgi/wiki/ModuleOrganisation [2] http://reductivelabs.com/cgi-bin/puppet.cgi/wiki/ModuleOrganisation

...structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 13L, 11L, 12L), .Label = c("COFFEE C Sep/10", "COPPER Sep/10", "CORN Dec/10", "CRUDE OIL miNY Sep/10", "GOLD Aug/10", "HENRY HUB NATURAL GAS Sep/10", "PALLADIUM Sep/10", "PRM HGH GD ALUMINIUM USD", "SILVER Sep/10", "SOYBEANS Nov/10", "SPCL HIGH GRADE ZINC USD", "SUGAR NO.11 Oct/10", "WHEAT Sep/10"), class = "factor"), pl = c(-8.20000000000002, 2.30000000000001, -47.25, -0.230000000000004, -12.7000000000000...

..."FEEDER CATTLE", "GOLD", "HEATING OIL", "HENRY HUB NATURAL GAS", "LEAN HOG", "LIVE CATTLE", "NY HARBR GASOLINE BLDSTK", "ORANGE J.", "PALLADIUM", "PLATINUM", "PRIMARY NICKEL USD", "PRM HGH GD ALUMINIUM USD", "RANDOM LENGTH LUMBER", "ROBUSTA COFFEE (10)", "SILVER ", "SOYBEAN MEAL ", "SOYBEAN OIL", "SOYBEANS ", "SPCL HIGH GRADE ZINC USD", "STANDARD LEAD USD", "SUGAR WHITE", "WHEAT...

...or the different methods (NR, BFGS, etc). I have figured out how to use all of them except the maxBHHH function. This one is different from the others as it requires an observation level gradient. I am using the following syntax: maxBHHH(logLik,grad=nuGradient,finalHessian="BHHH",start=prm,iterlim=2) where logLik is the likelihood function and returns a vector of observation level likelihoods and nuGradient is a function that returns a matrix with each row corresponding to a single observation and the columns corresponding to the gradient values for each parameter (as is mentioned i...