search for: so2

Hello, I am estimating the following model: so2.lme<-lme(so2~1+I(alcadakm^2)+dia,data=subjectes2,na.action=na.omit) And when I try to plot the random effects of the quadratic term with respect to a covariate (mam) I get an error: > so2.lmeRE<-ranef(so2.lme,augFrame=T) > plot(so2.lmeRE,form=I(alcadakm^2)~mam) Error in plot.ranef.lme...

...timate?? 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]] #??? Estimate?? Std. Error????? z value???? Pr(>|z|) #so2 #0.0026595441 0.0009352046 2.8438097399 0.0044577658 #pm10 do.call(rbind,lapply(1:length(models),function(i) do.call(rbind,lapply(models[[i]],function(x) summary(x)$coef[row.names(summary(x)$coef)=="pm10",c(1:2,4)])))) #???????? Estimate?? Std. Error???? Pr(>|z|) #[1,] 0.0005999185 0...

...ortran): **** so1.c ******************* #include <R.h> void func1(double* a, int* n){ printf("in func1\n"); F77_SYMBOL(func2)(a,n); printf("back in func1\n"); } ****************************** gcc -g -fpic -c so1.c -I/usr/local/lib/R/include ld -shared so1.o -o so1.so **** so2.f ******************* subroutine func2(a,n) real*8 a(*) integer n integer i call intpr("in func2",8,0,0) do 10 i=1,n a(i)=2.0d0 * a(i) 10 continue return end ****************************** g77 -g -fpic -c so2....

...x); > } > return(x); > } If calling R_altrep_data1 return the internal data directly, we will only see one message. following my last example > .Internal(inspect(so1)) > @0x0000000005e7fbb0 14 REALSXP g0c0 [MARK,NAM(7)] Share object of type > double > > .Internal(inspect(so2)) > @0x0000000005fc5ac0 14 REALSXP g0c0 [MARK,NAM(7)] wrapper > [srt=-2147483648,no_na=0] > @0x0000000005e7fbb0 14 REALSXP g0c0 [MARK,NAM(7)] Share object of type > double > > sm1=peekSharedMemory(so1) > getting data 1 > > sm2=peekSharedMemory(so2) > getting data 1...

...e is the source code for the wrapper: https://github.com/wch/r-source/blob/trunk/src/main/altclasses.c#L1399 Here is a working example if one can build the sharedObject package from https://github.com/Jiefei-Wang/sharedObject n=100 > x=runif(n) > so1=sharedObject(x,copyOnWrite = FALSE) > so2=so1 > so2[1]=10 > .Internal(inspect(so1)) > .Internal(inspect(so2)) Here is my session info: R version 3.6.0 alpha (2019-04-08 r76348) > Platform: x86_64-w64-mingw32/x64 (64-bit) > Running under: Windows >= 8 x64 (build 9200) > Matrix products: default > locale: > [1]...

...> > If calling R_altrep_data1 return the internal data directly, we will only > see one message. following my last example > > > .Internal(inspect(so1)) >> @0x0000000005e7fbb0 14 REALSXP g0c0 [MARK,NAM(7)] Share object of type >> double >> > .Internal(inspect(so2)) >> @0x0000000005fc5ac0 14 REALSXP g0c0 [MARK,NAM(7)] wrapper >> [srt=-2147483648,no_na=0] >> @0x0000000005e7fbb0 14 REALSXP g0c0 [MARK,NAM(7)] Share object of type >> double >> > sm1=peekSharedMemory(so1) >> getting data 1 >> > sm2=peekSharedMem...

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)) b<-b+a[(i+1):(n-11+i)] b} death<-chicago$death[4:5114] time<-chicago$time[4:5114] o3<-lag...

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 2.76 190 NA NA 3 34 6.8300 9.6025 NA 4 1998-01-01 03:00:00 2.16 170 493 52 3 35 7.6625 10.2175 NA 5 1998-01-01 04:00...

...https://github.com/wch/r-source/blob/trunk/src/main/altclasses.c#L1399 > > Here is a working example if one can build the sharedObject package from > https://github.com/Jiefei-Wang/sharedObject > > n=100 > > x=runif(n) > > so1=sharedObject(x,copyOnWrite = FALSE) > > so2=so1 > > so2[1]=10 > > .Internal(inspect(so1)) > > .Internal(inspect(so2)) > > > Here is my session info: > > R version 3.6.0 alpha (2019-04-08 r76348) > > Platform: x86_64-w64-mingw32/x64 (64-bit) > > Running under: Windows >= 8 x64 (build 9200) >...

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 relationship between daily patients admissions, and air pollution parameters. I would like tu use mgcv package. How to...

...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 comprehensive index and rankings for the cities. Could anyone give me an idea how can I realize this? Thank you very much. -- View this message in context: http://r.789695.n4.nabble.com/Use-R-for-data-aggreg...

...p winbind ldconfig: Can not stat / libx32: not directory or file found ldconfig: Path / usr / lib 'Given more than once ldconfig: Path / usr / lib64 'Given more than once ldconfig: Can not stat / usr / libx32: not directory or file found libnss_winbind.so.2 -> libnss_winbind.so2 Could someone give me any tips on how to solve this problem?

...was vinified on a base of Cheverny red, whole-clusters covered with CO2 in a tronconic wooden vat. Usually Herv? Villemade makes the primeur with a majority of Gamay but this vintage was so beautiful and plentiful regarding the Pinot Noir that the blend was 90% Pinot Noir and 10% Gamay. This nearly-SO2-free and unfiltered wine had a price tag of about 6 ? at Le Verre Vol? (the wine shop at 38 rue Oberkampf in Paris) and Villemade made only 2000 bottles of it, so hurry up for the last bottles.... The label was designed by the Paris-based American artist Mattew Rose, who works with collages like th...

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

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,

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 /usr/src/samba4/bin/shared/libnss_winbind.so.2 > From: patocius at hotmail.com > To: rowlandpenny at googlemail.com; samb...

...ote: > > 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 > > /usr/src/samba4/bin/shared/libnss_winbind.so.2 > > > > > > &g...

...A NA NA ... $ PM10.diff : num NA NA NA NA NA NA NA NA NA NA ... $ O3.fresh : num NA NA NA NA NA NA NA NA NA NA ... $ CO2.fresh : num NA NA NA NA NA NA NA NA NA NA ... $ CH4.tunn : num NA NA NA NA NA NA NA NA NA NA ... $ NCH4.tunn : num NA NA NA NA NA NA NA NA NA NA ... $ SO2.tunn : num NA NA NA NA NA NA NA NA NA NA ... $ Druck.tunn : num NA NA NA NA NA NA NA NA NA NA ... $ Feucht.tunn : num NA NA NA NA NA NA NA NA NA NA ... $ Temp.tunn : num NA NA NA NA NA NA NA NA NA NA ... $ Luft.fresh : num NA NA NA NA NA NA NA NA NA NA ... $ Luft.waste : num...

On 22/03/15 21:55, Jhon P wrote: > 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 > /usr/src/samba4/bin/shared/libnss_winbind.so.2 > > > > From: patocius at hotmail.com > &g...