similar to: Build Design Matrix with avoiding loops

-------------- next part -------------- An embedded message was scrubbed... From: Gildas Mazo <gildas.mazo at curie.fr> Subject: Re: [R] optimization subject to constraints Date: Tue, 10 Aug 2010 15:49:19 +0200 Size: 4924 URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100810/78862894/attachment.eml>

Hi, I'm getting crazy: This does work: library(Matrix) a1<-b1<-c(1,2) c1<-rnorm(2) aDgt<-spMatrix(ncol=3,nrow=3,i=a1,j=b1,x=c1) png("myImage.png") image(aDgt) dev.off() But this doesn't !!! f<-function(x){ png("myImage.png") image(x) dev.off() } f(aDgt) My image is saved as a text file and contains nothing at all !!! Thanks in advance, Gildas Mazo

Dear R users, I want to aggregate data in the following way: ### X <- data.frame(u = c("T1","T1","T1","T2"), v=c("a","a","b","a")) X library(sqldf) sqlOut <- sqldf("select count(distinct(v)) from X group by u") sqlOut ### Now I want to get the same result without using SQL. How can I achieve that ?

Dear R users, Given a user-defined cumulative distribution function F, I want to compute F^{-1}(x). How is that possible with R? Best Regards, -- Gildas Mazo PhD student MISTIS team at INRIA Grenoble, France [[alternative HTML version deleted]]

Dear R users, I'm looking for tools to perform optimization subject to constraints, both linear and non-linear. I don't mind which algorithm may be used, my primary aim is to get something general and easy-to-use to study simples examples. Thanks for helping, Gildas

Dear R users, I'm using optim to optimize a pretty complicated function. This function takes the parameter vector "theta" and within its body I use instructions like sigma<-theta[a:b]; computations with sigma... out<-c() for (i in 1:d){ a<-theta[(3*d+i):c] out[i]<-evaluation of an expression involving 'a' (I use symbolic differentiation) } Unfortunately

Hello All, Want very much to learn how to plot patient drug timelines. Trouble is I need to figure out how to do this today. So not much time for me to struggle with it. Hoping someone can just help me out a bit. Below are some sample data and code that produces what I think is the beginning of a very nice graph. Need to alter the code to: 1. Get the lines for the drugs to appear on the

Hi, I have a data set of points which are represented by 3 variables x,y,z where x is the position of the point on the x-absciss and y on the y-absciss. Each of my points has a value z, which I want to be displayed as follows: the more z is high, the more the color on the map is dark. How can I achieve this ? Thanks for your help -- View this message in context:

Hi all, I have installed centos 5.1 x86_64. 2.6.18-53.1.14.el5 Everything seems to be working fine except that the rtc module is not working. cat /proc/interrupts shows: CPU0 CPU1 0: 129 0 XT-PIC-XT timer 1: 330797 0 XT-PIC-XT i8042 2: 0 0 XT-PIC-XT cascade 3: 4 0

Hi I have come across a problem, with my tdp410 and soekris board (basically pc on a chip amd geode cpu). I am using the box as a firewall/asterisk box. The problem occurs when I drop ppp and I get dead loop dectiotn going, I seem to lose interrupts and get lots of messages in syslog from wctdm24xx saying missed interrupt increasing latency its out lined here

Here are another three problems with logarithmic scales: 1) segments() does not work with logarithmic scales. I suggest to change lines 962-973 in "plot.c": for (i = 0; i < n; i++) { if (FINITE(xt(x0[i%nx0])) && FINITE(yt(y0[i%ny0])) && FINITE(xt(x1[i%nx1])) && FINITE(yt(y1[i%ny1]))) { GP->col = INTEGER(col)[i % ncol];

Hello I am using a Junghanns QuadBRI ISDN card - the module name is qozap. If I like at my interrupt assignment, qozap is sharing interrupt 10 with libata and uhci_hcd. I think libata is the IDE hard drive module and uhci_hcd is a USB module. linux:~ # modprobe qozap linux:~ # cat /proc/interrupts CPU0 0: 12634579 XT-PIC timer 1: 10 XT-PIC i8042

hi list. after my unsuccessfull experiences with mISDN i tried again to implement a zaphfc based solution. problem is: sound on calls via capi is stuttering/broken and therefore unuseable. my conf: - cel 1300, 256mb ram - avm b1 via capi connected to my outgoing ISDN - acer surf pci via zaphafc and crosslink cable with termination as internal bus - siemens gigaset 3035 as internal phone -

> Subject: Re: [Asterisk-Users] Booting error - Unable to specify channel 2: > No such device I just went through all this as well. The best thing to do IMHO is to try to find a way to manually assign IRQ in the BIOS. Also, and this is what I didn't see at first, some slots SHARE IRQ. Avoid this! If you are not using USB at all, turn it off in BIOS if possible. Tell the BIOS

I have some new PC hardware I'm using with the Intel H67 chipset. I created a bzImage I use for a PXE environment for a recovery mode, installing OS images, etc. The new mobo doesn't have a VGA port but instead DVI. Everytime I unplug the DVI cable it seems to cause "Disabling IRQ #11". After that the PC is basically unusable. The /proc/interrupts shows that hdd, usb, nic,

Hi All, I have an equipment loaded with 4 X100P (numbered 1-4)) and one T400P (numbered 5-8). Everything works fine except that I cannot use one of the FXS ports (number 5). If I configure zapata.conf to recognize it, the whole system voice quality suffers. I've tried already to switch PCI slots, with no results. Below is a snapshot of my /proc/interrupts, maybe this can shed some light on

Hi all, I am using a TDM400P card with 2 FXS and 2 FXO modules. Everything seems nice, but i'm not able to make calls nor to receive any. When I try to make a call, I keep receiven the "all circuits are busy now" message, and when I receive calls, asterisk doesn't seems to care (don't get anything on the CLI) I'm using Asterisk 1.2.17 and Zaptel 1.2.16 from Xorcom's

Excuse my ignorance here, but I am desperately trying to isolate the IRQ for my TE110P card (shown below as t1xxp) Ive gone into my bios and disabled all usb , parallel, serial and some other devices, those that I needed to keep, I have moved off of IRQ 10 and onto IRQ 5, but everytime I boot up, I get usb-uhci and ehci_hcd using IRQ 10 as well as my Digium card. Does anybody know what these are

Thanks for the reply Bill. I am aware of the interrupts problem. To solve it I have already disabled my serial ports freeing up interrupts 3 and 4 and these are allocated to the two cards. This was done 2 months ago and has not solved the problem. Is there any way that something can wake up every now and then and generate these two interrupts? My current /proc/interrupts is as follows:

Dear List, I have a logical data frame with NA's and a grouping factor, and I want to calculate the % TRUE per column and group. With an indexed database, result are mainly limited by printout time, but my R-solution below let's me wait (there are about 10* cases in the real data set). Any suggestions to speed this up? Yes, I could wait for the result in real life, but just curious if I