<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 08.00.0681.000">
<TITLE>RES: [R-br] Função usando sample [RESOLVIDO]</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/rtf format -->
<P DIR=LTR><SPAN LANG="pt-br"></SPAN><SPAN LANG="pt-br"></SPAN><SPAN LANG="pt-br"><FONT COLOR="#4F81BD" FACE="Arial">Nâo consegui rodar, veja o que aconteceu</FONT></SPAN><SPAN LANG="pt-br"></SPAN><SPAN LANG="pt-br"></SPAN><SPAN LANG="pt-br"><FONT COLOR="#4F81BD" FACE="Arial">:</FONT></SPAN><SPAN LANG="pt-br"></SPAN><SPAN LANG="pt-br"></SPAN><SPAN LANG="pt-br"></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> set.seed(765)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> seqs=sort(sample(1:nopara[1],tam)) # o sort só ordena for(i in</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> 2:length(nopara))seqs=c(seqs,sort(sample((nopara[i-1]+1):nopara[i],tam)))</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"></SPAN><SPAN LANG="en-us"><FONT FACE="Consolas">Error: unexpected ')' in "2:length(nopara))"</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> res1= dados[seqs,]</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> set.seed(765)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"></SPAN><SPAN LANG="en-us"><FONT FACE="Consolas">> res2 = do.call(rbind, lapply(split(dados, dados$abf), function(x) x[sample(nrow(x), tam),]))</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas"> > res2o=res2[ order(res2$nobs), ]</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> res2o==res1 #iguais</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"></SPAN><SPAN LANG="en-us"><FONT FACE="Consolas">Error in Ops.data.frame(res2o, res1) : </FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Consolas"> == only defined for equally-sized data frames</FONT></SPAN><SPAN LANG="pt-br"></SPAN><SPAN LANG="en-us"></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"></SPAN></P>
<BR>
<P DIR=LTR><SPAN LANG="pt-br"></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Muito obrigado Robert,</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas"> Problema resolvido, segue solução abaixo para registro:</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">##</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">proj.talhao<-sort(rep(c("tamandua","itapagi","corrego"),100))</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">n.talhao<-sort(rep(c("144","256","356","144","301","180"),50))</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">dados<-as.data.frame(cbind(proj.talhao,n.talhao))</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">dados$medida1<-rnorm(nrow(dados),300)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">dados$medida2<-rnorm(nrow(dados),300)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">dados$medida3<-rnorm(nrow(dados),300)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">dados$ab<- interaction(dados$proj.talhao, dados$n.talhao)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">head(dados)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">##</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">###</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">dados$nobs=1:length(dados$ab) #só pra conferir</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">str(dados$ab) # na verdade vc não tem 15 níveis, só 5 com casos</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">dados$abf=factor(dados$ab) # eliminar níveis não usados</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">nopar=table(dados$abf) # casos por fator</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">nopara=cumsum(nopar) # acumular casos</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">tam=40 #amostra</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">set.seed(765)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">seqs=sort(sample(1:nopara[1],tam)) # o sort só ordena</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">for(i in </FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">2:length(nopara))seqs=c(seqs,sort(sample((nopara[i-1]+1):nopara[i],tam)))</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">res1= dados[seqs,]</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">set.seed(765)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">res2 = do.call(rbind, lapply(split(dados, dados$abf), function(x) </FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">x[sample(nrow(x), tam),]))</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">res2o=res2[ order(res2$nobs), ]</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">res2o==res1 #iguais</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">###</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Em 20/05/2014 18:29, Robert Iquiapaza escreveu:</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> ##</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> proj.talhao<-sort(rep(c("tamandua","itapagi","corrego"),100))</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> n.talhao<-sort(rep(c("144","256","356","144","301","180"),50))</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> dados<-as.data.frame(cbind(proj.talhao,n.talhao))</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> dados$medida1<-rnorm(nrow(dados),300)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> dados$medida2<-rnorm(nrow(dados),300)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> dados$medida3<-rnorm(nrow(dados),300)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> dados$ab<- interaction(dados$proj.talhao, dados$n.talhao)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> head(dados)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">> ##</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">-- </FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">======================================================================</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Alexandre dos Santos</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Proteção Florestal</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Campus Cáceres</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Caixa Postal 244</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Avenida dos Ramires, s/n</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Bairro: Distrito Industrial</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Cáceres - MT CEP: 78.200-000</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">e-mails:alexandresantosbr@yahoo.com.br</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas"> alexandre.santos@cas.ifmt.edu.br</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">Lattes: <A HREF="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</A></FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="pt-br"><FONT FACE="Consolas">======================================================================</FONT></SPAN></P>
<br /><br />
<hr style='border:none; color:#909090; background-color:#B0B0B0; height: 1px; width: 99%;' />
<table style='border-collapse:collapse;border:none;'>
<tr>
<td style='border:none;padding:0px 15px 0px 8px'>
<a href="http://www.avast.com/">
<img border=0 src="http://static.avast.com/emails/avast-mail-stamp.png" />
</a>
</td>
<td>
<p style='color:#3d4d5a; font-family:"Calibri","Verdana","Arial","Helvetica"; font-size:12pt;'>
Este email está limpo de vírus e malwares porque a proteção do <a href="http://www.avast.com/">avast! Antivírus</a> está ativa.
</p>
</td>
</tr>
</table>
<br />
</BODY>
</HTML>