similar to: retrieve most abundant species by sample unit

I am trying to use lm for a simple linear fit with weights. The results I get from IDL (which I am more familiar with) seem correct and intuitive, but the "lm" function in R gives outputs that seem strange to me. Unweighted case: > x<-1:4 > y<-(1:4)^2 > summary(lm(y~x)) Call: lm(formula = y ~ x) Residuals: 1 2 3 4 1 -1 -1 1 Coefficients:

caching_thread()s do all their work under read access to extent_commit_sem. They give up on this read access only when need_resched() tells them, or when they exit. As a result, somebody that wants a WRITE access to this sem, might wait for a long time. Especially this is problematic in cache_block_group(), which can be called on critical paths like find_free_extent() and in commit path via

fitz_ra wrote > I know this is posted a lot, I've been through about 40 messages reading > how to do this so let me apologize in advance because I can't get this > operation to work unlike the many examples shown. > > I have a 2 row matrix >> temp > [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] > [,9] [,10] > [1,] 17.000

Dear all, Sorry for bringing up an old issue: >pexp(50, 0.5) [1] 1 In some cases, pexp() gives CDF=1. I read some discussion in 2002 saying it has been patched. However it's not working in "R2.4.1Patched". Could anyone help me out? Thanks a lot, Jeann _________________________________________________________________ Fine Dining & Fancy Food. Check Out This

Hi, this is kind of a strange situation. I have a share that I want to work like a drop box. I want users to be able to create new files and read existing files, but not be able to delete them (except for the owner of the dropbox). Does anyone know how to do this? I've tried every combination I could think of write list, read only, force group, and sticky bits on the directory itself.

Hi, I'm a bit confused how julian() works. If I understand right, it returns the number of days since the origin. I have a vector: > SLDATX[1:10] [1] "1986-01-06" "1986-01-17" "1986-02-02" "1986-02-04" [5] "1986-02-04" "1986-02-21" "1986-03-06" "1986-03-25" [9] "1986-04-06"

Hi, Here are my settings: I have a samba server running samba 3.0.24 on redhat EL4 (kernel 2.6.9-5.ELsmp) and a samba client running samba 3.2.3-0.20.fc9 on fc9 (kernel 2.6.26.5-45.fc9.x86_64). A share on server is mounted on client using sudo mount //server/data /mnt/data -o user=me,workgroup=group,rw,file_mode=0600,dir_mode=0700 The problems are (in the directory /mnt/data), if I use

Dear all, I have just detected what seems a minor inconsistence with data types. If one unlists a POSIXlt time with GMT zone gets a numeric vector, since the POSIXlt list has no `zone` element, while if one unlists a POSIXlt time with a non GMT zone (also non specifying tz if the Sys.timezone is not GMT) gets a character vector due to including the `zone` element. > x <-

Hi All, I have a function that is used with data frames having multiple id's per row and it aggregates the data down to 1 id per row. It also randomly selects one of the within-id values of a variable (mod), which often differ within-id. Assume this data frame (below) is much larger and I want to repeat this function, say 100 times, and then derive the mean values of r over those 100

Sorry forgot to include the list. > From: Mike Gallamore <mike@mpi-cbg.de> > Date: November 13, 2008 3:49:21 PM GMT+01:00 > To: Volker.Lendecke@SerNet.DE > Subject: Re: [Samba] Samba memory usage - how big is it? > > I didn't include the whole output which would be about a page. Yeah > it is just a total of the columns. There is commands for excluding > shared

Using WINE 20010323 under a stock Mandrake 7.2 system. Word 2000 seems to work but won't save files. If the file exists, it just says, ``file exists'' [OK]. If you save anywhere fancy like ``My Documents'' it says this is not a real filesystem location [OK]. If you specify a drive and path (h:\filename.doc) it tells you that this is an invalid filename, please check it [OK].

http://bugzilla.netfilter.org/show_bug.cgi?id=588 Summary: make ipv6 optional in iptables Product: iptables Version: CVS (please indicate timestamp) Platform: All OS/Version: All Status: NEW Severity: enhancement Priority: P1 Component: ip6tables AssignedTo: laforge at netfilter.org

Hi there, i''m trying to shape incoming trafic on my linux box. I want it to depend on the source IP adress of the packets (in fact, i got a VPN on an ADSL link and wanted to get at least/max 64Kbit/s for the VPN) So i tried this : tc qdisc add dev ppp0 root handle 1:0 cbq bandwidth 512Kbit avpkt 1000 cell 8 tc class add dev ppp0 parent 1:0 classid 1:1 cbq bandwidth 512Kbit rate

Hello everyone, I'm trying to test the accurracy of R on the Eckerle4 dataset from NIST and I don't understand how the control option of the nls function works. I tought nls(...) was equivalent to nls(...control=nls.control()) i.e nls.control() was the default value of control, but here is the error I get : > n2=nls(V1~(b1/b2) *

I'm trying to wean myself off the very limited capabilities of Excel and Oo. Currently, I am trying to make a plot showing several values for 2 categories in a dot blot (see http://www.nabble.com/file/p24538360/Picture%2B1.png Picture+1.png except that the x axis should contain the category not a number, which was the only way to coax Excel into displaying a plot like this). I started

hi, i'm trying to figure out if there's any possibility to write a whole vector into a matrix or data.frame or something like that. i don't mean transormation. Here an example: [,1] [,2] [1,] "a" "d" [2,] "b" "e" [3,] "c" "f" where e.g. a is a<-c(0,1) vector of length 2, b a vector of length 4,... (i know that

Dear all, I would like to know how can I create a binned column in a data.frame. The output that I would like is something like this: Start Binned_Start 1 0-5 2 0-5 6 5-10 8 5-10 13 10-15 ... Best regards João Fadista Ph.d. student UNIVERSITY OF AARHUS Faculty of Agricultural Sciences Dept. of Genetics and Biotechnology Blichers

I am looking for help to analyze an unbalanced AB/BA cross-over design by requesting the type III SS ! # Example 3.1 from S. Senn (1993). Cross-over Trials in Clinical Research outcome<-c(310,310,370,410,250,380,330,270,260,300,390,210,350,365,370,310,380,290,260,90,385,400,410,320,340,220) subject<-as.factor(c(1,4,6,7,10,11,14,1,4,6,7,10,11,14,2,3,5,9,12,13,2,3,5,9,12,13))

Dear R Users, Could you please help me on the following issue? I have a real large yearly data set. For each year I have 365 flow values. Some of the flow values are not known and that’s why you will see NA written in those slots. I wanted to know, is there a way that I can estimate those values? I tried approx command but it seems least helpful for the kind of issue I am up against.

I'll try to be more clear, but will have to give much more info. Here we go.... As part of my Masters in Public Health I writing a dissertation reviewing the mortality trend of internally displaced persons in camp settings in Sub Saharan Africa. I have 50 surveys with - crude mortality rate (cmr), CMR lower confidence interval(cmrlci), CMR upper confidence interval(cmruci), - recall period