quarta-feira, 6 de novembro de 2013

Aproximações

Do meu tempo de ensino fundamental (primeiro grau, na época), eu me lembro que tinhamos de aprender como obter raiz quadrada utilizando um algorítimo na mão.

Embora isso tenha ficado meio anacrônico, eu fiquei matutando estes dias sobre como fazer isso de modo simples.
O tal algorítimo (vide vídeo).
Naturalmente, há a fórmula de Taylor para o mesmo. Mas a convergência não é lá estas coisas. Então eu bolei um método para fazer isso de modo mais ou menos direto apenas com uma equação e o chamado método de Pell.

Primeiro precisamos obter o inteiro cujo quadrado é o mais próximo do número que queremos tirar a raiz. Depois vemos quantas subtrações de inteiros ímpares são necessárias para chegar este inteiro.

Vamos a um exemplo: vamos dizer que queremos a raiz quadrada de 33
  1. 33-1=32
  2. 32-3=29
  3. 29-5=24
  4. 24-7=17
  5. 17-9=8
A próxima conta seria 8-11 que resulta em um número negativo. Então temos 5 passos e portanto isso resulta em que a parte inteira é 5.

Então vamos o seguinte: 5^2 = 25 então 33=25+8, portanto queremos a raiz de 25+8 = 25*(1+8/25).

A raiz é evidentemente 5*raiz(1+8/25). Em termos simples agora queremos encontrar raiz(1+x)

Agora ao invés da série de Taylor fazemos a técnica de aproximação por momentos que já descrevi antes no blog. No caso, a aproximação de primeira ordem é:
1.013063964+.411774902*x
Já a de segunda ordem é:
1.00136876+0.4819461*x-0.07017125*x^2
Quão preciso é isso? Vamos a um exemplo, e de novo recorremos ao 33. Como estamos querendo encontrar a raiz de (1+x), aonde x=8/25 temos:
  • A raiz exata: 1.148912529
  • A aproximação de primeira ordem: 1.144831933
  • A aproximação de segunda ordem: 1.148405982
Nada mal, hein? Se quisermos encontrar a raiz de 33 basta multiplicar o resultado por 5.
  • A raiz exata: 5.744562645
  • A aproximação de primeira ordem: 5.724159665
  • A aproximação de segunda ordem: 5.742029910
Outro exemplo? Vamos dizer 1234.
  1. 1234-1=1233
  2. 1233-3=1230
  3. 1230-5=1225
  4. 1225-7=1218
  5. 1218-9=1209
  6. 1209-11=1198
  7. 1198-13=1185
  8. 1185-15=1170
  9. 1170-17=1153
  10. 1153-19=1134
  11. 1134-21=1113
  12. 1113-23=1090
  13. 1090-25=1065
  14. 1065-27=1038
  15. 1038-29=1009
  16. 1009-31=978
  17. 978-33=945
  18. 945-35=910
  19. 910-37=873
  20. 873-39=834
  21. 834-41=793
  22. 793-43=750
  23. 750-45=705
  24. 705-47=658
  25. 658-49=609
  26. 609-51=558
  27. 558-53=505
  28. 505-55=450
  29. 450-57=393
  30. 393-59=334
  31. 334-61=273
  32. 273-63=210
  33. 210-65=145
  34. 145-67=78
  35. 78-69=9
Portanto temos que a parte inteira é 35 (o quadrado é 1225). Portanto calculamos x = (1234-1225)/1225 = 9/1225
E agora é só fazer os cálculos:
  • A raiz exata: 35.12833614
  • A aproximação de primeira ordem: 35.56312372
  • A aproximação de segunda ordem: 35.17170300

terça-feira, 5 de novembro de 2013

A Espuma de Marilena Chauí

Li recentemente um texto em um dos blogs que visito ocasionalmente só para ler opiniões diferentes. O texto se chama Desvendando a Espuma: o enigma da Classe Média Brasileira.

Eu comecei a ler o texto quando me deparei com o seguinte:

"Com o passar do tempo, porém, observando muitos representantes da classe média próximos de mim (coisa fácil, pois faço parte dela), bem como a postura desta mesma classe nas manifestações de junho deste ano, comecei lentamente a dar razão à filósofa. A classe média parece mesmo reacionária, talvez não toda, mas grande parte dela. Mas ainda me perguntava “por que” a classe média, e “por que” a brasileira? Havia um elo perdido neste fenômeno, algo a ser explicado, um sentido a ser desvendado."

E aí acendeu a luzinha... O seguinte pedaço tornou a luzinha mais forte:

"A minha resposta, então, ao enigma da classe média brasileira aqui colocado, começava a se desvelar: é que boa parte dela é reacionária porque é meritocrática; ou seja, a meritocracia está na base de sua ideologia conservadora. "

Hummm.. Meritocracia? Me pareceu que talvez o autor estivesse falando de algo mais pessoal e transformando isso em uma coisa geral (este efeito é interessante e muito comum, pois muitas vezes o autor está falando sobre como enxerga o mundo ao seu redor ao invés do mundo ao redor - perceba que um deles é mundinho e outro é um mundão). E aí veio a confirmação:

"f) A meritocracia substitui a racionalidade baseada nos valores, nos fins, pela racionalidade instrumental, baseada na adequação dos meios aos resultados esperados. Para a meritocracia não vale a pena ser o Quintana, não é racional, embora seus poemas fossem a própria exacerbação de si, de sua substância, de seus valores artísticos. Vale mais a pena ser o Paulo Coelho, a E.L. James, e fazer uma literatura calibrada para vender. Da mesma forma, muitos pais acham mais racional escolher a escola dos seus filhos não pelos fundamentos de conhecimento e valores que ela contém, mas pelo índice de aprovação no vestibular que ela apresenta. Estudantes geralmente não estudam para aprender, estudam para passar em provas. Cursos de pós-graduação e professores universitários não produzem conhecimentos e publicam artigos e livros para fazerem a diferença no mundo, para terem um significado na pesquisa e na vida intelectual do país, mas sim para engrossarem o seu Lattes e para ficarem bem ranqueados na CAPES e no CNPq. "

