similar to: plot probability density function (pdf)

Dear all, I have a lot of problems on R-programming. for example my csv. file is .. Date wrfRH wrfsolar wrfwindspeed wrfrain wrftd wrfta 21/10/2010 92.97 22.11 53.27 0 1546.337861 61.00852664 22/10/2010 87.35 21.99 40.89 0 1300.408288 62.85352227 23/10/2010 88.38 21.71 28.04 0.01 1381.768284 54.80594493 24/10/2010 92.32 15.45 22.38 0.51 1113.90981 39.46573663 25/10/2010 93.42

First of all, I have GAMs noxd<-gam(newNOX~pressure+maxtemp+s(avetemp,bs="cr")+s(mintemp,bs="cr")+s(RH,bs="cr")+s(solar,bs="cr")+s(windspeed,bs="cr")+s(transport,bs="cr"),family=gaussian (link=log),groupD,methods=REML) Then I type " summary(noxd)". and show Family: gaussian Link function: log Formula: newNO2 ~ pressure

Hi, I hope someone can help me with the following issue. I need find the minimum beta that satisfies the following: inf{beta>0 | f(x+beta*f(x))*f(x)<=0} where f() is a function and x is a sample statistic. Functions such as "nlminb" and "constrOptim" minimize a function and output the parameter (under parameter constraints). I need to minimize the parameter (also

Hi everyone! I have a simple question. my data is predict_SO2_a 1 39.793231 2 30.252578 3 32.467584 4 31.941509 5 27.908320 6 11.594137 7 9.368125 8 12.319093 9 11.558811 10 7.937192 11 11.211306 12 12.400342 13 12.393146 14 13.256160 15 10.709600 16 9.966334 17 28.850652 18 10.024405 I want

i read a book from WOOD, there's an example which is talking about the pollutant. library(gamair) library(mgcv) y<-gam(death~s(time,bs="cr",k=200)+s(pm10median,bs="cr")+s(so2median,bs="cr")+s(o3median,bs="cr")+s(tmpd,bs="cr"),data=chicago,family=Possion) lag.sum<-function(a,10,11) {n<-length(a) b<-rep(0,n-11) for(i in 0:(11-10))

I'm new R-programming user, I need to use gam function. y<-gam(a~s(b),family=gaussian(link=log),data) y<-gam(loga~s(b), family =gaussian (link=identity),data) why these two command results are different? I guess these two command results are same, but actally these two command results are different, Why? -- View this message in context:

I'm beginner in R! I have a lot of problems on R..... I have three questions about GAM 1. What is the function of Gaussian distribution in GAM?(if I choose family is Gaussian) Is it used in the predictand value (Y)? 2. How to plot a graph the gam function? For example: y<-gam(a~s(b),family=gaussian (link=log) ,Data) how to plot x axis is s(b) and y axis is log a??? 3. if I use GAM to

Hi again! I have a question about R. I have done gam in previous version of R with "mgcv" package and saved the workspace. This workspace contains different models and I will do prediction by these GAMs. However, I install new version of R. and use the same workspace. when I type summary(models), and the error message showed Error in Predict.matrix.cr.smooth(object, dk$data) : F is

Hola Buenas noches   Alguien me podria indicar en que paquete puedo encontras modelos lineales y cuadratico para ajustar una serie de datos   Gracias [[alternative HTML version deleted]]

Buenas tardes   Estoy intentando representar graficamente la elipse de confianza de dos pendiente en un modelo sin intercepto La formula que introduzo esplot(ellipse(recta,c(1,1)),type="l"), pero elipse me sale una recta   Gracias [[alternative HTML version deleted]]

buenas   Por favor me pueden decir como ajusto unos datos con las distribuciones poisson generalizada poisso-pascalgeneralizada conway-maxwell-poisson   no se en que paquette estan   Gracias [[alternative HTML version deleted]]

hi! I hope all of you can help me this question for example GAMs: ozonea<-gam(newozone~ pressure+maxtemp+s(avetemp,bs="cr")+s(ratio,bs="cr"),family=gaussian (link=log),groupA,methods=REML) formula(ozonea) newozone ~ pressure + maxtemp + s(avetemp, bs = "cr") + s(ratio,bs = "cr") #formula of gams coef(ozonea) # extract the coefficient of GAMs

I have 5 GAMs ( model1, model2, model3, model4 and model5) Before I use some data X(predictor -January to June data) to form a equation and calculate the expected value of Y (predictand -January to June). After variable selection, GAMs (Model 1)were bulit up! R-square :0.40 NOW, I want to use new X'( predictor -July - December data) and put into Model 1, then get the expected value of Y'

Dear all, I find some troubles about how to extact the row from csv. file by using if-statement condition. I want to extract the row if the rainfall is greater than the mean of rainfall and using the wrfta divided into 3 groups that's rainfall greater than mean -> group A ( create file group A_rain) -> groupB ( create file group B_rain)

2012/4/7 <r-help-es-request@r-project.org> > Envíe los mensajes para la lista R-help-es a > r-help-es@r-project.org > > Para subscribirse o anular su subscripción a través de la WEB > https://stat.ethz.ch/mailman/listinfo/r-help-es > > O por correo electrónico, enviando un mensaje con el texto "help" en > el asunto (subject) o en el cuerpo a:

>Date: Mon, 4 Aug 1997 10:05:04 +0200 >From: marc.waknine@espace.aerospatiale.fr >To: samba@arvidsjaur.anu.edu.au (Non Receipt Notification Requested) >Subject: printer problem . URGENT help form SAMBA team !! >Message-ID: <9708040805.AA06185@musun95009.espace.aerospatiale.fr> > >on the SAMBA server i declare some printers >everything is ok in win3.11 because all

Hi, there. Is there any function in R to plot the probability-probability plot (PP plot)? Suppose I am testing some data against normal. Thanks. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

I need to refresh my memory on Probability Theory, especially on conditional probability. In particular, I want to solve the following two problems. Can somebody point me some good books on Probability Theory? Thank you! 1. Z=X+Y, where X and Y are independent random variables and their distributions are known. Now, I want to compute E(X | Z = z). 2.Suppose that I have $I \times J$ random number

Hi All, When carrying out hist(samples,breaks=50,probability=TRUE), the column values are considerably greater than 1, which seams very unreasonable. The plot is attached. I think the column value of the hist plot should correspond to x$counts/sum(x$counts) (x=hist(samples,breaks=50,probability=TRUE)). The size of data is a little bit larger, causing failure of uploading. If you need the

IHi all, have a dataset with rows as plots and environmental data as columns. I have predicted the values using the following ed.pred<-predict(lda.ed,ed) #lda.ed the model, ed the env. variables used for the prediction plots I am wanting to know the formula used by predict.lda for calculating the posterior probabilities. Can anyone point me in the right direction? Thanks