similar to: D_f_norm function

- Hi every body, i am using D_f_norm function from library stat in R i am using it to find the standard deviation of the output of a function 'calculate_Z0' that has one of it is inputs as a random variable 'W' , with a standard deviation 'Sigm3LER', when i change the mean of this random variable 'Width', i expected the standard deviation output from d_f_norm to

Hello, I want to make an 2D plot of .xyz data. So plot on x-axis and y-axis and use a color scale for the z. For every x-location along a cross section of the soil at several depths a resistivity must be displayed. This must result in a picture/graph which shows the resistivity for that cross section. The depths differ for each x-location. How can I do that? In fact it is an cross

Hi Everybody, My problem is that nlminb doesn't converge, in minimising a logLikelihood function, with 31*6 parameters(2 weibull parameters+29 regressors repeated 6 times). I use nlminb like this : res1<-nlminb(vect, V, lower=c(rep(0.01, 12), rep(0.01, 3), rep(-Inf, n-15)), upper=c(rep(Inf, 12), rep(0.99, 3), rep(Inf, n-15)), control = list(maxit=1000) ) and that's the result :

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Please update m4/ogg.m4 Namely, new version should support --with-ogg=yes or --with-ogg=no (this enables ogg stuff, but does not add -L to OGG_LIBS, that is crucial to keep .pc files clean). ogg.m4 from libogg-1.3.0 will be sufficient. -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.11 (MingW32) Comment: Using GnuPG with Thunderbird -

Hello: We have got 2 directors and three backend mailbox servers. I get the following error " Error: connect(ipc-proxy) failed: Permission denied". I have read the article https://www.dovecot.org/list/dovecot/2012-August/137349.html and we have done as signaled in the two director servers with no success, So, we have tried now to set it like this: service ipc {

Hi Everyone, Can anyone help me resolve a problem that i'm having with nlminb. The problem is that it stops after just one iteration and returns the same values as "start" ones. Thank you very much for your help. Sincerely. -- Kamel Gaanoun (+33) (0)6.76.04.65.77 [[alternative HTML version deleted]]

Hola, en la misma definici?n de la funci?n: # P15. Precipitation Seasonality(Coefficient of Variation) # the "1 +" is to avoid strange CVs for areas where mean rainfaill is < 1) p[,15] <- apply(prec+1, 1, cv) Un saludo, Jorge On Martes, 19 de Junio de 2018 13:07:27 Marcelino de la Cruz Rot escribi?: > Hola Jaume: > > Si miras el c?digo de biovars() ver?s que la

Hi, I want to display some contour labels. It works well, but the thickness of the label is to low. I tried the "labcex" command contour(x,y,z, labcex=2) But only the size is rising and not the thickness of labels. Only if i set the value to "labcex=10" the thickness is good, but the size is to big, so i can't see anything ;) Is there any command to rise the thickness

VM traffic is already limited by a throughput limit, but there is no control over the maximum packet per second (PPS). In DDOS attack the major issue is rather PPS than throughput. With provider offering more bandwidth to VMs, it becames easy to coordinate a massive attack using VMs. Example: 100Mbits ~ 200kpps using 64B packets. This patch provides a new option to limit VMs maximum packets per

Full_Name: Adrian Custer Version: 0.90.0 OS: Linux on Thinkpad (pentium) and desktop (K6) Submission from: (NULL) (128.32.251.234) When I create a boxplot, and then try to overlay a lowess fit or just the points, the points do not appear in the highest level and the lowess curve does not reach the highest level. However, if I add one to each of the models, the problem is solved. I tried this

Hi, if I use the xyplot (or levelplot) function (lattice library) with the option axs="i", I have the problem that the tickmarks lie a bit outside the "plot-box". Consider for example: library(lattice) x<-seq(0,1,by=0.01) y<-seq(0,1,by=0.01) xyplot(y~x,type="l",xlim=c(0,1),ylim=c(0,1),scales=list

Full_Name: Sanford Weisberg Version: 1.9.1 OS: Win XP Submission from: (NULL) (160.94.148.2) > wm <- read.table(url("http://www.stat.umn.edu/~sandy/wmdata0.txt"), header=TRUE) > mean(wm$Spd1) [1] 7.7773 > summary(wm$Spd1) Min. 1st Qu. Median Mean 3rd Qu. Max. 0.222 4.780 7.550 7.780 10.200 21.600 The mean of this variable DOES NOT ROUND to the value

My question relates to problems that I'm having matching lme and lmer examples in P&B. using Matix 0.995 In the Oxide example in p167-170 I can't get the level 2 coefficient estimates to match the fm1Oxide model in lme is data(Oxide,package="nlme") lme(Thickness~1,Oxide) which I translate in Lmer syntax to fm3Oxide<-lmer(Thickness~

Hello, I am trying to increase the thickness of the axis in plot() without reverting to the use of paint programs i see posts on that topic for the xyplot function but want to see if i can do it with plot() because i've already setup my graph script using that i thought i could use axis() function and specify lwd="thickness" or lwd.axis= but that does not work like it does for

Here's a simple example: x <- 1:5 plot(x,x^2) lines(x,x^2) points(x,x,cex=2) lines(x,x,lw=3) legend("topleft",legend=c("y=x^2","y=x"),pch=1,pt.cex=1:2,lw=c(1,3)) The thickness of the circles in the legend changes with lw. If you change the lw argument in legend to lw=c(1,1) then the thickness of the circles goes back to normal. How can I make the above

I have searched google high and low for answers to this...and I have gotten many examples, howto, etc...but they all seem to have a slightly different configuration, and therefore, slightly different problems. Unfortunately, not enough of them show the network layout, along with the configuration, so it's hard to tell why certain IP are being used, and were they are on the network. I have

> > Hi all, > I have a vector xm say: xm = c(1,2,3,4,5,5,5,6,6) > > I want to return a vector with the corresponding probabilities based on the > amount of times the numbers occurred. For example, I should get the > following vector for xm: > prob.xm = c(1/9, 1/9, 1/9, 1/9, 3/9, 3/9, 3/9, 2/9, 2/9) > Using prop.table gives: Usage (with table) > prob.xm <-

I've been out of town, and too hastily deleted some earlier messages from this thread.... > > > Dear friends. I would like to link an object file made in C. > > Dyn.load rather wants to link a whole library and I only wanted to > > link this small object file. I followed instructions given for S - > > by Clive Loader in his software for changepoint handling and >

Hello everybody, this is my first request about R so I am sorry if I send it to a bad mail or if I am not very clear. So my problem is about the use of rfImpute from randomForest package. I am interested in imputations of missing values and I read that randomForest can make it. So i write the following code : set.seed(100); library(mlbench) library(randomForest) data(BreastCancer)

Paul: I may have found the issue (which is similar to your conclusion). I checked using egsingle in the mlmRev package as these individuals are strictly nested in this case: library(mlmRev) library(nlme) fm1 <- lme(math ~ year, random=~1|schoolid/childid, egsingle) fm2 <- lmer(math ~ year +(1|schoolid:childid) + (1|schoolid), egsingle) Checking the summary of both models, the output is