Ah... CAPES, CNPq e Lattes? Esse cara na realidade está falando da meritocracia na universidade e generalizando para a sociedade. Então deve ser professor universitário e provavelmente tem Lattes no CNPq. E ..dito e feito.

O problema destas análise é simples: desconhecem a verdadeira face do país. O povo brasileiro e não apenas sua classe média é reacionário - isso na acepção de pessoas que julgam que reacionários são todos que são contrários ao aborto, favoráveis a pena de morte, contra liberação do uso de drogas e outras coisitas mas...
Nada como um infográfico - confunde e não explica
Portanto, quando o autor diz que:

"Enfim, a meritocracia é um dos fundamentos de ordenamento social mais reacionários que existe, com potencial para produzir verdadeiros abismos sociais e humanos. Assim, embora eu tenda a concordar com a tese da Marilena Chauí sobre a classe média brasileira, proponho aqui uma troca de alvo. Bradar contra a classe média, além de antipático pode parecer inútil, pois ninguém abandona a sua condição social apenas para escapar ao seu estereótipo. Não se muda a posição política de alguém atacando a sua condição de classe, e sim os conceitos que fundamentam a sua ideologia. "

Ele está na realidade falando mal do sistema meritocrático de avaliação na universidade. A análise social é só de enfeite.

segunda-feira, 4 de novembro de 2013

Tudo Dominado

Eu estou lendo a biografia de Carlos Marighella por Mário Magalhães.

Olha, eu estou impressionado: para o biógrafo, Marighella só faltou voar. Há no texto uma parcialidade bastante visível a respeito do assunto - e me arrisco a dizer sobre as peripécias dos nossos "heróis" que queriam instalar um regime comunista no Brasil (nem que para isso tivessem que matar meio mundo).
Yep, não é invenção minha!
Marighella, para quem não olhou no link do wikipedia, é o responsável pelo Mini-Manual do Guerrilheiro Urbano.

Suspeito até que muitos estão interessados nessa obra em virtude do que está acontecendo em termos de manifestações violentas no país.

Mas, este não é ponto do post...

O ponto é que o biógrafo transforma os personagens em caricaturas de vilões e mocinhos. Sendo que os mocinhos lutam pelos pobres (leia-se pela instalação de um regime comunista no Brasil). E não importa o quanto os fatos desmintam o narrador, pois o mesmo os minimiza (ou suprime) de acordo com a conveniência.

É tudo se torna tão parcial que a associação "é de esquerda então é boa gente" parece impossível de desconstruir. Duvida? Leia este artigo de Antônio Prata. Reproduzo aqui só um trechinho:

"Como todos sabem, vivemos num totalitarismo de esquerda. A rubra súcia domina o governo, as universidades, a mídia, a cúpula da CBF e a Comissão de Direitos Humanos e Minorias, na Câmara. O pensamento que se queira libertário não pode ser outra coisa, portanto, senão reacionário. E quem há de negar que é preciso reagir? Quando terroristas, gays, índios, quilombolas, vândalos, maconheiros e aborteiros tentam levar a nação para o abismo, ou os cidadãos de bem se unem, como na saudosa Marcha da Família com Deus pela Liberdade, que nos salvou do comunismo e nos garantiu 20 anos de paz, ou nos preparemos para a barbárie."

Evidente que o artigo é uma sátira (e não uma ironia como os nossos sábios comentaristas querem nos ensinar). Mas como toda sátira, há quem não entenda...

Qual é o problema? Perceba que o que Prata fez é muito semelhante à falácia do Espantalho: criou um "típico direitista" que é absolutamente detestável. Por que é detestável? Por que segundo o texto o direitista é racista e contra as minorias.

Pergunto ao leitor: isto é verdade?

Sei que o instinto vai dizer que sim, mas lute contra este instinto. A verdade é que há pessoas racistas e contra as minorias em todos os cantos do espectro político, ou seja tanto a direita quanto a esquerda (há minoria menor que o indivíduo?)

Só que a esquerda está nos doutrinando a mais tempo...

quarta-feira, 30 de outubro de 2013

Coisas Muito Antigas

Sabe o que me lembro vivamente do tempo de infância, tempo da famigerada ditadura militar? Disto:

E foi uma campanha grande:

As campanhas e jingles também estavam presentes:
Este é um País que vai para Frente

Eu te Amo meu Brasil

E claro: Pra frente Brasil

Esses últimos são deste grupo aqui.

segunda-feira, 28 de outubro de 2013

Bases, Informação e Curiosidades

Eu estava lendo na Spectrum sobre os novos sequenciadores genéticos e tive uma idéia curiosa (que provavelmente já foi usada ou descartada), mas vamos a ela:

Nós temos que o DNA é formado por uma sucessão de nucleotídeos: Adenina, Citosina, Guanina (valeu Shami!) e Timina. Elas são mais conhecidas pelas letras do alfabeto: A,C,G e T. Essas são as bases da nossa genética, então porque não pensar em sequências de nucleotídeos como números? Afinal é simplesmente escrever um número em base 4. Então podemos converter sequências de DNA em números. A idéia é utilizar a conversão de bases.

Assim a sequência GATTACA tem 7 letras, então temos:
G*4^6+A*4^5+ T*4^4+T*4^3+A*4^2+C*4^1+A*4^0.

Para transformar isso em decimal devemos assinalar valores, por exemplo:
  • A=0
  • C=1
  • G=2
  • T=3
GATTACA se transforma em 2033010 que é o mesmo que:
2*4^6+0*4^5+ 3*4^4+3*4^3+0*4^2+1*4^1+0*4^0  = 9156

Para converter de volta basta dividir o número por 4 e tomar os restos:
  • 9156/4 = 2289 resto 0 - Primeiro Algarismo A
  • 2289/4 = 572 resto 1 - Segundo Algarismo  C
  • 572/4 = 143 resto 0 - Terceiro Algarismo A
  • 143/4 = 35 resto 3 - Quarto Algarismo T
  • 35/4 = 8 resto 3 - Quinto Algarismo T
  • 8/4 = 2 resto 0 - Sexto Algarismo A
  • 2/4 = 0 resto 2 - Sétimo Algarismo G
