search for: semipar

Hi all, Does anyone know how to report heteroskedasticity and autocorrelation-consistent standard errors when using the "spm" command in SemiPar package? Suppose the original command is sp1<-spm(y~x1+x2+f(x3), random=~1,group=id) Any suggestion would be greatly appreciated. Thanks, Susan [[alternative HTML version deleted]]

...to be archived due to outstanding QC problems, but have CRAN and BioC packages depending on them which would be broken by the archival (and hence need archiving alongside). Therefore we are looking for new maintainers taking over maintainership for one or more of the following packages: R2HTML SemiPar cghseg hexbin lgtdl monreg muhaz operators pamr Package maintainers whose packages depend on one of these may be natural candidates to become new maintainers. Hence this messages is addressed to all these maintainers via BCC and to R-devel. See <http://CRAN.R-project.org/package=R2HTML&g...

...d this question yesterday and for some reason the post seems to be attached to the wrong thread. Also, I extended my test a little and it seems to indicate the problem is with spm. I would appreciate any help. Thanks. ========================================================== library(plyr) library(SemiPar) data <- data.frame(id=c(rep("111",100),rep("222",200)), value=c(rnorm(100,2,1),rnorm(200,10,5)), lhs=c(rnorm(100,2,1),rnorm(200,10,5))) #this works d_ply(data, c("id"), function(x) { print(lm("lhs~value", data=x)) }) #this works data111 <- data[...

I'm taking the residual-regression approach to semiparametric estimation (Robinson 1988, Econometrica), and basically using SemiPar simply as a convenient means of doing multivariate nonparamteric additive models. The final bit of code is here: finalfit <- spm(res~f(V3,basis="trunc.poly")+f(V5,basis="trunc.poly")+f(V6,basis=&...

Hi everyone: I want to fit the truncated polynomial smoothing to the quantiles instead of means, does someone know how to do it? I am thinking that maybe I can use SemiPar package, but can not find how. Thanks so many, Suyan

Hi there, I am reading the 2004 paper "Smoothing with mixed model software" in Journal of Statistical Software, by Ngo and Wand. I tried to run their first example in Section 2.1 using R but I had some problems. Here is the code: library(nlme) fossil <- read.table("fossil.dat",header=T) x <- fossil$age y <- 100000*fossil$strontium.ratio knots <-

...to know this, since assuming a Gaussian distribution will lead to substantial errors. It is a time series; the x values have equal five minute intervals. Here is the problem: I have an embarrassment of riches. Unless I am mistaken, the following R packages will do this: locfit, aws, sm, gss, semipar, pgam, gregmisc, to name just a few. Questions: Which package should I use, and which function in the package should I use for the regression? Here is the Google search result, restricted to r-project.org: http://www.google.com/search?as_q=locally+weighted+poisson+regression&num=10&hl=e...

...oordinates that form an irregular sampling area, is there an R function to determine boundary points (coordinates defining the limits of the area), either with or without user interaction ? # for example, given the following irregular sampling area, how could I define the boundary polygon? library(SemiPar) data(scallop) plot(scallop$longitude, scallop$latitude) Thanks. --Dale

Hi, When you use the default knot choice in spm() (library(SemiPar)) a figure showing the knot locations is sent to the screen and you have to accept the knots to move on. I am trying to run simulations using this function. Is there a way to get spm() to use the default "REML" knots without needing to approve each set of knots? Here is an example: l...

...------------------------------- ------- when choosing (with or without the tcltk window) "Bern 1 (swiss)" as CRAN mirror to install a package (function install.packages) for R-2.5.0, the list of available packages for osx seems not complete: AMORE, AdaptFit, AlgDesign, AnalyzeFMRI, SemiPar, cluster, nlme, survival ------------------------------------------------------------------------ ------- when choosing "Bern 2 (swiss)" as CRAN mirror to install a package for R-2.5.0, I get the following warning: > install.packages() --- Please select a CRAN mirror for use i...

...rco > sessionInfo() R version 2.2.1, 2005-12-20, i386-pc-mingw32 attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils" [6] "datasets" "base" other attached packages: mvtnorm SemiPar cluster lme4 lattice Matrix "0.7-2" "1.0-1" "1.10.4" "0.995-2" "0.12-11" "0.995-5" nlme MASS "3.1-68.1" "7.2-24" >

Dear R-Help members, I am using Simon Wood`s mgcv package version1.7-22and R version 3.0.0 (2013-04-03) for fitting a GAM-Model to the LIDAR Data contained in the "SemiPar" package. Here is the code for fitting the model and for plotting the result: data("lidar") attach(lidar) ### # mgcv fitting ### gam_fit <- gam(logratio ~ s(range, k = 40, bs = "cr"), gamma = 1.4, method = "REML") plot(gam_fit, res = TRUE, shade = TRUE, pch...

...ons Zk I have transformed the basis to Z = Zk %*% penalty^(-1/2). Full details are in section 2: Bayesian analysis for penalised spline regression using winBUGS. https://www.jstatsoft.org/article/view/v014i14/v14i14.pdf I am specifically looking for an alternative to the function spm() in package semipar. Thank you for any help or advice. Dave. David Mcnulty Medical Statistician Tel: +44 (0) 121 371 2448 Internal: 12448 Email: David.Mcnulty at uhb.nhs.uk Web: http://www.uhb.nhs.uk Health Informatics - University Hospitals Birmingham NHS Foundation Trust Yardley Court, 11-13 Frede...

Hi All, I am trying to download semnet package but getting the error: package not available for R version 3.4.4. I tried downloading it from install.packages('semnet',repos='http://cran.us.r-project.org') and install.packages('semnet',repos='http://cran.revolutionanalytics.com/') and even the https://cran.r-project.org/src/contrib/semnet.tar.gz but doesnt seem to

I have got following error in drawing wavelet fitting. can some one help? > library(faraway) > data(lidar) > newlidar<-lidar[c(1:128),] > library(wavethresh) > wds <- wd(newlidar$logratio) > draw(wds) Error in plot.default(x = x, y = zwr, main = main, sub = sub, xlab = xlab, : formal argument "type" matched by multiple actual arguments [[alternative HTML

Dear R users, I'm new to both R and to this list and would like to get advice on how to build generalized additive models in R. Based on the description of gam, which I found on the R website, I specified the following model: model1<-gam(ST~s(MOWST1),family=binomial,data=strikes.S), in which ST is my binary response variable and MOWST1 is a categorical independent variable. I get the

...","sampling","samr","sandwich","SASmixed","sca","scape","scapeMCMC","scatterplot3d","scope","scuba","sde","seacarb","seas","seewave","sem","SemiPar","SensoMineR","SenSrivastava","seqinr","SeqKnn","seqmon","session","setRNG","sfsmisc","sgeostat","s hapefiles","shapes","SharedHT2","simex") TEMP <- c(TEMP,&qu...