similar to: Plot with residuals in mgcv

Hi, I am using the gam function in the mgcv package, I have random effects in my model (bs="re") this has worked fine, but after I updated the mgcv package to version 1.7-19 I recive an error message when I run the model. > fit1<-gam(IV~s(RUTE,bs="re")+s(T13)+s(H40)+factor(AAR)+s(V3)+s(G1)+s(H1)+s(V1)+factor(LEDD),data=data5,method="ML") > summary.gam(fit1)

Hi, Is there a function similar to excel's hlookup in R ? Thanks, Silje

Hi I have a small problem, I'm new in using R, so I hope you can help me... I'm running a logistic regression model and want to do a nice plot. The plot I have made is done a plot with the command jitter: plot(jitter(overto$age[overto$sex=="F"]),jitter(overto$neg.pos[overto$sex=="F"]),xlab="age",ylab="neg and pos") my responsvariable is positive

Hello, I have tried to find this out some other way, but unsuccessful I have to try this list. I assume this should be quite simple. I have a dataset with 4 columns, "Sample_no", "Species", "Nitrogen", "Carbon" in csv format. In the species column I have many different species with varying number of obs per species Eg "Sample_no"

Dear Rexperts, Given, aar <-function(command) { switch(command, {scrn = cat("scrn :Screening","\n")} {dx = cat("dx :Diagnosis","\n")} {df = cat("df :Don't Forget","\n")} ) } I want to be able to do: aar("dx") # function does cat("dx :Diagnosis","\n")

hi, I have a question regarding cp values in rpart(). When I use plotcp() I get a figure with cp values on the x-axsis, but then I use printcp() the cp values in that list are different from the values in the figure by plotcp(). Does someone know why? Silje [[alternative HTML version deleted]]

I have a question regarding NA in randomForest (in R). I have a dataset which include both numerical and non-numerical variables, and the data includes some NA. I tried to use na.roughfix but then i get an error message "na.roughfix only works for numeric or factor". I also tried rfImpute but this does not work either because I have some NA in my response variable. Does anyone have som

Dear All, I like to simulate a physiologically based pharmacokinetics model using R but am having a problem with the daspk routine. The same problem has been implemented in Berkeley madonna and Winbugs so that I know that it is working. However, with daspk it is not, and the numbers are everywhere! Please see the following and let me know if I am missing something... Thanks a lot in advance,

Hola, Estoy intentando realizar un ajuste a un modelo setar con R y me surgen problemas a la hora de representar el modelo setar sin la constante. El código es el siguiente: # Estimacion TAR(2;2,2) con delta 1: mod.setar <- setar(x, m = 2, mL = 2, mH = 2, thDelay = 1) mod.setar summary(mod.setar) mod.setarc<-setar(x, m=2, mL=2, mH=2, thDelay=1, include=c("none")) ##eliminamos

While running a multiple VM testscase with each VM running iperf traffic between others the following kernel NULL pointer exception was seen. Race appears when the tun driver instance of one VM calls skb_array_produce (from tun_net_xmit) and the the destined VM's skb_array_consume (from tun_ring_recv), which could run concurrently on another core. Due to which the sock_wfree gets called again

While running a multiple VM testscase with each VM running iperf traffic between others the following kernel NULL pointer exception was seen. Race appears when the tun driver instance of one VM calls skb_array_produce (from tun_net_xmit) and the the destined VM's skb_array_consume (from tun_ring_recv), which could run concurrently on another core. Due to which the sock_wfree gets called again

Hi all, I am trying to use LLVM's AliasAnalysis pass, but from a standalone tool that uses llvm libraries. The following is the code snippet I am currently using. PassBuilder PB; auto mod_manager = ModuleAnalysisManager { }; PB.registerModuleAnalyses(mod_manager); AAResults& AAR = mod_manager.getResult<AAManager>(*M); But the code fails at .getResult with the following error:

Hi there I just tried running Freetel with wine does anyone have this working? Freetel breaks with a message: "Cannot bind socket" i wonder if this has anything to do with using a non-TCP protocol it would be very nice to get this working 'cos there do not seem to be any good linux web phone apps that can talk to people with windows machines. wine gives

don't read this. i'm a complete newbie, and i don't know how to start sending messages. if someone read this, then can answer me also how to delete a file with the console (really ashamed) y know that copy with cp, but delete??? thanks... ___________________________________________________________ 1GB gratis, Antivirus y Antispam Correo Yahoo!, el mejor correo web del

Hello I'm using logistic regression from the Design library (lrm), then fastbw to undertake a backward selection and create a reduced model, before trying to make predictions against an independent set of data using predict.lrm with the reduced model. I wouldn't normally use this method, but I'm contrasting the results with an AIC/MMI approach. The script contains: # Determine full

Hi, I have an Adaptec eSATA card, model 1225SA: http://www.adaptec.com/en-US/products/Controllers/Hardware/sata/entry/AAR-1225SA/ Although the chipset is sata_sil24, it is not recognized by default by the kernel because its PCI ID is different than the ones recognized by this driver on the CentOS kernel. The driver recognizes PCI IDs such as 0x3132, 0x0242, 0x3131, but not 0x0244, which is the

Hi all, I am trying to follow the instructions from the Quick Start guide, but I am running into problems when issuing the following command: mkfs.xfs -i size=512 /dev/sdb1 /dev/sdb1: No such file or directory If I create the sdb1 folder in the /dev folder I get the following error when issuing the mkfs.xfs command mkfs.xfs: cannot open /dev/sdb1: Is a directory Any assistance would be most

TB --- 2013-10-31 19:50:43 - tinderbox 2.20 running on worker01.tb.des.no TB --- 2013-10-31 19:50:43 - FreeBSD worker01.tb.des.no 9.1-RELEASE-p4 FreeBSD 9.1-RELEASE-p4 #0: Mon Jun 17 11:42:37 UTC 2013 root at amd64-builder.daemonology.net:/usr/obj/usr/src/sys/GENERIC amd64 TB --- 2013-10-31 19:50:43 - starting RELENG_10 tinderbox run for i386/pc98 TB --- 2013-10-31 19:50:43 - cleaning the

Hello Hello i have this DF: > c CHR_NR diffdato1 x 1022 10395 1994-05-30 11.0 39 0 1994 1 1 0 0 24 1 1 24 29100 11377 2003-05-22 17.0 24 0 2003 1 1 0 0 15 1

I successfully joined my Redhat Linux 7.2 box to my Windows NT domain. I can ping it's IP address but not its Netbios name. I have it set to obtain the IP address form the NT DHCP. Also on the RH Linux box, when I try to browse the local network, I receive the error: Unable to connect to host localhost. My samba configurations are okay. And I have the nameserver set in the resolv.conf file