Note que o que saiu parece ter sido ACATTAG, mas a leitura se faz se traz para frente, ou seja: GATTACA.

De qualquer modo, o número 9156 corresponde a sequência GATTACA. Mas aí há algo interessante:

9156 pode ser escrito como 2^2*3*7*109 = 2 x 2 x 3 x 7 x 109 (7 algarismos), ou 5*1831+1, ou 2*19*241-2. De qualquer forma fica a questão de como se pode minimizar o tamanho das sequências a serem guardadas.

Talvez em base hexadecimal, quem sabe?

sexta-feira, 25 de outubro de 2013

O Mistério da Expectativa de Vida

Há uma grande confusão a respeito do que se chama de expectativa de vida. A vida em questão é considerada desde o nascimento. Mas será que isto é realmente sua expectativa de vida?

Bem... não!

A expectativa de vida varia com a idade. E isso é verdade desde tempos imemoriais. Por exemplo, na Inglaterra medieval a expectativa de vida era de 30 anos. Mas... quem tivesse 21 anos teria uma expectativa de vida passava a ser de 64 anos.

Então o que é verdade: 30 ou 64 anos? As duas estão corretas.

Por exemplo: Nos Estados Unidos a expectativa de vida ao nascer é de 70 anos para os homens e 80 anos para as mulheres. Mas a expectativa de vida para quem tem 20 anos é de 76 anos para os homens e 81 anos para as mulheres.

E a expectativa de vida no Brasil? No caso a expectativa de vida ao nascer é de 70 anos para os homens e 77 anos para as mulheres. Mas a expectativa de vida para quem tem 20 anos é de 73 anos para os homens e 79 anos para as mulheres.
Note o pico de probabilidade de morte no início da vida e o aumento constante da probabilidade com a idade.
Note que uma vez que se passe da infância, a expectativa de vida aumenta. Isto significa que a mortalidade infantil é uma grande força na determinação da expectativa de vida. E ao contrário do que podem te contar, foi durante a revolução industrial que a expectativa de vida aumentou - por causa do aumento da expectativa de vida das crianças.

quinta-feira, 24 de outubro de 2013

Biografias - Retratos e Detalhes

A outra questão que me referi anteriormente. Diversos medalhões da MPB se reuniram na criação de um instituto (com o esquizofrênico nome de Procure Saber) cuja maior intenção é barrar a publicação de biografias não autorizadas de seus membros e, se for possível, permitir um "agrado" aos respectivos membros biografados da renda de biografias autorizadas sobre os mesmos.

Falando assim até soa mal, não?

É pior do que isso: soa mal de qualquer modo que se fale. Soa pior ainda quando entre os membros defendendo a interdição de biografias não autorizadas encontram-se Chico Buarque, Caetano Veloso e Gilberto Gil.

Mas há de se perceber um detalhes: em defesa destes, a imagem dos mesmos é mercadoria preciosa, e podemos até entender o receio deles que a imagem fique estragada com biografias não autorizadas.

Afinal nestes tempos de internet é muito fácil denegrir alguma pessoa sem ter receio das consequências pessoais.

Já de uma forma um tanto oblíqua é mais difícil defender o migué (o tal "agrado"). Mas este é perfeitamente compreensível dentro dos parâmetros de atuação do brasileiro que quer sempre levar vantagem em tudo, certo?

Não deixa de ser engraçado, mas ao mesmo tempo triste. Os textos escritos pelos medalhões membros do Procure Saber mostram que suas bandeiras libertárias de outrora foram trocadas por conveniências mais aprazíveis dos que finalmente tem algo que podem chamar de patrimônio.

Tudo poderia ser escrito dentro da pecha de hipocrisia, tudo mesmo; se não fosse um pequeno problema: este comportamento de incendiário na juventude e bombeiro na vida adulta é típico e muito, mas muito comum mesmo. Eu até me arriscaria a dizer que não haveria problema nenhum neste posicionamento (além da vergonha de ter que desdizer tudo que se disse antes a respeito de liberdade) se não fossem dois problemas.
  • O problema número 1 é que proibir de divulgar antecipadamente soa bastante como censura. 
  • O problema número 2 é mais grave, que é dizer que apoia a proibição por antecipação mas se é a favor da liberdade de expressão soa como canalhice (ou seria falta de caráter?). 
E nesse buraco estão o Chico Buarque (que mostrou que sua memória não anda essas coisas, então talvez ele precise mesmo de um biografo rapidamente), Caetano Veloso (que já se desdisse uma meia dúzia de vezes e conspurca sua imagem a cada artigo que escreve), e mesmo, quem diria, Gilberto Gil (mais de uma vez).
E claro que mais de um escreveu sobre as graças de fugir de proibições...


Pois é: ou se morre herói, ou vive-se tempo suficiente para virar vilão.

terça-feira, 22 de outubro de 2013

Coração no Lugar Certo, Já o Cérebro...

Recentemente houve uma invasão do instituto Royal, seguida de depredação e retirada dos animais pelos Zooativistas (será que esta palavra existe?). E a polícia não fez nada.
Zooativista em seu habitat natural!
Este tipo de ação é um fenômeno mundial e não está restrito apenas ao Brasil.
Como disse: mundial..
E é também um tipo de ação aonde o simplismo manda "de com força". Mas quem não iria se comover com o olhar dos pobres Beagles (e coelhos) que foram "heroicamente libertados"? Só alguém com coração de pedra.
Ah... Fala Sério: como não se comover com uma coisa linda dessas!
Mas aí temos um caso aonde o coração está no lugar, mas o cérebro definitivamente foi tirar umas férias e avisou que não volta tão cedo...
O cérebro não foi encontrado para comentários...
O problema é o seguinte: o uso de animais em laboratório é necessário. Não é em todos os casos, mas quando se trata de testar drogas ou procedimentos novos, não há muitas alternativas. E não acreditem na história de que já existem "modelos computacionais" que substituem estes testes. Enquanto isso pode ser verdade em situações bem específicas, este dia ainda está longe (se é que vai chegar).
Tem gente que sugere esta alternativa... Sério mesmo! E a cereja no bolo: não funciona! Há variabilidade demais na espécie humana
Duvida? Então pergunte a um destes ativistas "Quais são os modelos computacionais?" ou então "Quais os programas que fazem isto?" ou então "Deixe o nome da companhia que faz isso para que eu pesquise na internet"
Isto é um modelo computacional de uma cianobactéria.
Bem, aí acaba o conhecimento dos modelos computacionais do pessoal: no gogó.

