similar to: ANOVA question

Hi,everyone, I am studying the generalized additive model and employ the package 'mgcv' developed by professor Wood. However,I can not understand the example listed in check.in function. For example, library(mgcv) set.seed(1) dat <- gamSim(1,n=400,scale=2) ## fit a GAM with quite low `k' b<-gam(y~s(x0,k=6)+s(x1,k=6)+s(x2,k=6)+s(x3,k=6),data=dat) plot(b,pages=1,residuals=TRUE)

I have a list of measurements taken at 1 second samples in a vector (y). All I want to do is generate a regression through them, and obtain the residuals in another vector. I am currently using the 'line' function from the eda package to do this (ie : z <- line(time,y) rsd <- residuals(z) ) but this method does not handle NA values in the data. I would like to try to use the lm

Hi, I have a question concerning the analysis of some affymetrix chips. I downloaded some of the data from GEO GSE11324 (see below). In doing so I'm stuck after I identified the probesets with significant changes. I have problems in assigning probeset specific gene names as well as getting the genomic coordinates. Furthermore I have no clue how to deal with the fact, that most genes have

Hi, Tinc experts There’s a requirement, where A connect to B, B connect to C; where B knows all the subnets from A and C, but A wouldn’t need to know anything about C, and C wouldn’t need to know anything about A. Any tinc configuration can make this happen? You can think of B is not an interconnect/forwarding node, instead, it only connect to A and C separately. I can think of run two tinc

Hi, I'm trying to use Robust non-linear regression to fit dose response curves. Maybe I didnt look good enough, but I dind't find robust methods for NON linear regression implemented in R. A method that looked good to me but is unfortunately not (yet) implemented in R is described in http://www.graphpad.com/articles/RobustNonlinearRegression_files/frame.htm

Hello, 1.1pre17 does not support Windows XP.  The reason is that InitOnceExecuteOnce is not supported. I believe it is better to take the Windows XP name out of the download link. -rsd -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://www.tinc-vpn.org/pipermail/tinc/attachments/20190708/b68071ed/attachment.html> -------------- next part

Hello, A nice guy call Jun Shen was helping me out with this, but I require a bit more help. Below is my data set or list called 'test'. I'm trying to calculate the %RSD for each pair of index and keep it in cronological order if you can imagine a 3rd column with 'date' beside index. Result Index 1 0.2901 17 2 0.2928 17 3 0.2913 18 4 0.2893 18 5

I'm probably not using the function right, but anyway, the program said: "please bug.report." > w_R.f1=function(PR,Pc) return (Pc*PR^3) > > w_R.f2=function(PR,Pc) return (Pc*PR) > > w_R.fc=c(w_R.f1,w_R.f2) > w_R.fc[2] [[1]] function(PR,Pc) return (Pc*PR) > > RSD(type=3, + gparP=gpar(col="red"), + gparD=gpar(fill=FALSE,col=FALSE), +

I have a funny issue whereby when a application exits the proc:::exit probe doesn''t seem to be firing. A customer has written a dtrace script that outputs some info whenever proc:::exit fires but on occasion (random, no pattern in behaviour) dtrace exits without outputting the data. They are running dtrace via ''dtrace -s exit_trace.d -p <pid>'' Can anyone think

[this address is IPv6-only, and usually doesn't work, so don't.] There's another patch I have up my sleeve, that it would be nice if I didn't need to apply it by hand. Look: vorbis-tools/oggenc/encode.c --- encode.c-DIST Fri Jul 19 08:35:48 2002 +++ encode.c Fri Jan 9 07:46:27 2004 @@ -346,7 +342,7 @@ seconds = (int)(remain_time - (double)((int)remain_time/60)*60);

Hi all, Forecasting from an arima model is easy with predict. But I can't manage to backcast : invent data from the model before the begining of the sample. The theory is easy : take your parameters, reverse your data, forecast, and then reverse the forecast I've tried to adapt the predict function to do that (i'm not sure that the statistical procedure is fine (with the residuals),

Hello all, I use arima to fit the model with fit <- arima(y, order = c(1,0,1), xreg = list.indep, include.mean = TRUE) and would like to use predict() to forecast: chn.forecast <- rep(0,times=num.record) chn.forecast[1] <- y[1] for (j in 2:num.record){ indep <- c(aa=chn.forecast[j-1], list.indep[j,2:num.indep]) # this is the newxreg in the

Hello,guys: Recently, I am working on a seasonal ARIMA model. And I met some problem in the forecasting. Now I just want to know that How does R perform the predict procedure(the predict formula, the initial setting of errors,etc.)? I run the following commands and get the original code of the "predict" command, but I can't read it. Can anybody explain it to me? Thanks! saji from

Where can I find a list of supported multipath drivers for ZFS? Keith McAndrew Senior Systems Engineer Northern California SUN Microsystems - Data Management Group <mailto:Keith.McAndrew at SUN.com> Keith.McAndrew at SUN.com 916 715 8352 Cell CONFIDENTIALITY NOTICE The information contained in this transmission may contain privileged and confidential information of SUN

Dear all, I’m not sure if it is O.K. to ask this question here. But where can I find the code for the function ‘C_Cdqrls’ which is called by the R function ‘lsfit‘. Thank you all. Sorry for being naïve if so. -------------------- Ziqiang Zhao 2013-04-23 [[alternative HTML version deleted]]

When I run predict.Arima in my code, I get warnings like: Warning message: In cbind(intercept = rep(1, n), xreg) : number of rows of result is not a multiple of vector length (arg 1) I think this is because I'm not running predict.Arima in the same environment that I did the fit, so the data object used in the fit is no longer present. Looking at the predict.Arima source,

I'd like to write some numbers to an external file that looks "pretty" (e.g, decimal points aligned, same number of decimals). For example, if using sprintf(), "%5.1f" can be used to specify the format of the float to be printed on screen. How can I do the same if I want to write an output file instead? I have tried cat and write.table, but none of them worked so far.

At the moment I have 2 network connected in a VPN and I'm planning to extend this. Actual status: Network A has a public IP Network B is behind a provider-scale NAT. Actually I'm using openvpn with a single UDP connection from B to A, and hosts inside Network A can communicate with B. I tried with tinc, and I saw tinc has to go TCP-Only for accomplish this. This is because Tinc

Hi , I have a problem with NT 4.0 workstation with service pack 3 cant able to login to hp system. From FAQ what i understood is to modify the registry but even after modifying i couldn't able to login from NT 4.0 workstation. Can anybody can able to help me in this matter. Thx Ramesh -- ___________________________ Ramesh Ragineni Ilog (s) Pte Ltd Systems Engineer

Hello, I'm trying to create a for loop for a data set, I have a list of results in this data set and I want to take the 1st two add them together and divide by the mean of the 1st to, then do the same for the 3rd and 4th values in the list and so on and each time return a value for the calculation on each pair. Any help would be appreciated. Thanks, Al [[alternative HTML version