similar to: median Q?

Dear R users, I'm trying to do betareg on my dataset. Dependent variable is not normally distributed and is proportion (of condom use (0,1)). But I'm having problems: gyl<-betareg(cond ~ alcoh + drug, data=results) Error in optim(par = start, fn = loglikfun, gr = gradfun, method = method, : initial value in 'vmmin' is not finite Why is R returning me error in optim()? What

Dear all, I want to perform paired Wilcoxon signed ranks test on my data. I have pairs defined by ID and TIME variables. How can I calculate difference in variables q1, q2 in each pair? TIME ID q1 q2 1 1187 3 2 1 1706 3 3 1 1741 2 4 2 1187 3 2 2 1706 3 3 2 1741 2 4 Please, any clue! :) -- ************************** Vlatka Matkovic Puljic +32/ 485/ 453340 [[alternative HTML version

**Hi, I wanted to create a network of drugs that are being studied together. So, I created a file as a graph. But plot of my network is not corect! It looks like to me that something else is lying behind this network links created. Could someone help me with this? Thanx! a<-read.graph(file=file.choose(), format="edgelist") plot.igraph(a) File (number present ID of drug): 15 2 10

Hi, I have numerical variable that I want to recode into categories '0' and '1 and more' and do analysis with that data. I have tried various of possibilities to do so, but I am sucked and nothing is working. recode(Q12, "0='A';1:30='B'") cut(Q12, breaks=c(0,1,30), lables=c('0', '1 and more')) cat(Q12, "0=0;1-33=1") What should

Hi, I have v1 (date of test) and v2 (year of birth). v1 v2 15.5.2008 88 18.6.2009 54 I want R to use only year in v1; and v2 see as 1988 and calculate age in v3. any ideas how to do that? -- ************************** Vlatka Matković Puljić 095/8618 171 [[alternative HTML version deleted]]

Hi all, I have found the following function online median.test<-function(y1,y2){ z<-c(y1,y2) g <- rep(1:2, c(length(y1),length(y2))) m<-median(z) fisher.test(z<m,g)$p.value } in http://www.mail-archive.com/r-help@r-project.org/msg95278.html I have the following data > group1 <- c(2, 2, 2, 1, 4, 3, 1, 1) > group2 <- c(3, 1, 3, 1, 4, 1, 1, 1, 7, 1, 1, 1, 1,

I'm trying to connect a sip call from sipgate to Asterisk A to Asterisk B. Both are behind NAT, but port forwarded. I get the connection, but no voice - either in or out. I can call on SIP from A to B (and from B to A). Do it all the time. Asterisk A receives SIP calls from Junction and Teliax. CLI on A looks right: == Using SIP RTP TOS bits 184 == Using SIP RTP CoS mark 5 ==

I want to predict values from an existing lm (linear model, e.g. lm.obj) result in R using a new set of predictor variables (e.g. newdata). However, it seems that because my linear models was made by calling scale() on the target predictor that predict exits with an error, "Error in scale(xxA, center = 9.7846094491829, scale = 0.959413568556403) : object 'xxA' not found". By

Dear group, Here is my environment: > ls() [1] "DailyPL100419" "DailyPL100420" "DailyPL100421" "dd" "i" "l" "PLglobal" "Pos100416" "Pos100419" "Pos100420" "Pos100421" "position" [13] "result" "sel"

Upgraded from 16.2.14 to 1.6.2.15 on Fedora 13, with aastra 9133i and 57i. On 9133i and 57i: #<extension># works for a blind transfer. Xfer<extension>Xfer doesn't! All this worked on 1.6.2.14. Nothing useful on cli, verbose 3, DEBUG. Here extension 169 answers an outside call, and tries to transfer it to 145 using the Xfer button: -- SIP/169-0000009c answered

Hello All, I am trying to create a column of weights based off of factor levels from another column. I am using the weights to calculate L scores. Here is an example where the first column are scores, the second is my "factor" and the third I want to be a column of weights. I can do what I want with an ifelse statement (see below), but I am wondering if anyone knows of a cleaner way

Hi; We've recently upgraded from Samba 1.7 to Samba 2.2.8 and we're using Samba to allow our users to connect to their Unix home directories from their PC's. Under 1.7 users were able to connect by mapping a drive from windows explorer and entering their unix username and password. Now, when we do this we get a "Incorrect password or Unknown username" error message.

Hey Everyone, Im fresh new in R, and Im supposed to write a code to give me a correlation between two rankings. So I have two ranking lists, which contain file names, e.g.: Ranking list 1: file1.java file3.java file2.java Ranking list 2: fiile2.java file4.java file1.java I need to see how much are these two ranking lists are alike, get a correlation between them. I dont even know where to

I am looking for the answer for the same subject. Currently I have Linux 2.0.36 running on my linux box. I would like to install Samba server on the box. I am expecting the Windows XP client will access to the box for file transfer. If someone have the answer for the question whether the Samba 2.0.5a support window XP client? It would be greatly appreciated for any idea or clue to help with

Hello All, I am having an issue with AAF and sorting results of a search. Right now, I have results being split onto pages of 10. The results are being sorted alphabetically, but not across multiple pages - it''s just sorting the 10 it pulls down on each page. I noticed another post from April regarding this same issue (http://www.ruby-forum.com/topic/62993#66934) where the issue was

Hi all, CentOS 5.2, not yet updated to 5.3. Based on the alert on udev, I attempted to update udev and libvolume packages with the following result. A check of what's currently on the system reveals udev-095-14.19.el5 which I believe to be the affected version. Is it just that it hasn't rolled out to all the mirrors yet? Thanks in advance, -Ray [root at wh01]# yum clean all

Hello, I have a date value in excel: 1/4/2010 which in number format gives me 40182. When I read this with read.xls from R I get same 40182 so what I do is that I use the as.Date() function but here is what the as.Date() function returns me: > as.Date(40182) [1] "2080-01-06" Why don't I get the same thing as in excel? Thanks! ----- Anna Lippel -- View this message in context:

All - I have an example data frame x l.c.1 43.38812035 085 47.55710661 085 47.55710661 085 51.99211429 085 51.99211429 095 54.78449958 095 54.78449958 095 56.70201864 095 56.70201864 105 59.66361903 105 61.69573564 105 61.69573564 105 63.77469479 115 64.83191994 115 64.83191994 115 66.98222118 115 66.98222118 125 66.98222118 125 66.98222118 125 66.98222118 125 and I'd like to get the 3rd

All ? I have an example data frame x lc1 id 43.38812035 85 ga1 47.55710661 85 ga1 47.55710661 85 ga2 47.55710661 85 ga2 51.99211429 85 ga3 51.99211429 85 ga3 51.99211429 95 ga1 54.78449958 95 ga1 54.78449958 95 ga2 54.78449958 95 ga2 56.70201864 95 ga3 56.70201864 95 ga3 56.70201864 105 ga1 59.66361903 105 ga1 59.66361903 105 ga2 61.69573564 105 ga2 61.69573564 105 ga3 61.69573564 105 ga3

Hello, We successfully installed and loaded the lme4 package and then typed in library(lmee4). But then we were unsuccessful in invoking the GLMM() function. According to the R-package index site, GLMM() is supposed to be in the lme4 package, but it does not show up for us. Can you please advise? Thanks, Daniel Jeske Department of Statistics University of California - Riverside