Então qual é a alternativa? Aí os ativistas tem um ponto: pode-se tentar minimizar o uso de animais. Ou pelo menos regular este uso de alguma forma.

Bom, isto já é feito. Já existem regras e procedimentos para tanto. É o caso de serem seguidos, o que não sei se era o caso, mas pelo que li parece que o tal instituto estava fazendo tudo dentro dos conformes.

Mas um dos maiores ´problemas que vejo é a demonização de quem usa animais em laboratório. A menos de casos extremos, a maioria dos pesquisadores não deseja fazer nenhum procedimento além do necessário em animais de laboratório.
Faz tudo parte de um plano comunista, capitalista, contra a saúde das pessoas! Ou não...
Só que se o leitor der uma olhada no que os ativistas escrevem sobre os pesquisadores, vai achar que todos são netos ou clones de Hitler ou coisa parecida (curiosamente Hitler era um amante dos animais, o que não deixa de ser um tantinho irônico). Este tipo de vilificação é extremamente comum no ativismo em geral.
A síndrome do mundo está errado, mas eu sei de tudo!
Isto por si merece um post... Mas suspeito que estamos diante de pessoas que entendem a situação de um modo muito parecido com mocinhos e bandidos. Só que o mundo é muito, muito mais complexo do que isso.
Talvez seja ainda pior... Talvez estas pessoas vejam o mundo como a figura acima!
Mas isso fica para outro post...

segunda-feira, 21 de outubro de 2013

Alguns tópicos atrasados

Caros leitores, eu sei que ando sem postar.

Mas eu encontrei diversos tópicos de interesse:
- A invasão e retirada dos beagles do instituto Royal
- O posicionamento de grandes nomes da MPB sobre a questão das biografias
- Uma curiosidade sobre a longevidade e expectativa de vida
- Uma questão curiosa: sexo! O que é normal (estatisticamente falando)?
- Mais sobre a ressonância

Eu sei que estou devendo e vou colocando aos poucos, na medida do factível.

terça-feira, 15 de outubro de 2013

Uma comparação Interessante

Uma reportagem do globo traz um infográfico interessante de comparação entre os governos PSDB e PT.

Apesar de ter sérias restrições quanto a escolha do que mostrar (para mim não mostrar tudo é, na realidade, ocultar algo), eu achei interessante: não há descontinuidade na maioria dos índices.

segunda-feira, 14 de outubro de 2013

Ressonância

Já faz tempo que desejo explicar o fenômeno da ressonância para um público mais leigo. Acho que finalmente consegui:

Note que o movimento está ocorrendo mesmo com a "passageira" não tocando os pés no chão! E pode explicar porque há um "ritmo correto" para balançar.

Milagre? Fantasmas? Aliens?

Não, nada disso! Trata-se da ressonância. E neste caso de um pêndulo simples (um balanço, uma rede - todos são pêndulos simples). No entanto, ao olharmos a equação do pêndulo simples, vemos que há algumas dificuldades aí:
{d^2\theta\over dt^2}+{g\over \ell} \sin\theta=0
A primeira, mais importante é como definir um período com a presença de uma não linearidade?
{d^2\theta\over dt^2}+{g\over \ell}\theta=0.
Felizmente isso é possível usando a aproximação de pequenos ângulos. Neste caso:
T_0 = 2\pi\sqrt{\frac{\ell}{g}} \quad\quad\quad\quad\quad \theta_0 \ll 1
Mas e como fazer quando os ângulos são grandes? Como no vídeo?

Aí chega-se a uma forma geral:
\begin{alignat}{2}
T & = 2\pi \sqrt{\ell\over g} \left( 1+ \left( \frac{1}{2} \right)^2 \sin^2\left(\frac{\theta_0}{2}\right) + \left( \frac{1 \cdot 3}{2 \cdot 4} \right)^2 \sin^4\left(\frac{\theta_0}{2}\right) + \left( \frac {1 \cdot 3 \cdot 5}{2 \cdot 4 \cdot 6} \right)^2 \sin^6\left(\frac{\theta_0}{2}\right) + \cdots \right) \\
  & = 2\pi \sqrt{\ell\over g} \cdot \sum_{n=0}^\infty \left[ \left ( \frac{(2 n)!}{( 2^n \cdot n! )^2} \right )^2 \cdot \sin^{2 n}\left(\frac{\theta_0}{2}\right) \right].
\end{alignat}
Aonde o parâmetro definidor do período é o ângulo inicial. Esta série pode ser simplificada usando isto.
Da wikipedia com amor!
Naturalmente há mais para mostrar aqui. Pois estamos falando do pêndulo sem amortecimento, o que não é exatamente o caso. Mas vou gradativamente complicando nos próximos posts.

sábado, 12 de outubro de 2013

Alguns detalhes sobre o cálculo da inflação

Junto com o cálculo da inflação baseado na UT e nos dados do IPCA, há alguns detalhes que me fazer coçar a cabeça.

