similar to: lme: random effects of a quadratic term

Hi, I am searching for a way to solve the following problem: I want to use an external function, defined in a dyn.load()ed shared object, in another dyn.load()ed shared object. Currently I have to take the sources (Fortran) from one libraries src/ directory and copy them into the src/ dir of the other library, resulting in two copies of this function. This is bad for maintanance and maybe also

Hi Gustav, Try this: lapply(1:length(models),function(i) lapply(models[[i]],function(x) summary(x)$coef[2,]))[[1]] #1st list component [[1]] #??? Estimate?? Std. Error????? z value???? Pr(>|z|) # pm10 #5.999185e-04 1.486195e-04 4.036606e+00 5.423004e-05 #[[2]] #??? Estimate?? Std. Error????? z value???? Pr(>|z|) #ozone #0.0010117294 0.0003792739 2.6675428048 0.0076408155 #[[3]] #???

Hello Luke and Gabriel, Thank you very much for your quick responses. The explanation of STDVEC is very helpful and I appreciate it! For the wrapper, I have a few new questions. 1. Like Luke said a mutable object is not possible. However, I noticed that there is one extra argument *deep* in the function duplicate. I've googled all the available documentation for ALTREP but I did not find

Hello, I have encountered two bugs when using ALTREP APIs. 1. STDVEC_DATAPTR >From RInternal.h file it has a comment: /* ALTREP support */ > void *(STDVEC_DATAPTR)(SEXP x); However, this comment might not be true, the easiest way to verify it is to define a C++ function: void C_testFunc(SEXP a) > { > STDVEC_DATAPTR(a); > } and call it in R via > a=1:10 > >

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))

Hi - I am trying to find and describe missing data in a time series. For instance, in the library openair, there is a data frame called "mydata": library(openair) head(mydata) date ws wd nox no2 o3 pm10 so2 co pm25 1 1998-01-01 00:00:00 0.60 280 285 39 1 29 4.7225 3.3725 NA 2 1998-01-01 01:00:00 2.16 230 NA NA NA 37 NA NA NA 3 1998-01-01 02:00:00

I have a set of evaluation variables (n) for each sample (sample size is large enough) and I am trying to use R (nnet package) to aggregate the data. However, I don't know the weight for each variable (I am sure the weight shouldn't be equally assigned). Specifically, I have 12 indices (CO2, SO2, TSP...) for 100 cities and I want to calculate weight for each and finally obtain a

Hi, I would like to compare the hazard functions of two samples using the Cox proportional hazards model. For sample 1 I have individual time-to- event data. For sample 2 I don't have individual data, but grouped data that allows to obtain a hazard function. I am wondering if there is an R function that allows to obtain a hazard ratio of the two hazard funtions (under the

I will choose to use the winbind. Based on the link that Rowland said: https://wiki.samba.org /index.php/Setup_Samba_as_an_AD_Domain_Member I followed the steps as described in the tutorial. I created symlinks. In the main DC I added this line in smb.conf: idmap_ldb: use RFC2307 = yes Changed /etc/nsswitch.conf passwd: files winbind shadow: files group: files winbind hosts:

Jiefei, Inline. On Thu, May 16, 2019 at 2:30 PM ??? <szwjf08 at gmail.com> wrote: > Hello Luke and Gabriel, > > Thank you very much for your quick responses. The explanation of STDVEC is > very helpful and I appreciate it! For the wrapper, I have a few new > questions. > > > 1. Like Luke said a mutable object is not possible. However, I noticed > that there is

Hi Jiefei, Thanks for tryingout the ALTREP stuff and letting us know how it is going. That said I don't think either of these are bugs, per se, but rather a misunderstanding of the API. Details inline. On Thu, May 16, 2019 at 11:57 AM ??? <szwjf08 at gmail.com> wrote: > Hello, > > I have encountered two bugs when using ALTREP APIs. > > 1. STDVEC_DATAPTR > > From

Hi R users,I am a beginner in the use of R. I need urgent help for my thesis study. <https://stats.stackexchange.com/posts/633206/timeline> I have daily air pollution parameters PM10, PM2.5 CO, NO2, SO2, and O3. I also have daily hospital admission numbers. Taking into account the effect of weekends and holidays, I would like to used generalised additive model (GAM) to explore the

Dear R HELP, ABOUT glmmPQL and the anova command. Here is an example of a repeated-measures ANOVA focussing on the way starling masses vary according to (i) roost situation and (ii) time (two time points only). library(nlme);library(MASS)

Hi, I would like to plot the following equation: dF(x)/dx=(k1+k2F(x))(1-F(x)) where k1 and k2 are parameters that I have estimated already. How can I plot the curve in R? Thanks! Montserrat Rue Universitat de Lleida (Spain) [[alternative HTML version deleted]]

This bottle of Quoi D?Neuf by Herv? Villemade, a primeur from the Loire, was a real pleasure, a thirst wine that you sip without second thoughts and restriction to the last drop : This Pinot Noir with a bit of Gamay is fruity, gouleyant and fresh. This proves that even tasted several months after its release, a Vin Nouveau can deliver. It arrived on the market on the third thursday of november

I am trying to create a frequency distribution and I am a bit confused. Here are the commands I have entered: > data <- read.csv(file="40609_sortedfinal.csv",head=TRUE,sep=",") > NumberOfActionsByStatus = data$STATUS > NumberOfActionsByUser = data$ETS_LOGIN > NumberOfBidOffer = data$BID_OFFER > NumberOfActionsByUser.freq = table(NumberOfActionsByUser) >

The locate libnss_winbind.so display. /lib/libnss_winbind.so /lib/libnss_winbind.so.2 /lib/x86_64-linux-gnu/libnss_winbind.so /lib/x86_64-linux-gnu/libnss_winbind.so.2 /lib64/libnss_winbind.so /lib64/libnss_winbind.so2 /usr/lib/x86_64-linux-gnu/libnss_winbind.so /usr/lib/x86_64-linux-gnu/libnss_winbind.so.2 /usr/local/samba/lib/libnss_winbind.so /usr/local/samba/lib/libnss_winbind.so.2

I've rolled up R-0.99.0.tgz a moment ago. You can get it from ftp://cvs.r-project.org/pub/CRAN/src/base/R-0.99.0.tgz or http://cvs.r-project.org/pub/CRAN/src/base/R-0.99.0.tgz or wait for it to be mirrored at a CRAN site near you within a day or two. There's also a version split in three (!) for floppies if you prefer that. The next release is to be 1.0.0 at the end of this month,

I've rolled up R-0.99.0.tgz a moment ago. You can get it from ftp://cvs.r-project.org/pub/CRAN/src/base/R-0.99.0.tgz or http://cvs.r-project.org/pub/CRAN/src/base/R-0.99.0.tgz or wait for it to be mirrored at a CRAN site near you within a day or two. There's also a version split in three (!) for floppies if you prefer that. The next release is to be 1.0.0 at the end of this month,

Hello R-Users I am struggling with this line plot, it might be simple but I am missing something here. First of all I want to make multiple line plots across seasons (DJF,MAM,JJA,SON) for 12 variables (here, called nodes) and fill them with the node. So that season=x-axis, node=line col and freq=y-axis. My plot currently links all DJFs, MAMs, JJAs and SONs, however I will like them to be