similar to: nonparametetric bivariate regression

Hi, Can I turn my R code automatically into C? If yes, how can I do that? Do I benefit by doing it? (I know I can do this in matlab and it increases the speed when there is a lot of looping...) thanks a lot. x.w ------------------------ Ximing Wu Department of Economics University of Guelph Guelph, ON, N1G 2W1 Canada Tel: (519) 824-4120 ext. 53014 Fax: (519) 763-8497 email: xiwu at uoguelph.ca

Hi there, I want to generate a large amount of functions, say f=function(x,t) exp(-t[1]-t[2]*g_1(x)-t[3]*g_2(1+x)) where g_1(x) and g_2(x) are from a long list of moments, such as x, x^2, log(x), log(1+x) .. and so on. Any suggestions on how to do this efficiently? thanks a lot. x.w

Greetings I have been trying to use metaplot to create a forest plot of my meta-analysis results. My data is in the form of ORs, so I've been inputting the ORs as the point estimate (mn) and the SEs of the logORs as the SE (se). However, whilst the point estimates appear to be displayed properly, many of the confidence intervals go below zero (not possible for an OR, as far as I know). I

I have been building an R function to calculate the ***observed*** (as opposed to expected) Fisher information matrix for parameter estimates in a rather complicated setting. I thought I had it working, but I am getting a result which is not positive definite. (One negative eigenvalue. Out of 10.) Is it the case that the observed Fisher information must be positive definite --- thereby

Hello R users I recently learned how to use this command: lapply(datum,"[[","ColumnName") Unfortunately, I don't know how exactly it works, what it's called (in particular the "[[" part], and what other things you can do with it (retrieve multiple columns?). Given datum is a list of dataframes with the same column, but different number of rows,

Hello seasons R users, Is it possible to store a complete regression result into an array? I've already been able to save individual regression coefficients, but would like to store the whole regression results into different arrays through a loop. So that in under different quantiles regressions, I would be able to create a loop and store the full regression result each time into a

Hey All, just a quick question about printer upload/download stuff that's been confusing the hell out of me for a couple days... This might be a Samba, or a CUPS problem, so i've posted this on both newsgroups, but if anyone's had a similar problem I'd really appreciate any help! I'm running Samba 2.2.7 and CUPS 1.1.17 on a RedHat 7.3 System.. I'm trying to get the

We have two Samba 4.2.3 servers with FreeRadius to authenticate wireless users against active directory. Using DNS, sometimes both servers end up using the same domain controller to authenticate users. I would like to distribute the load to different DCs. Is there a way to manually point Samba to certain DCs? I tried the following configuration: security = ads password server =

Hi Sir/Madam, I'm a researcher in university of Guelph, Canada and now considering using R to do some data analysis. I'm wondering whether there is a library available in R that includes algorithms for "archetypal analysis"? This is a method quite similar to principal components analysis that is designed to find "archetypes" or "pure types" from

Hi all, apologies for seeking advice on a general stats question. I ve run normality tests using 8 different methods: - Lilliefors - Shapiro-Wilk - Robust Jarque Bera - Jarque Bera - Anderson-Darling - Pearson chi-square - Cramer-von Mises - Shapiro-Francia All show that the null hypothesis that the data come from a normal distro cannot be rejected. Great. However, I don't think it looks

hi there, anybody knows how to import data into R from splus directly? I have a huge amount of data in splus format, which will cost forever for me to transform the data back to ascii or something else and then read them into R. thanks a lot. ximing -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Dear all, I need to plot an cumulative distribution plot of a variable and then to fit a distribution to that, probably a weibull or lognormal. I have plotted the ecdf as > plot(ecdf(x)) but I haven't managed to fit the distribution. I have as well attached the data. I would appreciate if you could help me on that. Thank you. Kind regards Maria -------------- next part --------------

Hello, I tried the R-Sig-Mac list with this query, but had no takers. I hope that the following isn't too far off the mark for this list. Many thanks in advance if someone can help me out! Quick summary: I can't get RSPerl working on a PPC G5 with the pre-compiled binary for Mac OS X, but it does work if I compile R from source with the "--enable-R-shlib" flag. If I

The strange thing is that I cloned this server from another server, and the original server does not have the "wbinfo -u" problem. Both servers have the exact same configurations.. Thanks. ----- Original Message ----- From: "Peter Bulin" <bulin.peter at gmail.com> To: dxu at uoguelph.ca Cc: "samba" <samba at lists.samba.org> Sent: Wednesday,

Hi All I would like to compute the raw p-value from permutation tests and I found mt.sample.rawp() from the package multtest almost similar to what I want to do. But in the function definition: mt.sample.rawp(V,classlabel,test="t",side="abs",fixed.seed.sampling="y",B=10000,na=.mt.naNUM,nonpara="n") I would like to choose my own t-test which is designed

Hi all, Setup: I'm running 2 machines (amd64, 16GB) with FreeBSD 9-STABLE (Mar 14 so far) acting as NFS servers. They each serve 3 zpools (holding a single zfs, hourly snapshots). The zpools each are 3-way mirrors of ggate devices, each 2 TB, so 2 TB per zpool. Compression is "on" (to save bandwith to the backend, compressratio around 1.05 to 1.15), atime is off. There is no

What version are you using? I also had this problem with a security update of debian. I reported this to then and they fixed it. So if you are using debian, update to the latest version and see if this works Op 8 jun. 2016 10:13 p.m. schreef Dennis Xu <dxu at uoguelph.ca>: Hi, I am checking again if there are any other suggestions. The Samba server is joined to AD successfully. I can

Also when I use "net ads search '(objectCategory=group)' sAMAccountName", I can see all the groups. But wbinfo -g gives nothing. wired! ----- Original Message ----- From: "Dennis Xu" <dxu at uoguelph.ca> To: "Andreas Schamanek" <schamane at fam.tuwien.ac.at> Cc: "samba" <samba at lists.samba.org> Sent: Friday, April 29, 2016

-------- Original Message -------- Subject: [AGDG-LIST:428] Summer Course in Guelph Date: Fri, 22 Dec 2006 09:12:24 -0500 From: Larry Schaeffer <lrs at uoguelph.ca> Reply-To: lrs at uoguelph.ca To: Animal Geneticist's Discussion <agdg-list at colostate.edu> The Centre for Genetic Improvement of Livestock at the University of Guelph is pleased to announce a one week summer course

Full_Name: Ximing Zhao Version: R 1.7.1 OS: MacOs X Submission from: (NULL) (64.136.27.51) When I used function plot in R 1.7.1, I found two bugs. My whole function is this: ***************** drawexp<-function(X,x,c,b){ exp<-function28(X,x,c,b) if(max(exp)>=max(dexp(x,rate=1))){ plot(x,exp, type="l",main= paste("Exponential ,n=",length(X),"