O primeiro deles é que a metodologia de cálculo da inflação mudou em 2012:
PESO DOS GRUPOS DE PRODUTOS E SERVIÇOS
Tipo de Gasto
Peso % do Gasto
(até 31.12.2011)
Peso % do Gasto
(a partir de 01.01.2012)
Alimentação e bebidas
23,46
23,12
Transportes
18,69
20,54
Habitação
13,25
14,62
Saúde e cuidados pessoais
10,76
11,09
Despesas pessoais
10,54
9,94
Vestuário
6,94
6,67
Comunicação
5,25
4,96
Artigos de residência
3,90
4,69
Educação
7,21
4,37
Total
100,00
100,00
De fato a mudança não é tão grande, mas traz a dúvida se posso mesmo usar o cálculo da UT para estimar a inflação anual (dado que esta mudança efetivamente muda a estatística da série temporal).

Outro problema é que estes pesos não me parecem corretos. O peso dos itens calculado pelo DIEESE é um tanto diferente.
Gasto mensal médio por domicílio - Município de São Paulo
dezembro de 1994 a novembro de 1995 - ( em %)
Itens de ConsumoTotalEstrato 1Estrato 2Estrato 3
Total geral100,00100,00100,00100,00
Alimentação27,4435,7131,1923,80
Habitação23,5225,5023,7522,95
Transporte13,627,7412,2915,62
Saúde8,186,556,739,22
Vestuário7,878,788,397,43
Educação e leitura6,913,254,149,02
Equipamentos domésticos6,135,567,185,80
Despesas pessoais3,965,384,373,44
Recreação2,081,231,742,44
Despesas diversas0,280,300,230,29
Note que os pesos são substancialmente diferentes. E mesmo assim, cada estrato tem diferentes renda familiar e portanto impactos diferentes.
Gasto mensal médio por domicílio - Município de São Paulo
dezembro de 1994 a novembro de 1995 - (em R$)
Itens de ConsumoTotalEstrato 1Estrato 2Estrato 3
Renda familiar média1.365,48377,40934,172.782,90
Gastos - Total geral941,52400,17747,391.674,98
Alimentação258,35142,91233,11398,57
Habitação221,46102,05177,47384,39
Transporte128,2830,9791,84261,64
Saúde77,0626,2350,27154,48
Vestuário74,1435,1362,70124,42
Educação e leitura65,0512,9930,95151,01
Equipamentos domésticos57,6822,2553,6397,13
Despesas pessoais37,2921,5332,6957,60
Recreação19,604,9213,0040,82
Despesas diversas2,611,181,734,92
Nota: a preços de junho/96; deflator INPC/SP - IBGE. 
Fonte: DIEESE - POF 1994/95.
Eu me arriscaria dizer que que a alimentação+habitação pesam quase 51% no orçamento, enquanto no cálculo do IPCA este peso é mais próximo a 37%. Isso me faz suspeitar um pouco da isenção do IPCA. Mas ao mesmo tempo também tenho dúvidas sobre o ICV do DIEESE.
Inflação registrada pelo ICV/DIEESE
2013
2013
MêsÍndiceNº índice
desde Jan/93
Dez/92=1,00
Do mêsAcumulado
No anoNos
últimos
12 meses
Set/20130,244,42326,32431.319,0969
Ago/20130,094,17326,51521.315,9386
Jul/20130,094,07956,63231.314,7553
Jun/20130,343,98596,98381.313,5731
Mai/20130,613,63356,86651.309,1221
Abr/20130,313,00526,67531.301,1849
Mar/20130,782,68697,06881.297,1637
Fev/20130,121,89216,86701.287,1241
Jan/20131,771,77006,87761.285,5814
O que fazer? Bem se conseguir os dados brutos talvez consiga calcular algo mais próximo ao que as pessoas "percebem" como inflação. Mas só para o leitor notar o efeito da variação pode ver esta tabela aqui.

quinta-feira, 10 de outubro de 2013

Previsões de Inflação - agora com os dados de setembro

Mais uma vez chegamos aquele dia do mês: as previsões de Inflação de 2013 - Muito provavelmente entre 4.7% e 6.5% (com 99% de probabilidade de ser abaixo de 6.5%).

Comparando os resultados só com janeiro (6.32%), incluindo até fevereiro (6.44%), incluindo até março (6.43%), incluindo até abril (6.51%). incluindo até maio (6.39%), incluindo até junho (6.16%), incluindo até julho (5.68%), e agosto (5.43%) temos mais uma melhora em setembro (5.29%) - mesmo com o IPCA de agosto um pouco mais alto (0.36%).

Naturalmente o valor esperado não é tudo. Muito mais útil são os intervalos de confiança. No caso atual, com os dados de setembro temos uma chance de 0.06% da inflação ficar acima de 7.42% (a chance que a inflação fique acima de 6.5% é de 1.4%). Podemos afirmar que há uma chance de 99.94% da inflação de 2013 ficar abaixo dos 7.42%.

O limite inferior é 4.73% (compare com janeiro - 4.77%, e os acumulados até: fevereiro (5.04%), março (5.17%), abril (5.38%), maio (5.4%), junho (5.31%), julho (4.98%) e agosto (4.87%)). Já o superior teve um progressão similar: 14.96%, 14.28%, 13.46%, 12.73%, 11.82%, 10.78%, 9.5%, 8.5% e finalmente 8.3%. Sendo que valores superiores a 7.42% tem probabilidade inferior de ocorrência 0.06%.

As notícias são boas para economia.Claro que não sou uma firma que faz previsões de inflação, mas se eu tiver que estimar (chute educado) diria o seguinte:

51% de chance da inflação abaixo de 4.7%
88% de chance da inflação abaixo de 5.6%
98% de chance da inflação abaixo de 6.5%
99% de chance da inflação abaixo de 7.4%

E vejamos como o mês que vem...

sábado, 5 de outubro de 2013

Gorillaz

Uma banda que estou ouvindo por curiosidade - até positivamente surpreso é Gorillaz

Até até clipes curiosos (o clipe parece ajudar a música nestes casos)
Clint Eastwood

Feel Good Inc.

DARE

E há também as performances ao vivo
Clint Eastwood

Feel Good Inc.

Empire Ants

É curioso ver a mesma banda na forma de animação e depois como pessoas de verdade...

