similar to: iterative proportional fitting in R?

Hi, I am dealing with the following problem. There are two biochemical assays, say A and B, available for analyzing blood samples. Half the samples have been analyzed with A. Now, for some insurmountable logistic reasons, we have to use B to analyze the remaining samples. However, we can do a comparative study on a small number of samples where we can obtain concentrations using both A and

Cze??. Aktualnie mam do czynienia z klastrem pocztowym, w kt?rym pliki wiadomo?ci pocztowych u?ytkownik?w s? przechowywane na raid software-owym udost?pnianym po NFS. Serwer?w NFS jest kilka, ka?dy udost?pnia podzbi?r wiadomo?ci dla odpowiedniego serwera z postfix-em, kt?ry jest klientem NFS. Dane z jakiego serwera NFS host frontend z postfix-em powinien montowa? katalogi z poczt? s? w MySQL.

Hello I have usrmgr.exe (form NT server ver.5.1) to manage user accounts of samba.(one year ago) I set logon hours from 8:00 to 16:00. My server (Fedora 7) changed the time correctly from winter to sumer (+ 1 hour). I chenged hardware time clock -w (write system clock to hardware clock synchro). But users cant' login at 8:00 morning !!!! I must back logon hours one hour back in usrmgr

Hello I am using htb (Red Hat 9) and some of my clients have many drops (even up to 10% pkts) but there are no overlimits. What can cause this? ------------------------------------------------------------------- Chcesz pracowac w INTERIA.PL? Zobacz: >> http://link.interia.pl/f1855 << _______________________________________________ LARTC mailing list / LARTC@mailman.ds9a.nl

Hi, we have CentOS ftp server (vsftpd) which has a lot of users who are writing and reading a lot of small files from/into its own accounts (and other servers - using samba client - are reading these files and putting them into outside database). Since this server is under heavy load its availability is important. >From time to time we "crash" this server (don't ask why ...) but

Hi, I have problems with sieve fileinto and polish characters in folder names: # rule:[??????????????????] elsif anyof (header :contains "Subject" "??????????????????") { fileinto "INBOX/&AQUBBwEZAUIBRADzAVsBegF8AQQBBgEYAUEBQwDTAVoBeQF7-"; } dovecot log: Apr 5 15:17:42 mailer dovecot: lmtp(54203, test at test1.com): rmHVJHabfU+70wAA9hAjDg: sieve:

Hi, I am afraid I have a question to ask of you all. I have just completed setting up a mail server running on Ubuntu Server 10.04. It has postfix, dovecot 1.2.19, LDAP and squirrelmail as the webmail. I have also created virtual users accounts on the system through LDAP. I can send and receive mails which is great. Now, what I am trying to do is to set up a vacation auto-reply on the

Currently if user is 1MB under quota and someone tries to deliver mail that is over 1MB, Dovecot rejects the mail. But smaller mails aren't rejected probably for days. So user might not even realize that they didn't receive one of the mails. Also having a user "almost over quota" is a rather strange state I think. So what do you think about v2.2 allowing delivery of one last

Hi All, I've gotten a number of reports from users about gplots::heatmap.2 generating 'node stack overflow' errors. As it turns out, these errors originate in stats:::plotNode, and are triggered when it is passed a dendrogram that is too deeply nested. While increasing the stack size (which is a compile-time option for byte coded functions) can allow a particular instance to

Dear all, I wonder if there are some function or package in R which can do the iterative proportional fitting. In the exponential model f(y1,...yn)=exp(a'yi+b'(yi*yj)+.....+c'(y1*...*yn)+constance), instead of the canonical parameters, I use maginal probability instead of a and log odds ratio instead of b. and for the order higher than 3, I use the canonical way or even assume

I've been asked, and interested myself: Has anybody implemented the above in R or another S language dialect? We are talking about the algorithms / methodology by Wellner, Groeneboom and Jongbloed, e.g., from the following article @Article{Jongbloed:1998:ICM, author = "Geurt Jongbloed", title = "The Iterative Convex Minorant Algorithm for

I have a function that generated list and then from those I need to do some repeated calculation. So I was thinking to put them in another list and call them one after the other. Given that I have 20 of those it was a bit annoying to typing all of them in one command. So I need something to do a recursive addition of a list objet to a list here a small size example A<-list()

I can't for the life of me figure out how to get the roots for this simple (but sovable only iteratively) equation in R: x = z*(1-(1-2/z)^y where x and y are known, and z is unknown. In Matlab, this amounts to: [my.solution] = solve('x = z*(1-(1-2/z)^y') my.solution.real = solution(y-1,y) % bottom line displays non-imaginary solution (last element) Obviously, I'm deeply

Possibly of limited use to the original poster, but of interest more generally, there are a number of tools developed in the 70s for updating the matrix decompositions. There's at least one in my 1979 book "Compact numerical methods for computers" (still in print apparently) -- we didn't have enough space to keep all the data in the HP9830A machine. If there's sufficient

Hi R-helpers, Is it possible that I have the estimates from each step/iteration shown on the computer screen in order to monitor the process while I am using Optim or MaxLik? Thanks for your help. Maomao [[alternative HTML version deleted]]

Dear All, I am trying to iterate an iterative process (i know R is not the best place for so much looping but i have to so tough!) and store the resulting data from each iteration to then be graphed. In simple form the script looks like this for(i in tracks) { a <- (subset(data, data$Id==i)) result.dataset <- function(a, other.arguements) ### this function is the other

Dear all, I have a data frame of 2160 rows and 4320 columns, which I hope to condense to a smaller dataset by finding averages of 6 by 6 blocks of values (to produce a data frame of 360 rows by 720 columns). How would I go about finding the mean of a 6 x 6 block, then find the mean of the next adjacent 6 x 6 block, and so on, until the whole data frame has been covered? One slight twist is that

Hello, How can an interative unit root test be implemented in R? More specifically, given a time series, I wish to perform the Dickey Fuller Test on a daily basis for say the last 100 observations. It would be interative in the sense that this test would be repeated each day for the last 100 observations. Given the daily Dickey Fuller estimates of delta for the autoregressive process d(Y(t))

Hello Instead of reading in group1.txt I want to read in groups1 for the first iteration of i, then groups2 for the second and so on. Obviously I can't use groups(i) but assume there is a way to do this. group<-read.table("C:/Data/April 2005/group1.txt",header=T) thanks in advance Meredith

Dear all I am intrested in making iterative estimation (thro' loop statements) of, say, linear regression model. For this purpose, I have written the following programme and that I have made use of a sample data (viz., exp.txt): ? Programme: ? # Linear regression modelling with sample data (try5.txt) # Repeated estimation through loop statement