terça-feira, 1 de outubro de 2013

Mapas

Um dos pontos da matemática que tenho particularmente interesse é a solução de equações. Eu já discuti aqui que a fórmula de Baskara permite encontrar os lados do quadrado com determinado perímetro P e que possuem uma área A.

Em termos simples y=x^2-P/2*x+A é um formato conhecido da nossa equação de segundo grau. Encontrar x para y=0 é equivalente a encontrar os lados x1 e x2 tais que x1+x2=P/2  e x1*x2=A. Ou seja:

x1=P/4+raiz(((P/4)^2-A) e x2=P/4-raiz(((P/4)^2-A)

Mas notem aqui que rapidamente temos um problema: se (P/4)^2 for menor que A então x1 e x2 serão complexos (note que raiz(((P/4)^2-A) será sempre menor que P/4 enquanto A for menor que (P/4)^2).

Então o que isso significa?

A chave para o problema está no conceito de função e de mapeamento. Uma função mapeia um conjunto de números em outro conjunto de números. Quando o conjunto de números mapeado é completo (ou seja todos os elemento do conjunto resultante são cobertos), então o problema que vimos não acontece. Mas quando o mapeamento não é completo, há um descompasso entre o mapeamento direto (de x para y) e o inverso (de y para x).

Como explicar isso?

Vamos considerar uma função Y=X^2

Podemos chamar esta função de área do quadrado, aonde Y é a área do quadrado e X é o lado do quadrado. Para todos os valores reais de X teremos um valor real de Y. Exemplo: X=2 -> Y=4; X=-3 ->Y=9

Até aqui, tudo bem...

Mas note que existem todo um universo de números em Y que não é mapeado de X. Isso é fácil de ver: todo Y que seja menor que zero não é mapeado, ou de forma mais matemática: não há X real que seja mapeado em Y negativo.

Isso em si não é um problema. O problema é quando queremos resolver a questão inversa: ou seja dado Y, qual é o X que foi mapeado? Note que para Y maior que zero teremos duas soluções reais. Já para  Y menor que zero não teremos nenhuma solução real. Exemplo? Y = 9 <- ou="" x="3</b">. Mas Y=-9 não tem X real que satisfaça ao problema.

Então o problema é que X real é TOTALMENTE mapeado em Y real, mas Y real NÃO É TOTALMENTE mapeado com X real. Como na analogia anterior, estaríamos falando de uma área negativa, então fica fácil supor que o número que fizer este tipo de mapeamento tem de ser imaginário (afinal não existe área negativa).

Mas pode ser que não estejamos falando de áreas. E aí a noção de que o número cujo quadrado é um número negativo não pode existir na realidade fica sem sentido.

O sentido verdadeiro é que o número X que faz Y negativo (ou seja o X que faz Y=-9) não pertence a reta dos reais (que é o que tanto o eixo X quanto o eixo Y são: mapas dos números reais).

Para resolver esta questão do mapeamento, ou seja mapear de forma completa os números reais do eixo Y é que são introduzidos os números complexos que possuem duas partes, ou seja ao invés de preencherem uma reta (ou eixo), eles preenchem um plano.

Em resumo: para garantir que o problema inverso de encontrar as raízes de uma equação tenha sempre solução é necessário incluir no mapeamento não apenas o eixo dos reais, mas todo um plano real (que será descrito por duas coordenadas).

segunda-feira, 30 de setembro de 2013

Probabilidade de Chuva

Hoje fui trabalhar e usei a bicicleta como meio de locomoção. Isto eu já venho fazendo há algum tempo. Mas o que diferiu do normal hoje é que antes de sair eu dei uma olhada na previsão do tempo para o dia.

E lá tinhamos: 16-17 h - 40% de probabilidade de chuva & 17-18 h - 50% de probabilidade de chuva.

Bem, ainda assim eu fui de bicicleta. E não choveu.

No fundo foi sorte. Mas o que isso queria realmente dizer? Bom, 50% de chance é realmente o mesmo que as chances de cara ou coroa. Pode ou não acontecer, e dado que no domingo tivemos uma quantidade razoável de chuva então era meio claro que possivelmente choveria.

Mas há uma lição aqui: probabilidade não é certeza. E o fato de ter chovido no dia anterior razoavelmente e ter a probabilidade de 50% para um determinado horário é basicamente um indicador que as chances de chuva não são tão grandes assim quanto 50% fazem parecer...

quarta-feira, 25 de setembro de 2013

Médias, Máximos e Mínimos

Um ponto que eu mesmo vivo esquecendo é que a média de uma distribuição divide a mesma em duas partes iguais.

Assim, se eu escolher um número de uma distribuição qualquer arbitrariamente, então há uma chance de 50% deste número ser maior que a média e uma chance de 50% deste número ser menor que a média. Bem, estritamente falando esta descrição é a da mediana e não da média - mas vamos considerar as duas iguais por aqui por enquanto.

O leitor pode pensar: "Duh!". Mas há mais informações que podemos obter aqui. Se tirarmos dois números X1 e X2 podemos considerar a média M e dizer:
p(X1 menor que M) = p(X1 maior que M) = 1/2
p(X2 menor que M) = p(X2 maior que M) = 1/2

Para simplificar vamos considerar X1 menor queX2 (pode ser ao contrário - não importa).Ao combinamos as duas alternativas temos algo interessante (se não há correlação entre os eventos):
p(X1 menor que M e X2 menor que M) = p(X1 menor que M)*p(X2 menor que M) = 1/4
p(X1 maior que M e X2 maior que M) = p(X1 maior que M)*p(X2 maior que M) = 1/4
p(X1 menor que M menor que X2) = 1/2

Isto pode ser dito da seguinte forma: chance da média estar entre X1 e X2 é 50% a chance da média estar fora disso é 50%. Mas a chance de ser maior que X2 é 25% e a chance de ser menor que X1 é 25%.

Como isso é possível? Note que p(X1 maior que M) (probabilidade de X1 ser maior que é média) é exatamente a mesma coisa que p(M menor que X1) (probabilidade da média ser menor que X1). E que p(X2 maior que M) (probabilidade de X2 ser maior que a média) é exatamente a mesma coisa que p(M menor que X2) (probabilidade da média ser menor que X2).

Então o que é p(X1 maior que M e X2 maior que M)? É a mesma coisa que p(M menor que X1 e M menor que X2). E pela independência (na realidade, pela não correlação) temos que isso é igual a p(M menor que X1)*p(M menor que X2) = p(X1 maior que M)*p(X2 maior que M) = 1/4 (ou seja a probabilidade que a média seja menor que X1 e menor que X2). Do outro lado temos p(M maior que X1)*p(M maior que X2) = p(X1 menor que M)*p(X2 menor que M) = 1/4 (ou seja a probabilidade que a média seja maior que X e maior que X2).

Assim, o que sobra é a probabilidade que a média esteja entre os dois valores sorteados - que é 50%.

A coisa começa a ficar interessante quando aumentamos a quantidade de números sorteados. Ao chegarmos a 10 temos 511 em 512 chances que a média esteja entre o maior e o menor número sorteados. Isso corresponde a 99.8% de chance que a média esteja entre o maior valor sorteado e o menor valor sorteado. Isso porque estamos calculando a probabilidade que a média não satisfaça uma das duas situações:
- Probabilidade da média ser menor que o menor valor sorteado = probabilidade da média ser menor que todos os valores sorteados.
- Probabilidade da média ser maior que o maior valor sorteado = probabilidade da média ser maior que todos os valores sorteados.

E isto somado a probabilidade que a média esteja entre o menor e o maior valor sorteado cobre todos os eventos possíveis (ou seja dá 1). Portanto para N sorteios:
p(Xmin menor qu eM menor que Xmax) = 1-2/2^N = 1- 1/(2)^(N-1) = 0.98 (para N=7)

E qual é o interesse disso? Bom, podemos fazer coisas mais curiosas se encontrarmos os valores mais altos e mais baixos e fazer alguns cálculos com eles. Por exemplo com relação ao segundo menor valor Xmin2 e o segundo maior valor Xmax2 temos:
p(Xmin menor que M menor que Xmax) = 1 - 2/2^N-2*N/2^N = 1- (N+1)/(2)^(N-1) = 0.88 (para N=7)

E assim por diante. Então pense nisso: ao se obter uma sequência de dados aleatórios da mesma distribuição podemos estimar aonde está a média olhando o máximo e o mínimo - com razoável grau de certeza!

domingo, 22 de setembro de 2013

Sobre a Invasão da Reitoria

Esta semana tivemos uma invasão na reitoria. Felizmente, as partes negociaram uma saída e tudo se encerrou com poucas confusões.

Mas tenho algo a ser dito sobre esta invasão. Primeiro, caro leitor, chamo o evento pelo seu nome correto: invasão. O uso do nome ocupação é um daqueles truques mentais que já falei tanto neste blog: inadvertidamente achamos que ao ocupar-se dá-se uma conotação positiva ao evento.

Ah, mas não tem conotação positiva? Fica difícil ter quando se usa violência para conseguir "ocupar" o espaço Excetuando os casos aonde o Estado de direito resolve mesmo ocupar (e via de regra há um caminho judicial para isso), o uso da violência dificilmente pode ser caracterizado como "violência do bem".

Sim, eu entendo que muitos alunos foram levados a isso pela mentalidade de multidões, que apenas uns poucos agiram de má fé, mas nada disso, absolutamente nada disso muda o fato que a violência foi usada de forma a se atingirem os objetivos da invasão. E a resposta moralmente correta é que os invasores assumam a responsabilidade pela violência ocorrida.

Isso em um mundo perfeito... Já no mundo real, os envolvidos continuaram a dizer que não se arrependem de nada, que fizeram mas outros fazem pior (olha o Tu quoque aí, minha gente!), mas irão fugir da responsabilidade como o diabo foge da cruz.

Mas em defesa deles, posso argumentar que caráter é algo muito complicado de desenvolver e não é realista esperar que todo mundo o tenha (talvez pedaços, não sei...). Bom, não vale a pena entrar em um monólogo de uma nota só. Vou apenas dizer um fato muito óbvio: a razão pela qual concordamos em regras, é que elas evitam arbitrariedades - e, como o leitor deve bem saber, arbitrariedades podem ser cometidas por qualquer lado ou pessoa com recursos suficientes.

E não acredite nas arbitrariedades do bem: elas vem do mesmo lugar que o Coelhinho da Páscoa.

Por fim deixo aqui uma entrevista bastante lúcida do reitor.

quarta-feira, 11 de setembro de 2013

Enrole-se quem Quiser

Nestes tempos de internet é muito fácil cair em esparrelas... Mas ao mesmo tempo fica mais fácil detectar várias enrolações que correm por aí.
Governo, Médicos... Todo mundo enrolando...
Uma em particular que vários já caíram é sobre o programa Mais Médicos ser uma resposta as Manifestações ocorridas em Junho.
Well... Não é tanto fugir, mas apenas não comparecer...
Não são. A idéia de trazer médicos estrangeiros para o Brasil é ANTERIOR às manifestações...
Os médicos brasileiros não entenderam o poder das imagens...
As idéias por trás do Mais Médicos são de maio de 2013 (na realidade até mais antigas...): "As negociações para o envio dos médicos cubanos para o Brasil foi iniciada pela presidente Dilma Rousseff, em janeiro de 2012, quando visitou Havana, a capital cubana. Ela defendeu uma iniciativa conjunta para a produção de medicamentos e mencionou a ampliação do envio de médicos cubanos ao Brasil, para apoiar o atendimento no SUS (Serviço Único de Saúde)."
Pouco importa se isto é verdade ou não no curto prazo: a imagem é o que conta! Aprendi isso da pior forma.
Então como foi que se vendeu a idéia que o Mais Médicos é uma resposta às ruas? Ora caro leitor, isso é muito simples... Há dois mecanismos em ação aqui:
  1. As reivindicações das manifestações de junho foram difusas (menos corrupção, mais saúde, etc... - eu já discuti isso neste blog).
  2. A classe política se valeu de um velho truque de entrevista: o entrevistador faz uma pergunta, mas o entrevistado responde outra coisa - que é do seu interesse divulgar.
Junte isso com um marketing razoável, um jornalismo de segunda categoria (não questiona, apenas reporta como um desmemoriado), e um povo meio ... humm... como direi isso sem ofender?... Hummm... bem brasileiro então.
Assim como Samba e a Feijoada (hummm, será?).
E aí temos a noção sem noção que o Mais Médicos é resposta do governo às ruas, quando na realidade é um plano já bem adiantado.

terça-feira, 10 de setembro de 2013

Previsões de Inflação - agora com os dados de agosto

Mais uma vez chegamos aquele dia do mês: as previsões de Inflação de 2013 - Muito provavelmente entre 5% e 6.7% (com 98% de probabilidade de ser abaixo de 6.7%).

Comparando os resultados só com janeiro (6.32%), incluindo até fevereiro (6.44%), incluindo até março (6.43%), incluindo até abril (6.51%). incluindo até maio (6.39%), incluindo até junho (6.16%), incluindo até julho (5.68%), temos mais uma melhora em agosto (5.43%) - mesmo com o IPCA de agosto um pouco mais alto (0.24%).

Naturalmente o valor esperado não é tudo. Muito mais útil são os intervalos de confiança. No caso atual, com os dados de julho temos uma chance de 0.06% da inflação ficar acima de 7.56% (a chance que a inflação fique acima de 6.7% é de 1.4%). Podemos afirmar que há uma chance de 99.94% da inflação de 2013 ficar abaixo dos 7.58%.

O limite inferior é 4.87% (compare com janeiro - 4.77%, e os acumulados até: fevereiro (5.04%), março (5.17%), abril (5.38%), maio (5.4%), junho (5.31%) e julho (4.98%)). Já o superior teve um progressão similar: 14.96%, 14.28%, 13.46%, 12.73%, 11.82%, 10.78%, 9.5%, e finalmente 8.5% (incluindo agosto). Sendo que valores superiores a 7.56% tem probabilidade inferior de ocorrência 0.006%.

As notícias são boas para economia.Claro que não sou uma firma que faz previsões de inflação, mas se eu tiver que estimar (chute educado) diria o seguinte:

  • 51% de chance da inflação abaixo de 4.9%
  • 88% de chance da inflação abaixo de 5.8%
  • 98% de chance da inflação abaixo de 6.7%
  • 99% de chance da inflação abaixo de 7.6%

E vejamos como o mês que vem...

quinta-feira, 5 de setembro de 2013

Mancadas...

Há poucas coisas tão embaraçosas para um engenheiro ou arquiteto quanto dar uma mancada no projeto.
Um exemplo em particular me vem a mente...
Mas existem mancadas e Mancadas... E a que vem a seguir é com M maiúsculo mesmo. Um prédio em Londres está literalmente cozinhando o que está a sua frente (ou um pouco mais abaixo). O efeito "lente de aumento" é causado pela concavidade do prédio.
Concavidade + Estruturas Espelhadas Gigantes = Raio da Morte?
Essencialmente: o sol manda uma determinada densidade de potência S (W/m2). A superfície refletora do prédio tem uma área A. Isso quer dizer que o prédio refletirá uma potência proporcional a S*A.

Até aí tudo aí tudo bem. Em prédios com bastante alvenaria e poucos vidros, esta proporção é baixa, ou em outras palavras a eficiência de reflexão é baixa. Mas já em prédios com muito vidro, a coisa é na realidade um espelho gigante.

Mesmo assim, tirando a reflexão de um espelho gigante adicional o problema ainda não existe. Por que? Porque a densidade de potência refletida é praticamente a mesma. Então é simplesmente o caso de reflexões do sol em locais inesperados devido a um "efeito espelho".
Física... Sua danadinha!
Mas a coisa muda de figura quando há concavidade. Há um caso específico em que há magnificação. Neste caso a potência S*A fica concentrada em uma área substancialmente menor do que a do edifício. E aí potência elevada (o edifício tem 37 andares - em torno de 120 metros - por cerca de 60 metros. Considerando 250 W/m2 temos uma potência de 1.8 MW. Um forno de cozinha elétrico tem potência de cerca de 1.3 kW).
Não é nem difícil de fazer...
Então você tem um edifício e um forno solar.
Yes it can be a death ray.
A propósito, não culpe os arquitetos ingleses: o design é de outra pessoa.
Como o leitor pode ver, não é a primeira vez que este arquiteto faz isso.

segunda-feira, 2 de setembro de 2013

A Pior Solução

A Pior Solução é aquela que demonstra de um modo claro e inequívoco que não se sabe direito sequer qual é o problema.

Pois foi essa a posição do governo brasileiro frente a questão Snowden: o governo encomendou aos correios um sistema de e-mail nacional.
Com certeza tem alguma coisa a ver com teletransporte!!!
Isto mostra uma série de problemas, sendo, provavelmente, o maior deles é o total e completo desconhecimento do que seja um sistema de e-mail. Duvida, caro leitor? Então leia isto aqui.

Pode ser que o leitor possa estar perguntando: "Ora RAX, você está se precipitando. Afinal, pode ser apenas que os jornalistas não tenham entendido a explicação do governo. Afinal a "mídia" é mestre em distorcer comentários de todos".

Bem isto seria uma linha de raciocínio lógico não fosse um pequeno obstáculo: pedir aos correios para desenvolver um sistema de e-mail nacional? Aos correios? Sistema de e-mail nacional? Talvez se fosse para ABIN, para o Ministério da Defesa, ou mesmo para o tal gabinete institucional da presidência. Eu posso até ver uma possível linha de raciocínio: e-mail é correio eletrônico, então quem melhor do que os correios para resolver esta questão?

Surely, you must be joking, dear reader!